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Abstract 

The  synthesis  and  physical  properties  of  cyanurate  networks  formed  from  two  new  tricyanate 
monomers,  l,3,5-tris[(4-cyanatophenylmethyl]benzene,  and  3,5-bis[(4- 

cyanatophenylmethyl)]phenylcyanate,  are  reported  and  compared  to  those  of  and  1,1,1  -tris[(4- 
cyanatophenyl)] ethane  (also  known  as  ESR-255).  All  three  networks  possessed  somewhat 
different  aromatic  contents  and  cross-link  densities,  however  the  thermochemical  stability  of 
these  networks,  as  determined  by  TGA,  was  outstanding,  with  that  of  l,3,5-tris[(4- 
cyanatophenylmethyl)]benzene  being  among  the  best  known  for  organic  cyanate  esters  despite 
its  comparatively  high  segmental  flexibility.  Moreover,  the  moisture  uptake  of  cured  1,3,5- 
tris[(4-cyanatophenylmethyl)]benzene,  at  2.2%  after  96  hours  immersed  in  85  °C  water,  was 
comparatively  low  for  a  cyanate  ester  network  with  a  glass  transition  temperature  of  320  °C  at 
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full  cure.  When  cured  for  24  hours  at  210  °C,  the  dry  glass  transition  temperatures  of  the 
networks  ranged  from  245  to  285  °C,  while  the  wet  glass  transition  temperatures  ranged  from 
225  -  240  °C.  The  similarity  in  glass  transition  temperatures  resulted  from  a  lower  extent  of 
cure  in  the  networks  with  more  rigid  segments.  In  essence,  for  networks  with  very  high  glass 
transition  temperatures  at  full  cure,  the  process  conditions,  rather  than  the  rigidity  of  the  network, 
determined  the  attainable  glass  transition  temperature.  Because  networks  with  a  higher  extent  of 
cure  tend  to  exhibit  slower  long-term  degradation,  in  this  case,  the  networks  with  greater 
segment  flexibility  enabled  superior  performance  despite  exhibiting  a  lower  glass  transition 
temperature  at  full  cure.  These  results  illustrate  that,  in  contrast  to  the  prevailing  heuristics  for 
improving  the  performance  of  high-temperature  thermosetting  polymer  networks,  a  more  flexible 
network  with  a  lower  glass  transition  temperature  at  full  cure  can  offer  an  optimal  combination 
of  thermomechanical  and  thermochemical  performance. 

Keywords:  thermosetting  resin;  cyanate  ester;  cure  kinetics 
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Introduction 


Tightly  cross-linked  macromolecular  networks  with  maximum  continuous  use 
temperatures  above  175  °C,  formed  by  thermal  cure  of  polyfunctional  monomers,  represent  an 
extremely  important  class  of  polymeric  structural  and  adhesive  materials,  with  a  broad  range  of 
applications  in  the  electronics,  energy  storage,  and  aerospace  industries.  Among  these,  cyanate 
ester  resins1'3  (R-0-C=N)  made  from  trifunctional  monomers  have  been  studied4,5  as  potential 
materials  for  applications  ranging  from  aircraft  engine  turbine  brush  seals,6'8  to  nano-structured 
bio-sensors,9  due  to  their  excellent  high-temperature  mechanical  performance  combined  with 
ease  of  processing.  However,  cyanate  ester  monomers  also  possess  unusual  characteristics  that 
make  them  ideally  suited  to  the  study  of  high-performance  thermosetting  networks,  such  as  the 
highly  selective  formation  of  a  single  cross-linked  network  product10  (a  polycyanurate  network), 
no  need  to  balance  stoichiometry,  effective  control  of  polymerization  rate  through  a  wide  range 
of  available  catalysts,  low  shrinkage  on  cure  that  minimizes  residual  stresses,  ’  outstanding 
short-term  thermal  stability  of  the  network,  and  simple  verification  of  the  chemical  structure  of 
the  insoluble  cured  network  (including  quantitative  information  on  side  product  formation) 
through  FT-IR  spectroscopy.  As  a  result,  studies  of  the  relationships  among  structure, 
processing,  and  properties  involving  polycyanurate  networks  can  provide  significant  insights  into 
the  properties  of  thermosetting  networks  in  general. 

When  developing  new  high-performance  thermosetting  monomers,  it  has  generally  been 
assumed  that  the  glass  transition  temperature  ( Tg )  of  the  fully  cured  network  (herein  referred  to 
as  a  “fully-cured  Tg')  should  be  as  high  as  possible,  provided  that  an  acceptable  degree  of 
toughness  can  be  maintained.  For  single-components  networks  with  one  well-defined  type  of 
chemical  cross-link,  the  fully-cured  Tg  depends  on  the  density  of  cross-links  and  branch  points, 
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along  with  the  rigidity  of  the  network  segments,  but  not  on  the  details  of  the  cure  process.  As  a 
result,  the  “fully  cured”  Tg  provides  a  useful  quantitative  (though  limited)  indication  of  the 
relative  flexibility  of  network  segments.  For  cyanate  esters  such  as  Primaset®  PT-30  with  fully 
cured  Tg  s  as  high  as  400  °C,14  however,  complete  conversion  of  the  monomer  to  cyanurate  may 
require  undesirably  high  temperatures.15  Therefore,  the  network  properties  depend  on  the  “as 
cured”  rather  than  the  fully-cured  Tg.  In  many  thermosetting  materials,  vitrification  (the  point 
when  network  development  drives  the  glass  transition  temperature  above  the  cure  temperature) 
has  been  assumed  to  result  in  cessation  of  cure,  resulting  in  an  “as  cured”  Tg  that  is  roughly  the 
same  as  the  maximum  cure  temperature.  As  a  result  of  these  assumptions,  the  long-established 
approach  to  developing  high  Tg  thermosetting  resins  has  been  to  1)  make  the  network  segments 
as  rigid,  and  cross-link  densities  as  high,  as  constraints  imposed  by  brittleness  permit,  and  2) 
utilize  the  maximum  permissible  cure  temperature.  Even  though  cyanate  esters  have  been 
known  for  many  years  to  violate  the  assumptions  relating  the  Tg  and  maximum  cure 
temperature,16  the  approach  to  developing  cyanate  esters  has  been  the  same  as  the  one  outlined 
above  for  high-temperature  thermosetting  networks  in  general. 

Recently,  however,  it  has  become  more  widely  appreciated  that  in  thermosetting 
networks,  slow  conversion  of  monomer  continues  to  take  place  in  the  glassy  state,  and  thus  that 
“as  cured”  Tg  s  may  significantly  exceed  maximum  cure  temperatures.  The  kinetics  of  cure  for 
both  epoxy  and  cyanate  ester  resins  near  Tg  have  been  studied  recently.18,19  However,  many 
aspects  of  cure  well  below  Tg  remain  relatively  unexplored.  From  a  performance  standpoint, 
though,  the  highly  desirable  combination  of  good  thermomechanical  performance  (which 
demands  a  high  Tg )  and  affordable  processing  (which  requires  low  cure  temperatures)  means  that 
a  deeper  understanding  of  cure  well  below  Tg  will  become  increasingly  important. 
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Well  below  Tg,  large-scale  molecular  motion  in  the  developing  network  is  generally 
absent.  As  illustrated  in  Figure  1,  the  ability  for  cure  to  proceed  in  a  system  containing  reactive 
groups  only  at  chain  ends  is  likely  to  depend  strongly  on  the  degree  of  localized  motion  available 
to  these  chain  ends.  The  flexibility  of  the  individual  network  segments  thus  plays  a  key  role  in 
either  facilitating  or  hindering  the  formation  of  networks  with  very  high  levels  of  conversion. 
The  incorporation  of  flexible  chemical  linkages  that  enable  such  motions,  however,  tends  to 
result  in  a  lower  fully-cured  Tg.  Due  to  the  historical  considerations  outlined  earlier,  the  concept 
of  incorporating  more  flexible  network  segments  to  facilitate  cure  well  below  Tg  has  thus  far 
received  little  consideration  in  the  development  of  high-temperature  thermosetting  polymers. 
The  trade-offs,  however,  between  achieving  full  cure  and  maintaining  a  high  Tg  are 
technologically  quite  important  because  the  long-term  hydrolytic  stability  of  cyanate  esters,  for 
instance,  depends  critically  on  achieving  close  to  100%  conversion,  while  the  toughness  of 
cyanate  esters  has  been  shown  to  more  than  double  when  conversion  is  increased  from  95%  to 
near  100%.  Despite  their  importance,  these  trade-offs  have  not  been  explored  systematically  in 
many  types  of  thermosetting  resins. 

We  recently  compared  the  performance  of  the  commercial  polycyanurate  Primaset®  PT- 
30,  which  forms  a  highly  rigid  network,  with  a  more  flexible  analog,  and  found  that  many, 
though  not  all,  of  the  expected  disadvantages  associated  with  more  flexible  network  segments 
were  mitigated  by  the  relative  ease  of  achieving  close  to  complete  conversion  in  the  more 
flexible  network.  However,  we  also  found  that  the  thermochemical  stability  of  the  more 
flexible  network  was  significantly  reduced.  In  addition  to  degrading  application-related 
performance,  the  reduction  in  thermal  stability  prevented  a  clear  determination  of  the  fully  cured 
Tg  of  the  more  flexible  network.  A  thermochemically  more  stable  network  with  enhanced 
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segment  flexibility  was  thus  needed  to  systematically  understand  the  trade-offs  between  segment 
flexibility,  “as  cured”  and  fully  cured  Ts,  and  performance. 

a) 


b) 


Figure  1.  Schematic  illustration  of  the  envisioned  effect  of  network  segment  flexibility  on  the 
reactivity  of  cyanate  esters  in  the  glassy  state.  After  vitrification,  portions  of  the  network  away 
from  the  chain  ends  (indicated  by  the  darker  ellipses)  exhibit  highly  restricted  motion,  and  the 
ability  of  the  reactive  groups  to  move  into  the  correct  position  and  orientation  for  further  cure 
depends  on  the  local  motions  of  the  nearby  chain  segment,  as  illustrated  by  the  arrows,  a)  a 
portion  of  the  network  formed  from  the  rigid  monomer  ESR-255,  having  a  limited  range  of 
motion  b)  a  portion  of  the  network  formed  from  the  more  flexible  monomer  1  (see  Scheme  1) 
that  exhibits  a  greater  range  of  plausible  motion  due  to  a  larger  number  of  more  freely  rotating 
bonds. 
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As  reported  in  this  paper,  a  newly  synthesized  cyanate  ester  monomer  (1)  affords 
enhanced  network  segment  flexibility  with  no  loss  in  thermochemical  stability  compared  to  more 
rigid  analogs  (see  Chart  1  for  chemical  structures).  The  synthesis  of  this  monomer,  along  with 
the  synthesis  of  more  rigid  single-component  analogs,  thus  enabled  a  clear  comparison  of  the 
effect  of  network  segment  flexibility  on  both  extent  of  cure  and  Tg.  The  comparisons  indicated 
that  the  “as  cured”  Tg  (using  identical  cure  conditions)  was,  though  far  above  the  cure 
temperature  itself,  still  limited  to  a  maximum  value  determined  primarily  by  the  cure  conditions 
rather  than  by  the  rigidity  of  the  network  segments.  In  contrast,  the  fully  cured  Tg  of  the 
networks  clearly  decreased  with  increasing  segment  flexibility.  Therefore,  under  a  given  set  of 
cure  conditions,  near  complete  conversion  was  only  possible  in  networks  having  sufficient 
segment  flexibility.  Because  near  complete  conversion  is  often  a  requirement  for  acceptable 
long-term  performance,  in  this  case  the  more  flexible  network,  contrary  to  expectations  based  on 
the  prevailing  heuristics  described  previously,  may  provide  the  best  performance.  Therefore,  in 
addition  to  developing  a  new  cyanate  ester  monomer  with  superior  performance,  new  insights 
have  been  gained  into  the  role  of  network  segment  flexibility  and  its  usefulness  in  optimizing 
performance.  These  insights  will  greatly  aid  both  process  and  product  development  for  a  wide 
range  of  thermosetting  polymers  used  in  many  high-performance  electronics,  aerospace,  and 
energy  applications. 
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Chart  1.  Chemical  structures  of  new  tricyanate  esters  (1)  and  (2),  and  more  rigid  analogs  ESR- 
255  and  PT-30. 


PT-30 


Results  and  discussion 

Synthesis  of  the  new  monomers  1  and  2.  The  synthesis  of  tricyanate  1  is  shown  in 
Scheme  1.  The  preparation  of  1,3,5-triaroylbenzenes  such  as  5  has  been  typically  accomplished 
via  the  enaminoketone  generated  from  acetophenones  and  expensive  dialkyl  acetals  of 
dimethylformamide.23’24  In  this  case,  /?ara-methoxyacetophenone  was  condensed  with  methyl 
formate  using  potassium  fert-butoxide  as  base  to  give  the  potassium  salt  4  in  good  yield.25'27 
This  salt  was  then  cyclotrimerized  in  hot  glacial  acetic  acid  to  give  5  in  good  yield.28  Reduction 
of  the  carbonyls  by  sodium  borohydride  gave  the  diastereomeric  mixture  of  triol  6. 
Hydrogenolysis  of  6  by  catalytic  transfer  hydrogenation29"32  gave  the  unreported  symmetrical 
tribenzylbenzene  7.  Demethylation  of  7  by  molten  pyridine  hydrochloride  gave  the  triphenol  8.33 
Single-crystal  X-ray  analysis  proved  the  structure  of  8  (see  Supporting  Information  Section  SI). 
Finally,  8  was  cyanated  with  cyanogen  bromide  and  triethylamine  base  to  give  tricyanate  1  at 
30%  overall  yield  in  the  six  steps  from  /?ara-methoxyacetophenone.  NMR  spectra  for 
monomers  1,  2,  and  ESR-255  are  shown  in  Supporting  Information  Section  SI. 
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Scheme  1.  Chemical  synthesis  route  to  tricyanate  1. 


NH4HCO2,  Pd/C,  HOAc/^- 


91% 

pyridineHCI,  180  C 
90% 


6  R-i  =  OH,  R2  =  Me 

7  R,  =  H,  R2  =  Me 

8  R1f  R2  =  H 


The  synthesis  of  tricyanate  2  is  shown  in  Scheme  2.  5-Hydroxyisophthalic  acid  was 
heated  with  dimethylsulfate  and  potassium  carbonate  in  acetone  to  give  the  trimethyl  derivative 
9. 34  Hydrolysis  by  refluxing  in  aqueous  potassium  hydroxide  gave  the  diacid  10.35  Reaction 
with  thionyl  chloride  followed  by  Friedel-Crafts  acylation  of  anisole  gave  the  unreported 
diketone  11  in  52%  yield  over  the  two  steps.  Reduction  with  lithium  tetrahydroaluminate  in 
tetrahydrofuran  gave  the  diastereomeric  mixture  of  diol  12.  Hydrogenolysis  of  the  benzylic 
alcohols  in  refluxing  acetic  acid  with  palladium  catalyst  gave  the  dibenzyl  compound  13.36 
Deprotection  followed  by  cyanation  in  similar  fashion  to  that  described  above  gave  tricyanate  2 
in  8%  yield  over  the  seven  steps. 
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Scheme  2.  Chemical  synthesis  route  tricyanate  2. 


S02(M6)2 

k2co3 

acetone 


94% 


KOH 

MeOH/H20 

- ► 


95%  MeO 


C02H 


COoH 


52%  1 .  SOCI2 

2-steps  w  2.  MeOPh,  AICI3,  CS2 


In  addition  to  the  new  monomers  1  and  2,  the  rigid  cyanate  ester  known  as  ESR-255, 
originally  reported  by  Shimp  et  al ,,37  was  also  prepared.  These  three  compounds  form  a 
convenient  series  for  comparison  purposes:  all  are  composed  of  three  cyanated  phenyl  rings; 
ESR-255  has  no  flexible  methylene  linkages  between  rings,  2  has  two  flexible  methylene 
linkages  and  is  a  small-molecule  analog  of  the  oligomeric  commercial  Primaset®  PT-30  resin, 
while  1  has  three  methylene  linkages  and  thus  features  the  greatest  degree  of  network  segment 
flexibility. 

Thermomechanical  analyses.  By  avoiding  aliphatic  network  junctions,  themochemical 
stability  comparable  to  Primaset®  PT-30  was  maintained  in  the  more  flexible  network  formed 
from  cured  1,  as  illustrated  by  the  comparative  TGA  data  in  Figure  2.  Note  that  2  shows  similar 
thermochemical  degradation  behavior  to  that  previously  reported  for  PT-30,22’38  as  expected. 
The  char  yields  of  over  70%  at  600  °C  in  both  nitrogen  and  air  observed  in  cured  1  are  among 
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the  best  known  for  cyanate  esters,  though  many  other  recent  approaches39'41  may  ultimately 
provide  better  heat  and  flame  resistance  overall,  and  more  detailed  long-term  testing  will  be 
needed  to  determine  the  applicability  of  these  results  to,  for  example,  the  thermo-oxidative 
stability  of  composite  structures.  Despite  having  the  highest  aliphatic  content  and  largest 
molecular  weight  between  network  junctions  (see  Table  1)  of  all  of  the  materials  studied,  the 
network  formed  by  curing  1  was  the  most  chemically  stable  at  elevated  temperatures. 
Comparing  the  data  for  cured  1  and  cured  2  reveals  that  the  increased  phenyl  content  relative  to 
cyanurate  did  not  improve  the  char  yield. 


100%  - 

g>  80% 

c 

|  60%  - 

Cured  1 

a> 

*  40% 

Cured  2 

(/> 

03 

Cured  PT-30 

S  20%  - 

0%  - 

- 1 - 1 - 1 - 1 - 

350  400  450  500  550  600 

Temperature  (°C) 

Figure  2a.  TGA  of  cured  cyanate  ester  networks  featuring  only  phenyl  rings  and  methylene 
bridges,  under  nitrogen.  Note  that  data  for  PT-30  has  previously  been  published  (reference  22). 
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Figure  2b.  TGA  of  cured  cyanate  ester  networks  featuring  only  phenyl  rings  and  methylene 
bridges,  in  air. 


Table  1 

Chemical  Characteristics  of  Cyanurate  Networks  Relevant  to  TGA  Results 


Monomer 

Weight  Fraction  of  Pure  Cyanurate  Network 
Classified  as: 

Average 

Molecular  Weight 
Between  Network 
Junctions  (g/mol)d 

Char  Yield 
in  Air  at  600 
°C  (wt%) 

Phenyl3 

Aliphatic15 

Cyanurate0 

1 

0.65 

0.09 

0.27 

157 

74 

2 

0.60 

0.08 

0.33 

127 

73 

ESR-255 

0.60 

0.07 

0.33 

122 

59 

a.  Includes  phenyl  rings  (not  cyanurate  rings)  and  their  associated  hydrogens 

b.  Includes  aliphatic  carbons  and  their  associated  hydrogens 

c.  Includes  cyanurate  rings  and  their  neighboring  oxygen  bridges 

d.  Note  that  branch  points  in  the  monomer  count  as  network  junctions,  and  short  side-chains 
(e.g.  terminal  methyl  groups  not  along  a  path  between  junctions)  are  not  included 


Figure  3  shows  the  TGA  data  for  cured  ESR-255  compared  to  1  and  2.  Despite  having  as 
low  or  lower  aliphatic  content  than  either  1  or  2,  the  weight  loss  at  high  temperatures  in  ESR-255 
was  significantly  higher  than  that  of  both  1  and  2  under  both  nitrogen  and  air.  As  with 
previously  studied  flexible  tricyanate  esters,22  ESR-255  features  an  aliphatic  network  junction 
comprised  of  an  aliphatic  carbon  next  to  a  methyl  group.  Based  solely  on  empirical  observation 
of  both  these  and  previously  studied  cyanate  esters  without  considering  mechanistic  details,  we 
surmise  that  while  some  types  of  junctions  formed  by  tertiary  or  quaternary  aliphatic  carbons 
may  adversely  affect  thermochemical  stability,  the  presence  of  methylene  bridges  between 
phenyl  groups  does  not  significantly  detract  from  the  thermochemical  stability  of  highly  aromatic 
cyanate  esters.  In  the  specific  case  of  ESR-255,  the  presence  of  unreacted  cyanate  esters  groups 
at  high  temperatures  (as  detailed  in  later  sections)  may  also  negatively  affect  thermochemical 
stability.  The  aforementioned  conclusions  are  in  general  agreement  with  the  results  of  previous 
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comparative  studies  of  thermochemical  stability  in  cyanate  esters.1,42  The  outstanding 
thermochemical  stability  of  cured  1  also  means  that  any  significant  loss  in  mechanical  stiffness 
below  400  °C  in  the  fully  cured  network  may  be  attributed  definitively  to  a  comparatively  low 
value  of  the  glass  transition  temperature  brought  about  by  more  flexible  network  segments, 
rather  than,  as  was  the  case  with  previously  synthesized  flexible  cyanurate  networks,  some 
combination  of  mechanical  and  chemical  effects. 


Figure  3a.  Comparative  TGA  of  all  samples  synthesized  for  this  work,  under  nitrogen. 


Figure  3b.  Comparative  TGA  of  all  samples  synthesized  for  this  work,  in  air. 


Cure  Kinetics  and  Tg.  The  impact  of  selectively  incorporating  methylene  groups  on 
flexibility  and  the  subsequent  impact  on  ease  of  cure  in  cyanate  esters  can  be  seen  by 
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examination  of  non-isothermal  DSC  data.  Figure  4  compares  the  DSC  exotherm  of  monomers  1, 
2,  and  ESR-255  calculated  using  a  re-scanned  baseline  (that  is,  by  subtracting  the  heat  flow  seen 
on  the  second  heating  from  that  seen  on  the  first  heating).  Whereas  1  exhibited  a  symmetrical 
curve  with  an  enthalpy  of  cyclotrimerization  of  around  120  kJ/cyanate  eq.  (slightly  higher  than 
expected),  2  exhibited  a  more  asymmetric  curve  with  a  close  to  expected  enthalpy  of  106  kJ/eq. 
In  contrast,  ESR-255  showed  a  highly  asymmetrical  curve  with  an  “open”  tail.  Previous  work  on 
cyanate  ester  resins,16  as  well  as  our  own  verification  using  modulated  DSC  (see  Supporting 
Information  Section  S2)  indicated  that  the  open-ended  nature  of  the  exothermic  peak  was  not  due 
to  a  baseline  shift,  but  rather  resulted  from  continued  slow  cure  of  the  sample.  Although  the 
apparent  enthalpy  of  cure  of  ESR-255  based  only  on  the  non-isothermal  DSC  curve  was  a  lower 
than  expected  93  kJ/eq,  combined  FT-IR  and  isothermal  DSC  analysis  (see  Supporting 
Information  Section  S3),  determined  that  the  actual  enthalpy  of  cure  for  ESR-255  was,  as 
expected,  1 10  ±  7  kJ/eq.  Hence  ESR-255  attained  a  conversion  of  only  around  85%,  even  when 
heated  to  350  °C. 


Figure  4.  Comparative  non-isothermal  DSC  at  10  °C  /min.  of  tricyanate  monomers,  showing 
exothermic  region. 
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Insight  into  the  differing  non-isothermal  DSC  behavior  of  these  three  monomers  can  be 
gained  by  considering  the  evolution  of  the  sample  Tg  during  the  DSC  scan,  as  illustrated  in 
Figure  5.  The  Tg  of  the  sample,  which  changes  due  to  cure  during  the  DSC  scan,  is  estimated 
with  the  aid  of  the  diBenedetto  equation  from  conversion  data  computed  via  a  running 
integration  of  the  heat  flow.  Figures  5a-5c  show  the  sample  Tg  superimposed  on  the  scan 
temperature  and  the  heat  flow,  all  as  a  function  of  time  elapsed  from  the  start  of  heating  at  50  °C. 
The  diBenedetto  equation  previously  has  been  used  quite  successfully  to  predict  the  Tg  of  partly 
cured  cyanate  esters,43’44  most  likely  due  to  their  well-defined  cure  chemistry.  The  same  well- 
defined  cure  chemistry  of  cyanate  esters  enables  reliable  computation  of  the  conversion  from  the 
heat  flow  data.1  Details  of  the  determination  of  the  diBenedetto  equation  parameters  for 
monomers  1,  2,  and  ESR-255  are  provided  in  Supporting  Information  Section  S4. 

In  Figure  5a,  it  can  be  seen  that  the  Tg  rapidly  approaches,  but  never  actually  reaches,  the 
scan  temperature  for  monomer  1.  For  monomer  2,  as  seen  in  Figure  5b,  the  Tg  surpasses  the  scan 
temperature  (that  is,  the  sample  vitrifies)  but  then  the  Tg  stabilizes  as  conversion  nears 
completion,  remaining  only  modestly  higher  than  the  scan  temperature.  For  ESR-255,  as  seen  in 
Figure  5c,  even  after  vitrification,  the  cure  is  still  far  from  complete.  Additional  cure  drives  the 
Tg  upwards  quickly  enough  that  the  sample  is  pushed  deeper  into  the  glassy  state.  The 
characteristic  open-ended  “L”-shaped  DSC  trace  thus  appears  to  signify  a  sample  that  has 
vitrified  well  short  of  complete  conversion.  Figures  5b  and  5c  also  show  that,  as  has  been  well- 
documented  for  cyanate  esters,  cure  below  Tg  can  take  place  at  a  slow  but  significant  rate.  The 
particular  shape  of  the  Tg  curve  in  Figure  5c  after  vitrification  suggests  that  the  “as-cured”  Tg 
should  be  approximately  a  linear  function  of  the  cure  temperature  (with  a  slope  greater  than  one) 
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for  samples  that  have  vitrified  but  are  well  short  of  full  conversion,  as  we  recently  observed  for 


1  ^ 

dicyanate  esters. 


Figure  5a.  Comparison  of  sample  temperature  and  estimate  of  Tg  (with  heat  flow  also  plotted 
for  reference)  during  the  DSC  scan  of  1. 


Time  (min.) 


Figure  5b.  Comparison  of  sample  temperature  and  estimate  of  Tg  (with  heat  flow  also  plotted 
for  reference)  during  the  DSC  scan  of  2.  Note  that  the  error  in  the  estimated  Tg  is  comparatively 
large  for  2  due  to  greater  uncertainty  in  the  diBenedetto  parameters. 
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Figure  5c.  Comparison  of  sample  temperature  and  estimate  of  Tg  (with  heat  flow  also  plotted 
for  reference)  during  the  DSC  scan  of  ESR-255. 

The  behavior  of  ESR-255  seen  in  Figure  5c  has  significant  implications  for  the  study  of 
high-temperature  thermosetting  polymers.  For  one,  it  suggests  that,  if  following  the  standard 
practice  of  using  a  “residual  heating”  to  determine  the  extent  of  cure  by  DSC,  precautions  must 
be  taken  to  ensure  that  full  cure  of  the  sample  is  in  fact  achieved.  If  not,  then  a  technique  such  as 
the  combined  DSC/FT-IR  approach  we  describe  in  Supporting  Information  Section  S3  must  be 
utilized  if  true  conversions  are  to  be  measured.  As  the  ESR-255  data  illustrates,  the  errors  in 
conversion  caused  by  ignoring  incomplete  cure  may  be  substantial.  The  significant  amount  of 
cure  below  Tg  implies  that  in-situ  cure  of  vitrified  cyanate  esters  can  occur  during  standard 
characterization  procedures;  such  in-situ  cure  can  lead  to  incorrect  Tg  measurements,  and  may 
even  cause  the  spurious  appearance  of  two  separate  Tg  s  when  in  fact  only  one  is  present.45  Also, 
such  behavior  can  cause  an  “as-cured”  Tg  values  to  move  outside  the  range  of  measurement 
during  a  standard  DSC  scan,  leading  to  the  incorrect  conclusion  that  because  no  glass  transition 
was  observed,  the  “as-cured”  Tg  must  lie  outside  the  scan  range.  Such  phenomena  are  important 
to  consider  when  investigating  high  Tg  thermosetting  polymers  such  as  cyanate  esters. 
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Due  to  the  possibility  of  in-situ  cure,  dynamic  TMA  measurements  must  be  carried  out  at 
a  very  rapid  heating  rate  to  determine  the  “as-cured”  Tg  for  dry  samples  (shown  in  Table  2,  with 
plots  and  more  details  provided  in  Supporting  Information  Section  S5).  The  Tg  values  were 
similar  after  24  hours  of  cure  at  210  °C  under  a  nitrogen  atmosphere  with  no  added  catalyst,  with 
the  value  for  the  more  rigid  ESR-255  being  slightly  lower,  perhaps  due  to  its  markedly  slower 
cure  kinetics.  These  results  contrast  with  the  “fully  cured”  Tg  (determined  by  re-heating  samples 
previously  cured  to  350  °C),  which,  as  also  seen  in  Table  2,  varied  as  expected  with  the  rigidity 
of  the  network  segments  and  cross-link  density.  Using  the  diBenedetto  equation  to  convert  the 
“as  cured”  Tg  values  to  “as  cured”  conversions  confirms  the  trend  seen  in  the  non-isothermal 
DSC  data,  namely,  monomer  1  cured  the  most  readily  of  the  three  monomers  examined,  while 
ESR-255  cured  least  readily.  These  results  have  important  technological  implications. 
Specifically,  when  cyanurate  networks  are  not  fully  cured,  the  residual  cyanate  ester  groups  are 
prone  to  reactions  with  water  at  elevated  temperature  that  eventually  lead  to  formation  of 
volatiles,  which  in  turn  can  lead  to  blistering  and  loss  of  structural  integrity.20,46  Hence,  the  more 
rigid  ESR-255,  although  exhibiting  a  similar  “as  cured”  T*  provides  in  some  aspects  a  lower 
level  of  performance  compared  to  monomer  1  due  to  a  lower  level  of  conversion  (approx.  72% 
for  the  former  versus  about  96%  for  the  latter). 
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Table  2 


Results  of  Thermomechanical  Analysis  on  Cured  Cyanate  Esters 


Monomer 

“As  Cured”  Tg  (°C)a 

“Fully  Cured”  T„ 
(°C)b 

“As  Cured” 
Conversion11 

1 

280  ±  15 

338  ±  10 

0.96  ±  0.03 

2 

286  ±  15 

(>340) 

0.91  ±0.03 

ESR-255 

245  ±  15 

(>419)d 

0.72  ±  0.03 

a.  After  1  hr  at  150  °C  followed  by  24  hrs  at  210  °C,  under  nitrogen.  Tg  indicated  by  peak  in 
loss  component  of  stiffness,  corrected  for  thermal  lag  and  heated  at  50  °C  /  min.  to  avoid  in-situ 
cure. 

b.  Determined  on  heating  at  10  °C/min.,  subsequent  to  initial  temperature  ramp  to  350  °C 

c.  Based  on  diBenedetto  equation  (see  Supporting  Information  Section  S4) 

d.  Heated  to  decomposition,  which  takes  place  either  before  or  simultaneous  with  mechanical 
softening 


In  order  to  better  understand  the  reasons  for  the  lack  of  complete  cure  in  ESR-255,  a 
more  detailed  analysis  of  the  cure  kinetics  based  on  the  Kamal  model47  analyzed  via  a  modified 
form  of  Kenny’s  graphical  method48  was  undertaken  using  isothermal  DSC  data  for  monomers  1, 
2,  and  ESR-255,  with  baselines  determined  in  a  manner  similar  to  that  reported  by  Sheng  et  al.44 
Full  details  and  complete  results  are  provided  in  Supporting  Information  Section  S6.  The 
activation  energy,  the  graphical  determination  of  which  is  depicted  in  Figure  6,  was  only  slightly 
lower  for  monomer  1  than  for  monomer  2  or  ESR-255  for  the  auto-catalytic  portion  of  the  cure, 
suggesting  that  the  rigidity  of  the  ESR-255  monomer  was  not  a  major  impediment  to  cure  prior 
to  gelation.  Interestingly,  the  initial  rates  of  conversion,  as  seen  in  Figure  7  at  290  °C,  for 
example,  were  always  significantly  higher  in  monomers  1  and  2  than  for  ESR-255.  The  final 
conversions  achieved  were  also  always  much  lower  in  ESR-255,  despite  ample  time  to  reach 
complete  conversion  based  on  the  observed  rate  of  cure  at  the  higher  temperatures. 
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Figure  6.  Arrhenius  plot  of  auto-catalytic  rate  constant  k2  (prior  to  gelation)  for  the  cyanate 
ester  monomers  studied.  Note  that,  because  the  exponents  for  each  monomer  are  different,  a 
higher  k2  value  does  not  necessarily  translate  into  a  faster  reaction  rate. 


Time(s) 


Figure  7.  Conversion  as  a  function  of  time  for  cyanate  esters  cured  at  290  °C  under  nitrogen. 


To  help  explain  the  observed  conversion  behavior,  an  estimate  of  the  Tg  of  the  cyanurate 
networks  during  the  isothermal  DSC  cure  was  created  (with  the  same  procedure  based  on  the 
diBenedetto  equation  used  for  Figure  5),  and  is  shown  in  Figure  8.  Although  the  slower  cure 
kinetics  of  ESR-255  lead  to  an  initially  slower  increase  in  Tg,  in  all  cases,  the  rise  in  Tg  slows 
dramatically  as  the  Tg  reaches  the  cure  temperature  (290  °C),  causing  the  sample  to  vitrify. 
Since  cyanate  esters  possess  the  capability  to  cure  slowly  in  the  glassy  state,  the  Tg  continues  to 
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rise  slowly,  even  after  vitrification.  In  monomer  1,  vitrification  does  not  take  place  until 
conversion  is  nearly  complete,  thus  further  conversion  is  limited  by  the  low  concentration  of 
remaining  unreacted  groups  and  further  increases  in  Tg  are  small.  By  contrast,  in  ESR-255, 
vitrification  takes  place  well  short  of  full  conversion,  and  thus  the  rate  of  conversion  in  the 
glassy  state  appears  to  be  somewhat  higher.  This  result,  when  combined  with  a  higher  sensitivity 
of  the  Tg  to  small  changes  in  conversion  (as  seen  in  the  diBenedetto  parameters  for  ESR-255), 
leads  to  a  comparatively  rapid  rise  in  Ts,  at  least  initially,  after  vitrification  in  ESR-255. 


Time  (s) 


Figure  8.  Estimated  Tg  as  a  function  of  time  for  cyanate  esters  cured  at  290  °C  under  nitrogen. 

Although  the  data  shown  in  Figure  8  are  for  relatively  short  cure  times  and  involve  Tg 
values  only  modestly  higher  than  the  cure  temperature,  the  data  in  Table  2  confirm  that,  if  given 
sufficient  time,  Tg  values  well  in  excess  of  the  cure  temperature  are  possible.  As  with  the  short¬ 
term  cure  at  higher  temperatures,  however,  it  would  appear  that  the  Tg  values  tend  to  converge 
over  time.  The  faster  cure  and  somewhat  higher  “as  cured”  (for  210  °C  for  24  hours  under 
nitrogen)  Tg  values  for  monomers  1  and  2  than  for  ESR-255  could  be  due  to  differing  levels  of 
impurities  in  the  samples,  van’t  Hoff  purity  determinations  (see  Supporting  Information,  Section 
S7)  showed  a  substantially  higher  apparent  purity  for  ESR-255.  However,  a  small  sample  of 
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monomer  1,  after  purification  by  flash  chromatography  through  silica  gel  in  an  ethyl 
acetate/hexane  mixture,  showed  no  differences  in  either  cure  kinetics  or  apparent  “as  cured”  Tg 
despite  a  significant  increase  in  the  van’t  Hoff  purity  estimate  (details  provided  in  Supporting 
Information  Section  S8). 

A  more  interesting  explanation  for  the  data  in  Table  2  is  that,  paradoxically,  more  flexible 
network  segments  in  cyanate  esters  allow  for  higher  “as  cured”  Tg  due  to  their  ability  to  facilitate 
continued  cure  after  vitrification,  as  described  in  the  Introduction.  In  fact,  somewhat  higher  “as 

1  O 

cured”  Tg  values  for  more  flexible  networks  have  been  reported  in  both  dicyanates  and  other 
tricyanates15  under  similar  cure  conditions.  Because  the  differences  involved  are  fairly  small, 
more  intensive  measurements,  perhaps  in  combination  with  recent  developments  in  the 
molecular  modeling  of  cyanurate  networks,49’50  will  be  needed  to  determine  with  certainty 
whether,  and  under  what  conditions,  a  more  flexible  network  can  routinely  provide  a  higher  Tg. 
Despite  being  much  higher  than  the  cure  temperature,  the  “as  cured”  Tg  values  for  these  networks 
are  much  more  similar  than  the  corresponding  set  of  fully-cured  Tg  values.  These  results  imply 
that  the  parameters  of  cure  process,  rather  than  the  rigidity  of  network  segments,  are  the  primary 
determinants  of  “as  cured”  Tg  in  these  networks. 

Other  Physical  Properties.  The  lower  cyanurate  ring  density  of  monomer  1  would  be 
expected  to  lead  to  a  lower  moisture  uptake  compared  to  monomer  2  and  ESR-255,  based  on 
previous  correlations  developed  for  cyanurate  networks  at  high  conversion.51  Indeed,  as  Table  3 
shows,  such  was  the  case.  In  fact,  the  moisture  uptake  of  cured  1,  at  just  2.2%  by  weight,  is 
among  the  lowest  observed  for  a  cyanurate  network  with  a  correspondingly  high  fully-cured  Tg, 
though  similar  tricyanates  have  come  fairly  close.  Although  the  lower  moisture  uptake  in 
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combination  with  the  similar  dry  Tg  of  cured  1  compared  to  cured  2  and  ESR-255  might  be 
expected  to  lead  to  a  significantly  higher  “wet”  Tg  for  cured  1  compared  to  both  cured  2  and 
ESR-255,  previous  work  on  cyanate  esters  has  shown  that  the  drop  in  Tg  on  exposure  to  hot 
water  becomes  significantly  less  pronounced  as  conversion  decreases.  Therefore,  the  lower 
conversion  of  cured  ESR-255  is  expected  to  compensate  somewhat  for  the  higher  moisture 
uptake,  resulting  in  more  modest  differences  in  “wet”  Ts,  as  observed  (see  Table  3). 


Table  3 

Other  Key  Physical  Properties  of  Cured  Cyanate  Esters 


Monomer 

“As  Cured” 
Density  (g/cc) 

Theoretical 
Cyanurate  Ring 
Density  at  Full 
Cure  (mmol  /  cc) 

Moisture  Uptake 
(wt%) 

“Wet”  Tg  (°C) 

1 

1.171  ±0.003 

2.48  ±  0.01 

2.2  ±  0.2 

238  ±  12 

2 

1.297 

3.40 

3.3 

232 

ESR-255 

1.270 

3.33 

3.5 

224 

Taken  together,  these  results  have  important  implications  for  the  manner  in  which 
thermosetting  polymers  are  selected  for  high-temperature  applications.  Resins  are  typically 
evaluated  solely  on  the  basis  of  a  “cured”  Tg,  with  catalyst  packages  then  selected  in  order  to 
attain  sufficient  cure  for  a  given  set  of  processing  conditions.  In  such  a  scenario,  a  higher  fully- 
cured  Tg  is  most  often  deemed  a  superior  choice.  This  prevalent  heuristic,  however,  ignores  both 
the  steep  dependence  of  the  Tg  on  the  extent  of  cure,  and  the  effect  of  different  catalyst  packages 
on  long-term  performance  degradation.  For  cyanurate  networks  in  particular,  these  two  effects 
are  especially  pronounced  in  comparison  to  other  types  of  thermosetting  polymer  networks,  and 
thus  the  prevalent  heuristic  is  less  appropriate. 
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When  there  is  a  steep  dependence  of  Tg  on  the  extent  of  cure,  it  is  most  likely,  as  seen  in 
the  data  presented  herein,  that  the  Tg  will  reach  and  eventually  exceed  the  cure  temperature, 
causing  the  cure  process  to  slow  dramatically.  As  a  result,  the  “as  cured”  Tg  of  the  network  will 
be  limited  largely  by  the  cure  conditions  themselves  rather  than  the  fully  cured  Tg,  provided  that 
the  fully  cured  Tg  is  high  enough.  As  a  result,  systems  with  a  higher  fully  cured  Tg  will  exhibit  a 
lower  extent  of  cure,  which  in  many  cases  is  detrimental  to  long-term  performance.  The  extent 
of  cure  (and  hence,  the  attainable  Tg )  may  be  increased  by  substitution  of  a  more  active  catalyst 
package,  with  the  trade-off  that  more  active  catalyst  packages  may  also  accelerate  degradation. 
Under  such  circumstances,  the  best  performing  network  may  not  be  the  one  with  the  highest  fully 
cured  Tg.  The  optimal  network  may  well  be  one  in  which  the  fully  cured  Tg  equals  the  maximum 
attainable  Tg  for  a  given  set  of  cure  conditions  and  a  catalyst  package  with  acceptable 
degradation  penalties. 

As  an  example,  consider  the  networks  formed  from  the  three  monomers  synthesized  in 
this  study.  If  the  cure  process  were  constrained  to  a  temperature  of  210  °C  and  an  available  time 
of  24  hours  (such  constraints  may  be  imposed  by,  for  instance,  by  the  thermal  stability  of  tooling 
or  resin  transfer  equipment),  then,  for  systems  with  no  added  catalyst,  the  maximum  Tg  attainable 
appears  to  be,  on  average,  roughly  250  -  280  °C.  In  such  cases,  a  network  with  a  fully  cured  Tg 
of  just  over  280  °C  would  be  ideal,  because  such  a  system  would  exhibit  both  the  highest 
attainable  Tg  and  full  cure.  Networks  with  a  higher  fully  cured  Tg,  such  as  ESR-255,  will,  after 
cure,  attain  the  same  “as  cured”  Tg  (250  -  280  °C),  but  will,  as  a  result,  exhibit  significantly  less 
than  full  conversion,  leading  to  potential  long-term  degradation.  Therefore,  for  many  practical 
situations,  particularly  with  cyanurate  networks,  there  will  exist  monomers  that  form  networks 
that  are  “too  rigid”,  that  is,  Tg  of  the  fully-cured  network  is  prohibitively  high.  The  discovery  of 
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monomer  1,  which  exhibits  a  lower  Tg  in  the  fully  cured  state  with  no  loss  in  thermochemical 
stability,  and  with  a  low  water  uptake,  thus  represents  a  significant  step  forward,  because  it  will 
allow  for  the  formation  of  optimal  networks  under  a  wider  variety  of  processing  conditions 
(particularly  those  involving  lower  temperatures  or  milder  catalysts)  than  has  previously  been 
possible. 

Conclusions 

The  selective  use  of  methylene  spacers  in  the  network  segments  of  highly  aromatic 
thermosetting  cyanate  ester  resins  has  proven  to  be  an  especially  valuable  technique  for 
providing  controlled  segmental  flexibility  while  maintaining  thermochemical  stability  and  low 
moisture  uptake.  In  particular,  a  newly  synthesized  monomer,  l,3,5-tris[(4- 
cyanatophenylmethyl)]benzene,  exhibits  char  yields  near  75%  at  600  °C  in  air,  equivalent  to  the 
best  known  common  cyanate  ester  resins,  with  a  moisture  uptake  of  only  2.2%  after  96  hours 
immersed  in  85  °C  water.  The  dry  Tg  of  the  fully-cured  network  produced  from  this  monomer 
was  around  320  °C.  Because  fully-cured  networks  prepared  from  similar  monomers  with  fewer 
methylene  spacers  exhibited  significantly  higher  Tg  s,  the  networks  did  show  the  expected  effects 
of  higher  segmental  flexibility.  However,  the  Tg  of  all  networks  studied  was  similar  after  a  24- 
hour  cure  at  210  °C,  as  was  the  “wet”  Tg.  Therefore,  in  practical  curing  and  use  situations,  the 
Tg  s  of  these  networks  are  likely  to  be  primarily  limited  by  parameters  of  the  cure  process  itself, 
specifically,  the  dramatic  slowing  of  the  cure  rate  that  occurs  upon  vitrification  and  the  degree  of 
catalysis,  regardless  of  the  degree  of  flexibility.  In  addition,  when  Tg  s  of  the  fully  cured 
networks  are  sufficiently  high,  the  more  flexible  networks  will  exhibit  the  highest  conversions 
for  a  given  cure  process.  Because  near  complete  conversion  of  cyanate  esters  typically  provides 
networks  with  superior  physical  properties,  in  this  case  the  most  flexible  network,  with  its  high 
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thermochemical  stability  and  low  water  uptake,  may  to  be  the  most  desirable.  These  results 
illustrate  that,  in  contrast  to  the  prevailing  heuristics  for  improving  the  performance  of  high- 
temperature  thermosetting  polymer  networks,  a  more  flexible  network  with  a  lower  fully-cured 
Tg  can  offer  an  optimal  combination  of  thermomechanical  and  thermochemical  performance. 
Both  the  newly  synthesized  materials,  as  well  as  the  new  insights  into  structure-process-property 
relationships  that  they  illustrate,  are  expected  to  help  advance  the  state  of  the  art  in  high- 
performance  polymeric  materials  for  the  electronics,  aerospace,  and  energy  storage  industries. 
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Experimental  Section: 

Synthesis  of  the  monomers  1,2  and  ESR-255 

General  consideration  and  instrumentation.  The  melting  points  were  collected  on  a  Mel- 
Temp  II  from  Laboratory  Devices  (Holliston,  MA)  and  are  not  corrected.  All  NMR  data  were 
collected  on  a  Bruker  Avance  II  300  MHz  spectrometer  (' H  at  300  MHz,  13C  at  75  MHz). 
Nuclear  magnetic  resonance  data  (free-induction  decays)  were  processed  using  NUTS  software 
from  Acorn  NMR  (Livermore,  CA).  All  spectra  are  referenced  to  solvent  or  tetramethylsilane. 
The  l,l,l-tris(4-hydroxphenyl)ethane  was  purchased  from  TCI  America  (Portland)  and  all  other 
reagents  were  purchased  from  Sigma-Aldrich  (Milwaukee).  Elemental  analyses  were  performed 
by  Atlantic  Microlab,  Inc.  (Norcross,  GA). 

Synthesis  of  monomer  1 

4-Methoxybenzoylacetaldehyde  potassium  salt  hydrate  (4) 

O  OK 


A  round-bottomed  flask  (500  mL)  equipped  with  magnetic  stirring  bar,  addition  funnel  and  N2 
bubbler  was  charged  with  KOtBu  (9.3  g,  83  mmol,  1  equiv)  and  anhydrous  THF  (100  mL). 
After  complete  dissolution,  4-methoxyacetophenone  (12.45  g,  83  mmol)  was  added  and  after  20 
min,  methyl  formate  (4.98  g,  5.1  mL,  83  mmol,  1  equiv)  was  added  dropwise.  The  mixture  was 
stirred  at  rt  for  2  h  and  then  diluted  with  hexanes  (200  mL).  The  mixture  was  filtered  through  a 
medium  porosity  glass  frit.  The  product  was  dried  to  a  constant  weight  in  a  vacuum  oven  (10 
torr,  50  °C)  overnight.  The  yield  was  15.3  g  (85  %).  Mp  >  240  °C.  *H  NMR  (300  MHz, 
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DMSO;  3,  ppm):  9.34  (d,  J=  9.6  Hz,  1H),  7.68  (d,  J=  8.5  Hz,  2H),  6.85  (d,  J=  8.7  Hz,  2H),  5.37 
(d ,J=  9.5  Hz,  1H),  3.75  (s,  3H).  Elemental  analysis  calcd  for  Ci0H9KO3  •  0.4  H20:  C,  53.74;  H, 
4.42;  found:  C,  53.36;  H,  4.07%. 


1 ,3,5-Tri(4-methoxybenzoyl)benzene  (5) 

O  O 


A  round-bottomed  flask  (250  mL)  equipped  with  magnetic  stirring  bar  and  reflux  condenser  was 
charged  with  4  (15.1  g,  68  mmol)  and  glacial  HO  Ac  (50  mL).  The  mixture  was  refluxed  for  30 
min.  After  cooling  to  rt,  the  precipitate  was  filtered  on  a  coarse  porosity  glass  frit. 
Recrystallization  from  MeCN  gave  the  title  compound  as  clear,  pale  yellow  rhombs.  The  yield 
was  9.2  g  (82  %).  Mp  179-181  °C  [lit.23  177  °C],  NMR  (300  MHz,  CDC13;  S,  ppm):  8.29  (s, 
3H),  7.86  (d,  J=  8.9  Hz,  6H),  6.98  (d,  J=  8.9  Hz,  6H),  3.89  (s,  9H).  13C  NMR  (75  MHz,  CDC13; 
3,  ppm):  193.98,  163.93,  138.85,  133.34,  132.79,  129.33,  114.08,  55.73.  Elemental  analysis 
calcd  for  C3oH2406:  C,  74.99;  H,  5.03;  found:  C,  74.87;  H,  5.01%. 

Racemic  l,3,5-Tris[(4-methoxyphenyl)hydroxymethyl]benzene  (6) 

OH  OH 


A  round-bottomed  flask  (2  L)  equipped  with  magnetic  stirring  bar,  reflux  condenser  and  N2 
bubbler  was  charged  with  5  (20  g,  41  mmol)  and  EtOH  (1  L).  The  mixture  was  refluxed  and 
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pulverized  NaBH4  (4.7  g,  125  mmol,  3  equiv)  was  added  in  portions  over  2  h.  Reflux  was 
continued  8  h  resulting  in  a  deep  red  solution.  The  mixture  was  cooled  to  rt  and  quenched  with 
HOAc.  The  mixture  was  poured  into  H20  (2  L)  and  the  white  precipitate  was  collected  on  a 
medium  porosity  glass  frit.  Recrystallization  from  MeCN  gave  the  title  compound  as  a  white, 
microcrystalline  powder.  The  yield  was  18  g  (90  %).  Mp  181-183  °C.  1 H  NMR  (300  MHz, 
DMSO;  S,  ppm):  7.27-7.13  (m,  9H),  6.83  (d,  J=  9.4  Hz,  6H),  5.67  (d,  J=  3.7  Hz,  3H),  5.55  (d,  J 
=  3.7  Hz,  3H),  3.70  (s,  9H).  13C  NMR  (75  MHz,  DMSO;  S,  ppm):  158.00,  145.32,  137.82, 
122.39,  113.32,  74.01,  54.98.  Elemental  analysis  calcd  for  C30H30O6:  C,  74.06;  H,  6.21;  found: 
C,  74.04;  H,  6.26%. 


1 ,3,5-Tris[(4-methoxyphenyl)methyl]benzene  (7) 


A  round-bottomed  flask  (500  mL)  equipped  with  magnetic  stirring  bar,  reflux  condenser  and  N2 
bubbler  was  charged  with  6  (12  g,  25  mmol),  5  %  Pd/C  catalyst  (2  g),  ammonium  formate  (8.8  g, 
140  mmol,  6  equiv)  and  glacial  HOAc  (250  mL).  The  mixture  was  refluxed  for  5  h.  After 
cooling  to  rt,  the  mixture  was  filtered  through  diatomaceous  earth  to  remove  the  catalyst.  The 
filtrate  was  partitioned  between  H20  (500  mL)  and  Et20  (500  mL).  The  organic  layer  was 
separated  and  washed  with  H20  (500  mL)  and  then  brine  (500  mL).  Rotary  evaporation  of  the 
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volatiles  gave  a  colorless  oil  that  gradually  solidified  into  a  white,  waxy  solid.  The  yield  was 
9.87  g  (91  %).  Mp  55-57  °C.  lU  NMR  (300  MHz,  DMSO;  S,  ppm):  7.07  (d,  J=  8.6  Hz,  6H), 
6.85  (s,  3H),  6.81  (d ,J=  8.6  Hz,  6H),  3.76  (s,  6H),  3.69  (s,  9H).  13C  NMR  (75  MHz,  DMSO;  S, 
ppm):  157.48,  141.71,  133.09,  129.49,  126.58,  113.72,  54.90,  54.81,  40.16.  Elemental  analysis 
calcd  for  C30H30O3:  C,  82.16;  H,  6.89;  found:  C,  82.08;  H,  6.88%. 

1 ,3,5-Tris[(4-hydroxyphenyl)methyl]benzene  (8) 


A  round-bottomed  flask  (500  mL)  equipped  with  magnetic  stirring  bar  and  reflux  condenser  was 
charged  with  7  (9.87  g,  22  mmol)  and  pyridine  hydrochloride  (50  g,  434  mmol,  19  equiv).  The 
mixture  was  refluxed  for  4  h.  The  mixture  was  carefully  diluted  with  H2O  (300  mL)  whereby  a 
white  precipitate  formed.  The  precipitate  was  collected  on  a  medium  porosity  glass  frit.  The 
precipitate  was  dissolved  in  EtOAc  (100  mL)  and  washed  with  H2O  (100  mL)  followed  by  brine 
(100  mL).  After  drying  over  anhydrous  MgSCL,  the  solvent  was  rotary  evaporated  leaving  a 
white  solid.  Recrystallization  from  MeOHdrLO  gave  the  title  compound  colorless  plates.  The 
yield  was  8  g  (90  %).  Mp  205-207  °C.  NMR  (300  MHz,  DSMO;  S,  ppm):  9.32  (s,  30H), 
6.93  (d,  J=  8.6  Hz,  6H),  6.81  (s,  3H),  6.62  (d,  J=  8.4  Hz,  6H),  3.69  (s,  6H).  13C  NMR  (75  MHz, 
DMSO;  S,  ppm):  155.61,  142.24,  131.86,  129.82,  126.78,  115.46,  40.58.  Elemental  analysis 
calcd  for  C27H24O3:  C,  81.79;  H,  6.10;  found:  C,  81.48;  H,  6.16%. 

1 ,3,5-Tris[(4-cyanatophenyl)methyl]benzene  (1) 
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A  round-bottomed  flask  (250  mL)  equipped  with  magnetic  stirring  bar,  N2  bubbler  and  addition 
funnel  was  charged  with  8  (3  g,  7  mmol),  BrCN  (2.81  g,  26  mmol,  3.5  equiv)  and  acetone  (100 
mL).  The  mixture  was  cooled  in  a  -20  °C  bath  before  TEA  (2.1  g,  21  mmol,  3  equiv)  was  added 
dropwise  in  20  min.  Copious  solids  (TEA-HBr)  precipitated  during  the  addition.  After  2h,  the 
mixture  was  poured  into  cold  H2O  (250  mL)  and  extracted  with  Et20  (2X  100  mL).  The  extracts 
were  collected  and  washed  with  H2O  (200  mL)  followed  by  saturated  aqueous  NaHCCL  (100 
mL)  and  finally  brine  (100  mL).  After  drying  over  anhydrous  MgSCL,  rotary  evaporation  left  an 
off-white  solid.  Recrystallization  from  EtOH  gave  the  title  compound  as  white  needles.  The 
yield  was  1.97  g  (60  %)  Mp  97-99  °C.  NMR  (300  MHz,  CDC13;  S,  ppm):  7.21  (s,  12H),  6.81 
(s,  3H),  3.91  (s,  6H).  13C  NMR  (75  MHz,  CDCI3;  S,  ppm):  151.53,  141.19,  139.99,  130.85, 
127.75,  115.54,  109.05,  40.99.  Elemental  analysis  calcd  for  C30H21N3O3:  C,  76.42;  H,  4.49;  N, 
8.91;  found:  C,  76.19;  H,  4.37;  N,  8.89%. 

Synthesis  of  Monomer  2 
Dimethyl  5-methoxyisophthalate  (9) 


A  round-bottomed  flask  (1  L)  equipped  with  magnetic  stirring  bar  and  reflux  condenser  was 
charged  with  5-hydroxyisophthalic  acid  (18.2  g,  100  mmol),  acetone  (300  mL),  freshly 
pulverized  anhydrous  K2CO3  (40  g,  289  mmol,  2.8  equiv)  and  dimethyl  sulfate  (50  mL,  66.5  g, 
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520  mmol,  5.2  equiv).  After  refluxing  overnight,  the  mixture  was  cooled  and  poured  into  H20  (1 
L).  The  resulting  precipitate  was  collected  on  a  medium  porosity  glass  frit  and  air-dried  several 
hours.  Recrystallization  from  cyclohexane  gave  the  title  compound  as  a  white  microcrystalline 
powder.  The  yield  was  21.1  g  (94  %).  Mp  110-112  °C  [lit.53  110-111  °C].  Elemental  analysis 
calcd  for  CnHi205:  C,  58.93;  H,  5.39;  found:  C,  58.66;  H,  5.31%. 

5-Methoxyisophthalic  acid  monohydrate  (10) 


A  round-bottomed  flask  (500  mL)  equipped  with  magnetic  stirring  bar  was  charged  with  9  (21.1 
g,  94  mmol),  MeOH  (200  mL)  and  a  solution  of  KOH  (16  g,  285  mmol,  3  equiv)  in  H20  (50 
mL).  The  mixture  was  refluxed  for  6  h.  After  cooling  to  rt,  the  mixture  was  poured  into  H20 
(500  mL)  and  acidified  to  pH  3  by  dropwise  addition  of  cone.  HC1.  The  precipitate  was 
collected  on  a  medium  porosity  glass  frit.  Recrystallization  from  glacial  HOAc  gave  the  title 
compound  as  a  colorless  needles.  The  yield  was  19.1  g  (95  %).  Mp  270-272  °C  [lit.54  267-268 
°C],  lH  NMR  (300  MHz,  DMSO;  S,  ppm):  13.20  (bs,  2  C02H),  8.09  (s,  1H),  7.61  (s,  2H),  3.83 
(s,  3H).  13C  NMR  (75  MHz,  DMSO;  S,  ppm):  167.11,  159.37,  133.74,  122.56,  118.28,  55.65. 
Elemental  analysis  calcd  for  C9Hg05-H20:  C,  50.47;  H,  4.71;  found:  C,  50.92;  H,  4.63%. 

3 , 5  -Bis(4-methoxybenzoyl)anisole  (11) 

O  O 


11 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


32 


A  round-bottomed  flask  (200  mL)  equipped  with  magnetic  stirring  bar  and  reflux  condenser  was 
charged  with  10  (18.2  g,  85  mmol)  and  SOCI2  (25  mL).  An  N2  bubbler  piped  to  a  caustic  wash 
bottle  was  attached  and  the  mixture  was  gently  refluxed  for  5  h.  The  excess  SOCl2  was  distilled 
under  reduced  pressure  (10  torr)  leaving  an  off-white  solid  of  5-methoxyisophthaloyl  dichloride. 
The  latter  was  transferred  to  a  500  mL  round-bottomed  flask  equipped  with  magnetic  stirring  bar 
and  addition  funnel  with  carbon  disulfide  (50  mL).  The  mixture  was  cooled  in  an  ice  bath  before 
AICI3  (17.69  g,  134  mmol,  2  equiv)  was  added  in  one  portion.  The  addition  funnel  was  charged 
with  anisole  (29  g,  270  mmol,  4  equiv)  which  was  added  over  dropwise  over  30  min.  The 
cooling  bath  was  removed  and  the  mixture  was  stirred  at  rt  for  3  h.  The  mixture  was  carefully 
poured  onto  cracked  ice  (500  g).  The  mixture  was  extracted  with  EtOAc  (3X  100  mL).  The 
extracts  were  collected  and  washed  with  H20  (500  mL)  followed  by  brine  (500  mL).  After 
drying  over  anhydrous  MgS04,  the  solvent  was  rotary  evaporated  leaving  a  pale  yellow  oil  that 
eventually  solidified  under  vacuum.  Recrystallization  of  the  crude  product  from  EtOH  gave  the 
title  compound  as  small  colorless  needles.  The  yield  was  16.61  g  (52  %).  Mp  115-117  °C.  1 H 
NMR  (300  MHz,  CDC13;  S,  ppm):  7.84  (d,  J=  8.6  Hz,  4H),  7.62  (t,  J=  1.4  Hz,  1H),  7.49  (d,  J  = 
1.2  Hz,  2H),  6.96  (d,  J=  9.0  Hz,  4H),  3.91  (s,  3H),  3.88  (s,  6H).  13C  NMR  (75  MHz,  CDC13;  S, 
ppm):  194.57,  163.69,  159.71,  139.76,  132.72,  129.78,  123.36,  118.27,  113.88,  55.94,  55.66. 
Elemental  analysis:  calcd  for  C23H2oC>5:  C,  73.39;  H,  5.36;  found:  C,  73.65;  H,  5.23%. 

Racemic  3,5-Bis[(4-methoxyphenyl)hydroxymethyl)]anisole  (12) 

OH  OH 


12 
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A  round-bottomed  flask  (1  L)  equipped  with  magnetic  stirring  bar  was  charged  with  LAH  (3.3  g, 
87  mmol,  1.9  equiv)  and  anhydrous  THF  (500  mL).  An  N2  bubbler,  addition  funnel  and  reflux 
condenser  were  equipped.  The  mixture  was  stirred  in  an  ice  bath  and  a  solution  of  11  (16.61  g, 
44  mmol)  in  THF  (50  mL)  was  charged  to  the  funnel.  The  addition  was  made  as  a  gentle  stream 
over  1  h.  After  stirring  at  rt  for  3  h  the  reduction  was  complete.  The  mixture  was  carefully 
quenched  by  addition  of  H20  (3.3  mL)  followed  by  15  %  NaOH  (3.3  mL)  and  finally  H20  (9.9 
mL).  After  stirring  overnight,  the  mixture  was  filtered  through  a  medium  porosity  glass  frit.  The 
filtrate  was  rotary  evaporated  leaving  a  crude  solid.  Recrystallization  from  toluene  gave  the  title 
compound  as  a  white  powder.  The  yield  was  10.1  g  (60  %).  1 H  NMR  (300  MHz, 

CDCf/DMSO;  d,  ppm):  7.27-7.20  (m,  4H),  7.02-6.93  (m,  1H),  6.82-6.74  (m,  6H),  5.63  (s,  2H), 
4.22  (bs,  2  OH),  3.73  (s,  6H),  3.69  (s,  3H).  13C  NMR  (75  MHz,  CDCI3/DMSO;  S,  ppm):  159.43, 
159.35,  158.41,  158.40,  146.43,  146.37,  136.87,  136.85,  127.63,  127.59,  117.09,  116.99,  113.31, 
110.43,  110.31,  74.83,  74.81,  54.97,  54.93.  Elemental  analysis  calcd  for  C23H2405:  C,  72.61;  H, 
6.36;  found:  C,  72.88;  H,  6.42%. 

3,5-Bis[(4-methoxyphenyl)methyl)]anisole  (13) 

OMe 

13 

A  round-bottomed  flask  (500  mL)  equipped  with  magnetic  stirring  bar  and  reflux  condenser  was 
charged  with  a  mixture  of  12  (10.1  g,  26  mmol),  glacial  HOAc  (300  mL)  and  5  %  Pd/C  (2  g). 
The  mixture  was  vigorously  stirred  and  heated  near  boiling  under  an  H2  atmosphere  (1  torr). 
After  1 8  h,  the  mixture  was  cooled  to  rt  and  the  catalyst  was  filtered.  Rotary  evaporation  left  an 
oil  which  was  chromatographed  on  silica  gel  eluting  with  20  %  EtOAc  in  hexanes.  Rotary 
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evaporation  of  the  fractions  gave  a  colorless  oil.  Reduced  pressure  (0. 1  torr)  distillation  gave  the 
title  compound  in  analytically  pure  form.  The  yield  was  3.87  g  (41  %).  1 H  NMR  (300  MHz, 
CDC13;  3,  ppm):  7.08  (d,  J=  8.6  Hz,  4H),  6.81  (d,  J=  8.8  Hz,  4H),  6.62  (t,  J=  1.4  Hz,  1H),  6.53 
(d ,J=  1.4  Hz,  2H),  3.84  (s,  4H),  3.77  (s,  6H),  3.69  (s,  3H).  13C  NMR  (75  MHz,  CDC13;  3,  ppm): 
160.10,  158.21,  143.28,  133.33,  130.04,  122.21,  114.09,  112.35,  55.46,  55.31,  41.26.  Elemental 
analysis  calcd  for  C23H2403:  C,  79.28;  H,  6.94;  found:  C,  79.31;  H,  6.94%. 
3,5-Bis[(4-hydroxyphenyl)methyl)]phenol  (14) 


OH 

14 


A  round-bottomed  flask  (250  mL)  equipped  with  magnetic  stirring  bar  and  reflux  condenser  was 
charged  with  13  (3.87  g,  1 1  mmol)  and  pyridine  hydrochloride  (26  g,  222  mmol,  20  equiv).  The 
mixture  was  heated  by  heating  mantle  to  180  °C  for  5  h.  After  cooling  to  rt,  the  mixture  was 
diluted  with  H20  (250  mL)  whereupon  a  tan  precipitate  eventually  formed.  Filtration  followed 
by  recrystallization  of  the  solid  from  toluene/MeOH  gave  the  title  compound  as  a  white, 
microcrystalline  powder.  The  yield  was  2.5  g  (75  %).  Mp  183-186  °C.  *H  NMR  (300  MHz, 
CDC13;  3,  ppm):  9.17  (2,  20H),  9.1 1  (s,  10H),  6.97  (d,  J=  8.3  Hz,  4H),  6.67  (d ,J=  8.0  Hz,  4H), 
6.50  (s,  1H),  6.35  (s,  2H),  3.67  (s,  4H).  13C  NMR  (75  MHz,  CDC13;  3,  ppm):  157.33,  155.48, 
143.13,  131.35,  129.59,  119.72,  115.11,  113.06,  40.32.  Elemental  analysis  calcd  for  C20Hi8O3: 
C,  78.41;  H,  5.92;  found:  C,  78.31;  H,  5.78%. 

3,5-Bis[(4-cyanatophenyl)methyl)]phenylcyanate  (2) 
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A  round-bottomed  flask  (250  mL)  equipped  with  magnetic  stirring  bar  was  charged  with  14  (1  g, 
3.3  mmol),  BrCN  (1.38  g,  13  mmol,  3.9  equiv)  and  acetone  (50  mL).  The  mixture  was  stirred  in 
a  -20  °C  cooling  bath  while  TEA  (1.07  g,  10.6  mmol,  3.2  equiv)  was  added  dropwise.  Copious 
solids  (TEA-HBr)  precipitated  during  the  addition.  After  stirring  2  h,  the  mixture  was  poured 
into  H2O  (100  mL).  The  mixture  was  extracted  with  EtOAc  (2X  50mL).  The  extracts  were 
combined  and  washed  with  H2O  (100  mL)  followed  by  brine  (100  mL).  After  drying  over 
anhydrous  MgSCL,  the  solvent  was  rotary  evaporated  leaving  the  crude  product.  The  title 
compound  was  obtained  as  colorless  needles  by  recrystallization  from  z'-PrOH.  The  yield  was 
1.2  g  (95  %).  Mp  1 19-121  °C.  *H  NMR  (300  MHz,  CDCI3;  S,  ppm):  7.24  (m,  8H),  6.94  (s,  2H), 
6.92  (s,  1H),  3.99  (2,  4H).  13C  NMR  (75  MHz,  CDC13;  S,  ppm):  153.47,  151.86,  143.97,  138.55, 
130.97,  127.92,  115.89,  113.95,  108.89,  108.71,  40.82.  Elemental  analysis  calcd  for 
C23H15N3O3:  C,  72.43;  H,  3.96;  N,  11.02;  found:  C,  72.43;  H,  3.91;  N,  10.98%. 

Synthesis  ofESR-255 

1,1,1  -Tris(4-cyanatophenyl)ethane  (ESR-255) 


A  round-bottomed  flask  (250  mL)  equipped  with  magnetic  stirring  bar  and  addition  funnel  was 
charged  with  l,l,l-tris(4-hydroxyphenyl)ethane  (10.0  g,  32  mmol),  BrCN  (13.78  g,  130  mmol,  4 
equiv)  and  anhydrous  acetone  (100  mL).  The  mixture  was  cooled  in  a  -20  °C  bath  before  TEA 
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(11.55  g,  114  mmol,  3.5  equiv)  was  added  dropwise.  During  the  addition,  copious  solids 
(TEA-HBr)  precipitated.  Afterwards,  the  mixture  was  stirred  in  an  ice  bath  for  30  min  before 
pouring  into  crushed  ice  and  H20  (500  mL).  The  mixture  was  extracted  with  Et20  (3X  100  mL). 
The  extracts  were  collected  and  washed  with  H20  (100  mL)  followed  by  brine  (100  mL).  After 
drying  over  anhydrous  MgS04,  the  solvent  was  rotary  evaporated  leaving  an  off-white  solid 
(12.59  g).  The  crude  product  was  recrystallized  from  EtOH  to  give  the  title  compound  as  a  white 
powder.  The  yield  was  6.5  g  (52  %).  Mp  113-115  °C  [lit.55  104  °C],  ‘H  NMR  (300  MHz, 
CDC13;  8,  ppm):  7.28  (d,  J=  9.1  Hz,  6H),  7.17  (d,  J=  9.0  Hz,  6H),  2.22  (s,  3H).  13C  NMR  (75 
MHz,  CDC13;  8,  ppm):  151.48,  146.58,  130.57,  115.34,  108.65,  51.69,  30.92.  Elemental  analysis 
calcd  for  C23Hi5N303:  C,  72.43;  H,  3.96;  N,  11.02;  found:  C,  72.35;  H,  3.88;  N,  11.06%. 

Cured  network  preparation  and  characterization.  Cured  resin  samples  were  prepared 
by  melting  the  monomer  powder  at  115  -  130  °C  at  a  reduced  pressure  of  300  mm  Hg  for  30 
minutes,  then  pouring  into  either  a  silicone  casting  mold  (13  mm  diameter  x  3  mm  discs)  for 
multi-step  cure,  or  into  a  6  mm  diameter  x  1  mm  high  aluminum  pan  for  single-step  cure.  Multi- 
step  cure  involved  heating  in  an  oven  under  nitrogen  to  150  °C  for  1  hour,  followed  by  210  °C 
for  24  hours,  with  subsequent  cooling  to  below  150  °C  prior  to  de-molding.  Heating  ramps  for 
multi-step  cures  were  5  °C  /  min.  Additional  details  such  as  the  mold  fabrication  procedures 
have  been  published  elsewhere.56 

Thermomechanical  analysis  of  cured  samples  was  performed  using  a  TA  Instruments 
Q400  thermomechanical  analyzer  (TMA)  in  dynamic  (oscillatory  compression)  TMA  mode. 
The  oscillatory  compression  mode  provides  qualitative  information  on  the  storage  and  loss 
components  of  the  sample  stiffness,  in  addition  to  sample  displacement  as  a  function  of 
temperature.  Heating  and  cooling  rates  of  50  °C  /  min  were  used  for  “as  cured”  dry  samples, 
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with  10  °C  used  for  “fully  cured”  samples,  and  20  °C  /  min  used  for  “wet”  samples.  These  rates 
are  considered  optimal  for  each  given  situation,  and  represent  a  balance  between  minimizing  in- 
situ  cure,  thermal  lag,  and  (for  wet  samples)  in-situ  drying.  A  standard  thermal  cycling 
procedure,  using  limits  of  0  °C  and  200  °C,  was  used  to  determine  and  correct  for  the  thermal 
lag.  The  procedure  was  employed  prior  to  heating  to  350  °C  for  dry  samples,  and  after  heating  to 
350  °C  for  wet  samples.  The  mean  compressive  load  on  the  samples  was  0.1  N,  with  an 
oscillatory  force  applied  at  an  amplitude  of  0. 1  N  and  a  frequency  of  0.05  Hz.  Complete  details, 
examples,  and  the  rationale  for  the  cycling  procedures  have  been  explained  at  length  elsewhere.56 
Thermogravimetric  analysis  (TGA)  was  carried  out  using  a  TA  Instrument  Q5000  under  60  mL  / 
min.  of  either  nitrogen  or  air.  For  TGA  analysis,  ~  2  mg  of  uncured  powder  was  heated  at  10  °C 
/  min.  to  600  °C.  The  density  of  cured  samples  was  determined  via  neutral  buoyancy  of  13  mm  x 
3  mm  discs  in  CaCl2  /  deionized  water  mixtures,  as  described  elsewhere.57  Some  multi-step 
cured  and  post-cured  samples  were  placed  in  250  mL  of  deionized  water  at  85  °C  for  96  hours  as 
a  means  of  testing  the  effects  of  exposure  to  hot  water. 
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SI.  X-ray  Crystal  Structure  and  NMR  Data 

Figure  Sl-1  shows  the  X-ray  crystal  structure  of  l,3,5-tris(4-hydroxyphenylmethyl)benzene,  compound  8 
in  the  main  manuscript,  and  the  precursor  that  was  cyanated  to  produce  monomer  1  .  CCDC  821411 
contains  the  supplementary  crystallographic  data  for  compound  8.  This  data  can  be  obtained  free  of 
charge  from  The  Cambridge  Crystallographic  Data  Centre  via  www.ccdc.cam.ac.uk/data  request/cif,  by 
e-mailing  data  request(5)ccdc.cam.ac.uk.  or  by  contacting  CCDC  12  Union  Road,  Cambridge  CB2  1EZ,  UK; 
fax:  +44  1223  336033.  Figures  Sl-2  through  Sl-4  show  the  3H  (300  MHz,  CDCI3)  and  13C  (75  MHz,  CDCI3) 
NMR  spectra  of  monomers  1,  2,  and  ESR-255,  respectively. 
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Figure  Sl-1.  Thermal  ellipsoid  plot  of  l,3,5-tris(4-hydroxyphenylmethyl)benzene,  compound  8. 
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Figure  Sl-2.  1H  (upper)  and  13C  (lower)  NMR  spectra  of  l,3,5-tris(4-cyanatophenylmethyl)benzene, 
compound  1. 
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Figure  Sl-3.  1H  (upper)  and  13C  (lower)  NMR  spectra  of  3,5-Bis[(4- 

cyanatophenyl)methyl)]phenylcyanate,  compound  2. 
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Figure  Sl-4.  1H  (upper)  and  13C  (lower)  NMR  spectra  of  l,l,l-Tris(4-cyanatophenyl)ethane,  ESR-255. 
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S2.  Modulated  DSC  Data 


52.1  Introduction.  The  DSC  traces  for  ESR-255  (seen  in  Figure  5c  of  the  main  manuscript,  for  instance) 
show  an  apparently  "open-ended"  exotherm  associated  with  cure  of  the  cyanate  ester  groups.  Based 
only  on  conventional  DSC  scans,  however,  it  is  not  possible  to  differentiate  between  a  truly  open-ended 
exotherm  in  which  cure  continues  even  in  the  "flat"  portion  of  the  curve  and  a  simple  shift  of  the 
baseline.  Using  modulated  DSC,  however,  it  is  possible  to  distinguish  among  these  alternative 
possibilities.  Continuing  cure  of  the  cyanate  ester  groups  would  affect  only  the  non-reversing  portion  of 
the  modulated  DSC  curve,  while  a  baseline  shift  would  affect  the  reversing  portion  (and  possibly  the 
non-reversing  portion,  depending  on  the  cause). 

52.2  Experimental  Method.  Samples  of  ESR-255  (~3  mg)  were  heated  at  varying  rates  in  a  TA 
Instruments  Q200  Differential  Scanning  Calorimeter  in  modulated  DSC  mode  under  50  mL/min  flowing 
N2.  A  modulation  amplitude  of  1  °C  and  a  modulation  period  of  60  s  were  used  over  a  temperature 
range  of  50  °C  to  350  °C.  Since  quantitative  heat  flow  information  was  not  required,  the  data  shown 
have  not  been  corrected  using  a  re-scanned  baseline. 

52.3  Results  and  Discussion.  Figure  S2-1  shows  the  modulated  DSC  of  ESR-255  heated  at  a  mean  rate 
of  3  °C  /  min  (this  sample  was  melted  just  prior  to  starting  the  DSC  run).  Though  there  is  a  small  shift  in 
the  reversible  component  of  the  heat  flow  due  to  vitrification  (a  kind  of  "glass  transition  in  reverse"  that 
causes  a  step  decrease  in  heat  capacity),  the  "open-ended"  heat  flow  curve  is  clearly  due  to  a  continuing 
non-reversible  heat  flow.  In  other  words,  cure  of  the  cyanate  ester  definitely  continues  even  in  the 
"flat"  portion  of  the  curve.  Since  degradation  of  the  network  begins  to  have  a  significant  effect  on  DSC 
scans  above  350  °C  (an  upturn  in  the  non-reversible  heat  flow  due  to  this  effect  is  just  visible  at  the  end 
of  the  scan),  there  is  no  practical  way  to  measure  the  total  heat  of  cyclotrimerization  of  uncatalyzed 
ESR-255  (or  any  other  similarly  rigid  thermosetting  polymer  requiring  high  cure  temperatures)  using 
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traditional  DSC  techniques  alone.  Thus,  a  method  such  as  the  combined  DSC/IR  technique  presented  in 


Section  S3  is  required.  Since  the  heating  rate  used  in  Figure  5  in  the  main  manuscript  was  somewhat 
different  than  the  mean  heating  rate  used  for  modulated  DSC  (which  requires  much  slower  heating  in 
general),  the  effect  of  heating  rate  changes  was  checked.  Figure  S2-2  shows  a  modulated  DSC  scan  of 
ESR-255  that  is  identical  to  the  one  used  to  produce  Figure  S2-1,  except  that  the  mean  heating  rate  was 
increased  to  5  °C/min,  and  in  the  case  pre-melting  was  not  carried  out.  The  mean  heating  rate  used  is 
high  for  modulated  DSC,  and  hence  the  reversible  heating  curve  shows  some  noise  at  the  vitrification 
point,  however,  it  is  clear  that  the  "open-ended"  curve  is  still  due  to  continuing  non-reversible  heat 
flow.  Flence,  the  necessary  difference  in  heating  rates  did  not  significantly  affect  the  key  result,  which  is 
that  continuing  cure,  rather  than  a  baseline  shift,  is  responsible  for  the  "open-ended"  L-shaped 
exotherm  seen  in  ESR-255. 


Figure  S2-1.  Modulated  DSC  of  previously  melted  ESR-255  at  a  mean  heating  rate  of  3  °C  /  min.  Note 
that  the  shift  in  heat  flow  after  the  exothermic  peak  is  due  to  mainly  to  the  non-reversing  component. 
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Figure  S2-2.  Modulated  DSC  of  ESR-255  at  a  mean  heating  rate  of  5  °C  /  min. 
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S3.  Combined  FT-IR/DSC  Method  for  the  Quantitative  Determination  of  Absolute  Conversions  in 


Cyanate  Esters 

S3.1.  Introduction.  In  most  thermosetting  systems,  the  absolute  conversion  (a)  may  be  determined 
most  simply  through  the  use  of  DSC  by  making  use  of  the  formula 

a  =  (  AH0- AHr)  /  AH0  (S3-1) 

in  which  AH0  represents  the  heat  of  reaction  of  the  monomer,  determined  by  integrating  the  exothermic 
reaction  peak  seen  in  a  first-heating  DSC  scan,  and  AHr  is  the  residual  heat  of  reaction  of  a  partially  cured 
sample  determined  by  integrating  the  exothermic  reaction  peak  using  the  same  heating  program  and 
baseline  determination  technique  applied  in  the  calculation  in  AH0.  The  key  assumptions  underlying  the 
method  include  direct  proportionality  between  the  heat  released  and  the  extent  of  conversion  and  the 
attainment  of  complete  conversion  without  significant  side  reactions.  For  cyanate  esters,  the  former 
assumption  is  almost  always  valid,  since  cyanate  esters  tend  to  exhibit  a  characteristic  enthalpy  of  -110 
kJ/mol  eq.  on  conversion  of  a  cyanate  ester  group. S1  For  cyanate  esters  with  a  glass  transition 
temperature  significantly  in  excess  of  300  °C  or  with  significant  steric  hindrance  to  cyanurate  formation, 
however,  the  latter  assumption  may  be  problematic,  since  significant  thermal  degradation  of  cyanurate 
networks  is  often  detected  by  DSC  at  or  above  about  350  °C.  In  such  cases,  a  DSC  heating  scan  of  the 
monomer  terminated  at  350  °C  or  lower  will  not  lead  to  full  conversion  of  the  monomer,  while  one 
terminated  above  350  °C  will  detect  heat  from  both  conversion  and  thermal  degradation.  A  good 
indicator  that  the  assumption  of  full  cure  with  no  side  reactions  has  been  violated  is  the  presence  of  an 
"open  ended"  DSC  trace  in  which  the  heat  flow  signal  does  not  return  to  its  baseline  vvalue  prior  to 
termination  of  the  heating.  Note  that  many  DSC  traces  of  this  type  exhibit  an  "L-shaped"  heating  curve 
that  resembles  a  full  peak  with  a  shifted  baseline,  thus  care  must  be  taken  to  establish  the  actual 
baseline  (usually  by  re-heating  the  sample  at  the  same  rate  immediately  afterword  without  moving  or 
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handling  the  sample)  when  using  DSC  to  assess  the  heat  of  reaction  for  high-temperature  thermosetting 


monomers. 

In  addition  to  DSC,  infrared  spectroscopy  (FT-IR  or  ATR-FTIR)  has  often  been  utilized  to  assess 
the  conversion  of  cyanate  ester  to  cyanurate.  Infrared  spectroscopy  has  the  key  advantage  that  residual 
cyanate  ester  groups  are  always  visible,  due  to  the  presence  of  strong  characteristic  absorption  bands  in 
the  otherwise  sparsely  occupied  region  of  the  infrared  spectrum  between  2200  and  2300  cm"1. 
However,  infrared  spectroscopy  is  often  not  considered  a  quantitative  technique  unless  calibrated  by 
DSC  or  other  measurements,  because  the  absorption  intensity  for  cyanate  ester  peaks  does  not 
decrease  linearly  with  conversion.51  The  reason  for  the  nonlinearity,  however,  is  usually  a  change  in  the 
physical  state  of  the  sample  from  crystalline  solid  or  small  molecular  liquid  to  oligomeric  or  polymeric 
gel  to  rubbery  or  glassy  network  sol  and  finally  to  glassy  or  rubbery  network  with  residual  cyanate  esters 
as  dangling  chain  ends.  Thus,  over  limited  ranges  of  conversion  in  which  the  material  remains  in  the 
same  general  state,  such  as  near  the  end  of  conversion  in  which  the  system  consists  of  a  glassy  network 
with  cyanate  esters  as  unreacted  chain  ends,  there  does  exist  a  linear  relationship  between  the 
integrated  intensity  of  the  cyanate  ester  peaks  and  their  concentration.  For  most  thermosetting 
polymers  for  which  DSC  provides  an  absolute  measure  of  conversion,  DSC  data  points  are  simply  used  to 
calibrate  the  infrared  absorption  curve. 

For  high-temperature  cyanate  esters  that  cannot  be  fully  cured  without  thermal  degradation  in 
a  DSC  scan,  though,  it  is  still  possible  to  use  DSC  data  points  to  calibrate  an  infrared  absorption  curve,  by 
the  method  outlined  in  Section  S3. 2.  This  method  relies  on  the  existence  of  a  region  of  conversion  in 
which  the  infrared  absorption  is  a  linear  function  of  cyanate  ester  concentration  (with  the  linearity 
extending  to  100%  conversion)  and  that  can  be  obtained  by  heating  in  a  DSC  scan  under  conditions  that 
permit  quantitative  determination  of  the  heat  of  reaction.  The  heat  of  reaction  may  not  be  determined 
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with  quantitative  precision  if  either  the  heating  rate  of  the  DSC  scan  is  too  slow,  or  if  the  isothermal 


portion  of  a  DSC  scan  continues  for  so  long  that  the  conversion  rate  becomes  too  low  to  measure 
reliably.  The  essence  of  the  method  is  that  the  linearity  of  both  infrared  and  DSC  signals  allows  for  the 
determination  of  the  differential  of  the  heat  of  reaction  with  respect  to  infrared  absorption,  enabling 
the  conversion  of  the  residual  infrared  absorption  (which  is  known  absolutely)  to  an  absolute  residual 
heat  of  reaction  (AHr).  When  added  to  the  known  value  of  the  evolved  heat  of  reaction  for  a  given 
sample,  the  value  of  AH0  may  be  determined,  and,  in  conjunction  with  the  AHr  value  and  Eq.  S3-1,  the 
absolute  conversion  of  the  sample  may  be  computed.  The  procedure  is  outlined  in  Section  S3-2,  while 
validation  examples  (using  cyanates  that  can  be  cured  completely  in  the  DSC  and  for  which  AH0  is 
known)  are  presented  in  Section  S3-3.  In  Section  S3-4,  the  method  is  used  to  determine  the  AH0  value 
for  ESR-255,  which  cannot  be  determined  by  DSC  alone  as  noted  in  Section  S2. 

S3.2  Method.  Three  or  four  samples  cured  under  different  conditions  in  an  oven  so  as  to  span  a  range 
of  conversions  from  not  far  beyond  the  vitrification  point  to  near  full  conversion  were  prepared.  An 
ATR-FTIR  spectrum  of  a  small  piece  of  each  sample  was  collected  using  a  Thermo  Corporation  Nicolet 
6700  FT-IR  Spectrometer  in  attenuated  total  reflectance  mode  from  400  -  4000  cm"1  with  a  resolution  of 
4  cm"1.  Spectra  were  determined  by  averaging  32  scans.  After  being  scanned,  a  small  bit  of  each  piece 
(~  5  mg)  was  sealed  in  a  DSC  pan  and  heated  to  350  °C  at  10  °C/min,  followed  by  cooling  at  10  °C/min.  to 
100  °C  and  re-heating  to  350  °C  at  10  °C/min  in  a  TA  Instruments  200  DSC  under  50  mL/min.  of  flowing 
nitrogen.  The  DSC  heat  flow  was  determined  at  each  temperature  using  the  second  heating  scan  as  a 
baseline,  and  then  integrated  to  determine  a  residual  enthalpy  of  cure  (AHr).  In  order  to  estimate  the 
error  in  AHr,  a  baseline  was  selected  manually  based  on  the  shape  of  the  first  heating  curve  only  and 
used  to  perform  a  separate  integration.  The  difference  between  the  two  integrated  values  for  the  same 
curve  was  then  recorded  as  a  rough  error  estimate. 
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A  linear  baseline  connecting  the  points  at  2000  cm'1  and  2500  cm'1  was  chosen  for  each  FT-IR 
spectra.  The  absorbance  values  were  then  normalized  such  that  the  peak  absorbance  near  1500  cm'1 
(representing  phenyl  group  vibrations  and  thus  insensitive  to  extent  of  cure),  after  subtracting  the 
baseline,  was  set  to  an  arbitrary  value  of  1.  The  normalized  absorbance  values  in  each  spectra  (with  the 
baseline  subtracted)  were  then  integrated  (using  a  dx  of  1.9285  cm'1  =  an  arbitrary  value  of  1  for 
convenience  based  on  the  format  of  instrument  output)  between  the  apparent  edges  of  the  cyanate 
ester  peak  (generally  2100  cm'1  to  2400  cm'1,  but  for  the  re-analysis  of  older  data  for  BADCy  and  LECy  a 
range  of  2100  cm'1  to  2300  cm'1  was  used  to  avoid  a  noisy  region  in  the  spectrum).  In  order  to  estimate 
an  error  for  this  integrated  value,  alternative  baseline  points  of  1800  cm"1  and  2400  cm'1  were  chosen. 
These  points  represent  a  sort  of  "least  reasonable"  baseline;  in  contrast  to  the  "most  reasonable" 
baseline  points  of  2000  cm'1  and  2500  cm'1,  these  points  represent  what  we  believe  were  the  furthest 
limits  away  from  the  peak  that  one  could  choose  as  endpoints  for  a  baseline  without  being  completely 
"unreasonable".  The  cyanate  ester  was  then  re-integrated  (without  re-normalizing  the  spectra)  using 
the  alternate  baseline,  and  the  difference  between  the  results,  which  represents  a  measure  of  the 
sensitivity  of  the  integration  to  the  choice  of  baseline,  was  used  as  a  rough  estimate  of  the  error  in  the 
measurement.  It  should  be  pointed  out  that  the  potential  baseline  error  was  judged  to  be  by  far  the 
largest  source  of  error  in  the  FT-IR  measurements. 

Once  a  complete  set  of  FT-IR  peak  areas  and  AHr  values,  with  the  respective  error  estimates, 
was  obtained,  the  FT-IR  peak  area  was  plotted  as  a  function  of  AHr.  A  linear  regression  of  the  data  was 
then  performed  to  determine  the  x-intercept  of  the  best  fit  line.  In  some  cases,  the  errors  were  many 
times  larger  for  a  single  data  point  than  for  the  others;  in  such  cases,  the  data  point  with  the  large 
estimated  error  was  simply  excluded  from  the  regression  (we  did  not  believe  a  weighted  regression  was 
justified  because  the  error  values  were  only  rough  estimates  of  variance).  The  error  in  the  x-intercept 
was  then  estimated  in  two  ways:  1)  based  on  the  standard  error  of  the  regression  extrapolated  to  the  x- 
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intercept,  and  2)  based  on  a  sensitivity  analysis  to  the  estimated  variations  in  each  data  point,  choosing 


a  worst-case  scenario  in  which  the  variations  in  each  x-  and  y-  value  all  affected  the  x-intercept  in  a 
similar  manner  (i.e.  assuming,  for  instance,  all  values  were  one  estimated  y-error  lower  and  one 
estimated  x-error  higher  than  measured,  so  that  all  variations  increased  the  x-intercept  value).  The 
linearity  of  the  plots  in  general  was  also  noted.  If  the  assumptions  explained  in  Sections  S3.1  are  correct 
(namely,  FT-IR  and  DSC  signals  that  are  linearly  proportional  to  cyanate  ester  concentration  over  the 
conversion  range  from  post-vitrification  to  full  cure),  then  the  data  points  on  the  plot  as  constructed 
should  fall  on  a  straight  line. 

The  x-intercept  value  was  interpreted  as  follows:  a  positive  intercept  indicates  that  some 
residual  heating  would  be  observed  even  if  all  cyanate  ester  groups  were  consumed,  hence  it  would 
indicate  that  the  DSC  is  detecting  chemistry  other  than  cure  (network  thermal  degradation,  for 
instance).  A  negative  x-intercept  would  indicate  that  even  in  a  sample  with  no  residual  cure  detected  in 
the  DSC,  cyanate  ester  peaks  would  remain  in  the  FT-IR  scan,  and  hence,  that  the  residual  heating  scan 
to  350  °C  in  the  DSC  was  not  sufficient  to  achieve  full  conversion.  The  value  of  the  x-intercept  would 
then  represent  either  the  enthalpy  of  side  reactions  (for  positive  values)  to  be  subtracted  from  the  DSC 
integration  to  obtain  the  true  residual  enthalpy  of  cure,  or  (for  negative  values),  the  enthalpy  to  be 
added  to  the  DSC  integration  to  obtain  the  true  enthalpy  of  cure  of  all  remaining  cyanate  ester  groups  (if 
in  fact  they  could,  hypothetically,  be  cured  completely  without  side  reactions).  In  equation  form, 

AHr  =  AHr,DSC-AHxi  (S3. 2) 

Where  AFirDSc  is  the  apparent  residual  enthalpy  measured  by  DSC  and  AHxi  is  the  x-intercept  obtained  by 
linear  regression  of  a  plot  of  cyanate  ester  FT-IR  peak  intensity  (y)  vs.  apparent  DSC  residual  enthalpy  (x) 
for  a  series  of  calibration  samples.  It  should  be  noted  that  under  most  circumstances,  AHxi  would  be 
expected  to  show  no  significant  differences  from  zero.  Note  that  this  formula  also  applies  to  the 
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determination  of  AH0  by  DSC,  which  may  be  considered  a  special  case  of  the  determination  of  AHr  using 


a  monomeric  sample.  The  value  of  AHxi  is  thus  usefully  compared  to  the  estimated  error  in  the 
measurement  of  AH0  by  DSC. 

S3.3  Method  Validation.  For  validation  of  this  method,  we  utilized  previously  published  data  on  the 
commercial  monomers  BADCy  and  LECy  as-supplied  from  Lonza  and  cured  at  various  temperatures  for 
12  hours  in  order  to  produce  a  range  of  conversions.52  The  systems  were  catalyzed  with  2  phr 
nonylphenol  and  160  ppm  Cu  (present  as  copper(ll)  acetylacetonate).  In  this  validation  case,  the  DSC 
residual  heating  produced  a  closed  curve,  hence  the  absolute  conversions  were  available  from  the  DSC 
data  alone.  Thus,  the  combined  DSC/IR  technique  should  result  in  an  adjustment  to  the  residual  heat  of 
cure  that  is  not  statistically  significant. 

In  Figure  S3-1,  the  residual  heat  measured  by  DSC  (values  published  in  Reference  S2)  is  plotted 
against  the  FT-IR  peak  intensity  integrated  from  2100  -  2300  cm"1.  For  reference,  the  IR  spectra  are 
shown  in  Figure  S3-2  (these  data  were  discussed  in  Reference  S2  but  peak  integration  as  described  in 
Section  S3. 2  had  not  been  performed).  A  close  look  at  the  spectra  in  Figure  S3-2  shows  that,  particularly 
at  a  final  cure  temperature  of  150  °C,  the  slope  of  the  baseline  is  somewhat  uncertain.  As  a  result,  the 
characteristic  uncertainty  estimate  depicted  in  Figure  S3-1  is  comparatively  large.  It  is  also  easy  to  see 
why  such  estimates  of  uncertainty  are  important,  because  even  a  small  change  in  the  assumed  baseline 
makes  a  large  difference  in  the  integrated  intensity.  A  linear  regression  of  the  points  in  Figure  S3-1 
yields  an  x-intercept  of  -35  J/g,  both  with  and  without  the  150  °C  cure  point  included.  Based  on  the 
range  of  uncertainties  reported  and  the  small  number  of  data  points,  one  would  expect  the  true  error  in 
the  x-intercept  estimate  to  be  at  least  35  J/g,  so  that  the  result  is  not  significantly  different  than  zero. 
Figures  S3-3  and  S3-4  present  the  analogous  data  for  catalyzed  LECy.  As  with  the  BADCy,  the  baseline 
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Residual  Heat  (J/g) 


Figure  S3-1.  DSC-IR  linear  regression  of  catalyzed  BADCy,  showing  the  relationship  between  residual 
heat  seen  in  the  DSC  scan  and  integrated  FT-IR  peak  intensity.  The  expected  x-intercept  is  zero. 


Figure  S3-2.  Comparative  FT-IR  spectra  of  catalyzed  BADCy  as  a  function  of  final  cure  temperature. 


15 


Normalized  Absorbance 


Figure  S3-3.  DSC-IR  linear  regression  of  catalyzed  LECy,  showing  the  relationship  between  residual  heat 
seen  in  the  DSC  scan  and  integrated  FT-IR  peak  intensity.  The  expected  x-intercept  is  zero. 


Figure  S3-4.  Comparative  FT-IR  spectra  of  catalyzed  LECy  as  a  function  of  final  cure  temperature. 
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of  the  150  °C  cured  sample  showed  considerable  uncertainty.  The  x-intercept  of  the  DSC/IR  comparison 


plot  is  an  insignificant  5  J/g  with  the  150  °C  data  included  and  29  J/g  without  it.  Again,  given  the 
uncertainties  and  the  small  amount  of  data  collected,  the  x-intercept  does  not  appear  significantly 
different  than  zero. 

Taken  together,  the  validation  cases  show  that,  in  general,  the  expected  relationship  between 
DSC  and  FT-IR  data  does  exist.  However,  these  data  clearly  show  the  importance  of  having  a  precisely 
determined  baseline.  They  also  indicate  that  in  order  to  obtain  a  more  reliable  linear  fit  of  the  data,  four 
or  five  data  points  will  be  preferable  to  only  three,  since  a  substantial  fraction  of  the  data  may  suffer 
from  significant  uncertainties  associated  with  the  assumed  baselines.  Having  four  or  five  data  points 
would  allow  for  meaningful  weighted  regressions  (having  only  three  data  points  means  that  in  many 
cases,  weighted  regressions  will  effectively  include  just  two  data  points). 

S3.4  Results  for  ESR-255.  Having  shown  with  closed-curved  DSC  data  that  (with  a  margin  of  error  of  a 
few  per  cent)  the  DSC/IR  technique  enables  the  determination  of  absolute  conversion,  the  technique 
was  applied  to  ESR-255.  Figure  S3-5  shows  the  DSC/IR  correlation  plot  while  the  individual  FT-IR  spectra 
are  shown  in  Figure  S3-6.  In  this  case,  the  samples  analyzed  were  all  uncatalyzed  ESR-255  cured  for  30, 
60,  120,  and  240  minutes  at  250  °C  in  an  oven.  As  a  further  check  on  the  method,  we  carried  out  an 
isothermal  cure  at  250  °C  for  240  minutes  in  the  DSC  on  a  small  sample  of  uncatalyzed  ESR-255  powder, 
and  then  subjected  this  in-situ  DSC  cured  sample  to  a  residual  heating.  While  the  residual  heating  after 
240  minutes  is  computed  as  normal,  the  residual  heats  at  30,  60,  and  120  minutes  were  computed  by 
adding  the  peak  integral  from  the  residual  heating  scan  to  the  integrated  heat  flow  between  the  time  of 
interest  (e.g.  30  minutes)  and  the  end  of  isothermal  cure  at  240  minutes.  These  two  methods  (oven  and 
DSC)  should  give  equal  residual  heating  values.  We  performed  this  particular  experiment  simply  as  a 
check  on  this  newly  developed  method;  we  would  not  anticipate  performing  such  a  check  routinely. 
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Figure  S3-5.  DSC/IR  correlation  plot  for  uncatalyzed  ESR-255  using  conventional  oven  cured  samples  as 
well  as  a  sample  cured  in-situ  in  the  DSC  as  a  special  check  of  the  method. 
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Figure  S3-6.  Individual  FT-IR  spectra  of  uncatalyzed  ESR-255  cured  in  an  oven  at  250  °C  for  the  times 
indicated. 
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As  seen  in  Fig  S3-5,  the  oven  cured  samples  tended  to  show  slightly  lower  residual  heats  than 


comparable  samples  cured  in  the  DSC.  This  could  be  due  to  a  larger  expected  temperature  overshoot  in 
the  oven  as  compared  to  the  precisely  controlled  DSC.  The  x-intercept  values  were  -129  ±  28  J/g  for  in- 
situ  DSC  cured  samples  and  102  ±  30  J/g  for  oven-cured  samples.  Note  that  the  error  estimates  are  from 
the  equally  weighted  linear  regressions  and  do  not  take  into  account  the  error  in  the  individual 
measurements,  which  is  of  a  similar  magnitude.  Hence,  the  total  uncertainty  in  the  x-intercept  is  likely 
around  50  J/g.  The  difference  between  the  methods  is  not  significant,  however,  the  intercept  itself 
does  appear  to  be  significantly  different  from  zero  at  about  90%  confidence.  To  obtain  the  correct  AH0, 
we  thus  subtracted  the  x-intercept  value  of  the  oven-cured  samples  (since  those  samples  were  used  for 
FT-IR),  or  -100  J/g  from  the  value  obtained  by  DSC  (730  J/g)  to  yield  a  AH0  or  830  J/g  or  (keeping  only 
two  significant  figures,  110  kJ/mol),  with  an  uncertainty  of  about  7  kJ/mol.  It  is  worth  noting  that  110 
kJ/mol  is  in  fact  the  expected  value  for  cyanate  esters,  and  matches  the  range  reported  by  us  using  DSC 
only  for  uncatalyzed  BADCy  and  LECy. 
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S4.  Determination  of  diBenedetto  Equation  Parameters 


S4.1  Introduction.  For  thermosetting  polymers  with  a  well-defined  network  structure  for  a 
given  conversion,  such  as  cyanate  esters,  the  diBenedetto  equation  provides  a  key  link  among  physical 
properties  by  relating  the  glass  transition  temperature  (Tg)  of  the  network  to  the  conversion  (a).sl  The 
equation  takes  the  form: 

(Tg~Tgo)  _ 

(Tgoo- T^o)  l— (l— A)cr 

in  which  Tg0  and  Tg«  represent  the  glass  transition  temperatures  of  the  monomer  and  fully  cured 
network,  respectively,  and  in  which  A  is  effectively  an  adjustable  parameter  (0  <  A  <  1)  that  in  theory  is 
determined  by  the  heat  capacity  characteristics  of  the  monomer  and  cured  network.sl  According  to  eq. 
S4-1,  the  glass  transition  temperature  of  a  cyanate  ester  network  is  independent  of  its  cure  history,  a 
fact  that  has  been  verified  experimentally  for  at  least  two  different  cyanate  ester  monomers.S3S4  As  a 
result,  for  any  cure  process,  be  it  a  complex  part  cured  in  an  oven  or  a  simple  non-isothermal  DSC  scan, 
the  network  Tg  at  any  location  and/or  instant  in  time  may  be  computed  as  long  as  the  conversion  at  the 
same  location  and/or  instant  is  known.  For  carefully  performed  DSC  scans,  an  instantaneous  measure  of 
conversion  is  always  available  by  integrating  the  heat  flow  and  dividing  by  AH0.  Thus,  even  in  systems 
such  as  ESR-255,  the  developing  Tg  may  be  determined  once  a  technique  such  as  the  combined  DSC/IR 
method  described  in  Section  S3  has  been  used  to  correctly  measure  AFIq  and  the  diBenedetto  equation 
parameters  have  been  computed. 

Techniques  for  the  determination  of  diBenedetto  equation  parameters  for  thermosetting 
polymers  may  be  divided  into  two  general  cases.  In  the  first  case,  which  applies  to  most  lower  service 
temperature  cyanate  esters  such  as  BADCy,  and  also  applies  to  monomer  1  (that  is,  l,3,5-tris(4- 
cyanatophenylmethyljbenzene),  the  parameters  Tg0  and  Tg*o  are  easily  determined  by  a  method  such  as 
DSC,  and  only  the  parameter  A  must  be  determined  experimentally.  In  such  cases,  even  a  single  Tg 
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measurement  in  combination  with  a  conversion  measurement  on  a  partly  cured  system  allows  for  a 


determination  of  X,  with  more  than  one  point  needed  for  an  estimate  of  uncertainty.  In  practice, 
however,  the  values  of  X  that  we  have  determined  for  a  variety  of  cyanate  esters  always  fall  in  the  range 
of  0.25  to  0.45,  and  regardless  of  how  many  data  points  are  used,  it  appears  that  an  uncertainty  of  0.05 
to  0.1  is  commonly  encountered  due  primarily  to  uncertainties  in  measuring  both  Tg  and  a.  Moreover, 
in  many  cases  a  fairly  narrow  range  of  conversions  (beyond  the  vitrification  point)  is  actually  of  practical 
interest,  thus  there  is  not  much  to  be  gained  from  performing  many  measurements  of  Tg  at  low 
conversion.  Instead,  a  few  Tg  and  a  pairs  at  high  conversion  may  be  used  to  determine  X  through  either 
a  non-linear  curve  fitting  approach,  or,  the  following  convenient  form  may  be  used  in  a  linear  regression 
(with  intercept  fixed  at  zero) 


S4-2 


For  cyanate  esters  with  high  glass  transition  temperatures,  however,  Tg«  may  not  be  directly 
measureable,  because  the  network  may  undergo  thermal  degradation  before  the  thermodynamic 
transitions  associated  with  Tgoo  can  be  observed  in  an  actual  heating  experiment.  As  a  result,  multiple 
parameters  appearing  in  the  diBenedetto  equation  may  need  to  be  determined  simultaneously.  In  such 
cases,  having  many  measurements  available  is  advisable  since  the  optimization  needed  is  both  multi¬ 
variate  and  nonlinear.  At  present,  the  best  approach  to  performing  this  kind  of  optimization  with 
experimental  data  on  thermosetting  polymers  remains  an  open  question.  As  in  the  simpler  problem, 
data  at  lower  conversions  does  not  add  much  to  the  robustness  of  the  optimization.  Thus,  when 
collecting  data,  if  the  conversions  are  not  controlled  (as  with  a  scheduled  cure  experiment),  it  may  be 
necessary  to  collect  a  large  number  of  data  points  just  to  ensure  that  enough  of  these  points  are 
available  at  high  conversions  to  enable  a  useful  determination. 
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S4.2  Methods  for  Determination  of  diBenedetto  Parameters.  Because  the  available  quantities  of  the 


three  monomers  varied  (2  was  produced  only  in  low  yield  after  many  synthesis  steps,  while  ESR-255  was 
produced  in  high  yield  from  one  step),  and  because  the  three  monomers  differed  in  availability  of  fully- 
cured  glass  transition  temperature  data,  we  chose  to  estimate  diBenedetto  parameters  based  on 
separate  approaches  judged  best  for  each  monomer.  The  level  of  information  obtained  thus  differed 
quite  a  bit  from  monomer  to  monomer,  with  most  parameters  for  2  available  as  fairly  broad  ranges 
while  parameters  for  ESR-255  were  all  available  as  point  estimates  with  clearly  determined 
uncertainties. 

S4.2.1  for  monomer  1  (l,3,5-tris(4-cyanatophenylmethyl)benzene).  The  glass  transition  temperature 
of  the  monomer  was  determined  by  DSC  on  heating  from  -70  °C  to  100  °C  at  10  °C  /  min.  The  glass 
transition  temperature  of  the  fully  cured  polymer  was  also  determined  DSC  scans  on  the  second  heating 
scan  of  a  cycle  from  100  °C  to  350  °C  at  10  °C  /  min.  All  DSC  scans  were  performed  on  a  TA  Instruments 
Q200  under  50  mL/min  of  flowing  nitrogen.  To  determine  a  value  for  X,  the  DSC  scans  of  samples 
previously  heated  to  250  °C,  270  °C,  and  290  °C  for  30  minutes  were  utilized.  These  samples  were 
heated  to  350  °C  at  10  °C  /  min  as  part  of  the  standard  procedure  for  chemical  kinetics  analysis  (see 
Section  S6),  a  technique  which  simultaneously  measured  the  conversion  of  the  as-cured  samples.  The 
observed  heating  scans  showed  no  residual  cure  below  the  glass  transition  temperature,  so  no 
adjustment  for  in-situ  cure  was  needed.  The  DSC  scans  of  samples  isothermally  cured  at  lower 
temperatures  showed  residual  cure  but  no  apparent  glass  transition  temperature  on  residual  heating. 
These  scans  were  used  to  exclude  a  range  of  possible  values  for  X  on  the  basis  that  a  glass  transition 
would  have  been  visible  had  X  fallen  within  the  excluded  range.  For  instance,  when  heated  to  230  °C 
for  30  minutes,  a  conversion  of  around  0.82  for  monomer  1  resulted.  At  a  conversion  of  0.82,  a  value  of 
0.25  for  X  would  have  resulted  in  a  glass  transition  temperature  of  159  °C  for  the  partly  cured  1,  whereas 
a  value  of  0.35  for  X  would  have  resulted  in  a  glass  transition  temperature  of  187  °C  (based  on  the 
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known  Tg0  and  Tgoo  values,  see  Table  S4-1).  Since  a  Tg  was  definitely  not  observed  between  153  °C  and 
178  °C  in  the  residual  heating  scan  under  these  conditions,  a  value  of  0.25  for  X  may  be  excluded, 
whereas  a  value  of  0.35  for  X  cannot  be  excluded  (on  the  basis  of  this  particular  scan). 

The  individual  estimates  for  X  based  on  single  heating  scans  (found  by  solving  the  diBenedetto 
equation  algebraically)  were  compared  with  the  best  fit  value  obtained  by  a  weighted  least  squares 
method  in  which  X  was  varied  in  increments  of  0.01,  then  in  increments  in  0.001  to  refine  the  estimate. 
The  weighting  factor  for  each  point  was  determined  by  the  inverse  square  of  the  estimated  uncertainty 
in  the  experimentally  observed  Tg  value.  For  further  comparison,  a  linear  regression  of  the  data  using 
Eq.  S4-2  was  also  performed.  Lastly,  the  values  of  X  excluded  on  the  basis  of  a  lack  of  observed  glass 
transition  temperature  signals  in  DSC  scans  were  checked  against  the  other  estimates  for  consistency. 

S4.2.2  for  monomer  2  (3,5-Bis[(4-cyanatophenyl)methyl)]phenylcyanate).  The  glass  transition 
temperature  of  the  uncured  monomer  was  determined  using  the  technique  described  in  S4.2.1.  The 
fully  cured  glass  transition  temperature  was  apparent  neither  from  DSC  nor  from  TMA  measurements, 
however,  because  2  is  a  close  structural  analog  of  REX-371  (whose  fully  cured  Tg  was  reported  at  370 
°C)S1,  as  well  as  a  small,  model  oligomer  of  PT-30  (fully  cured  Tg  reported  at  400  °C)S5  the  fully  cured  Tg 
may  be  assigned  a  probable  value  of  385  ±  15°C.  Based  on  this  information,  constraints  on  the  value  of 
X  could  be  determined  from  the  residual  heating  portion  of  isothermal  DSC  cure  experiments  in  which 
no  glass  transition  temperature  was  observed,  as  described  in  S4.2.1,  but  no  glass  transition 
temperatures  of  partly  cured  samples  were  available  for  an  estimate  for  X.  Instead,  the  values  of  X  for 
the  analogous  monomers  1  and  ESR-255  were  simply  averaged  to  determine  a  rough  estimate,  and  this 
estimate  was  checked  for  consistency  with  the  constraints  obtained  from  the  DSC  data.  Due  to  the 
many  uncertainties  in  the  estimation  of  diBenedetto  parameters  for  this  monomer  and  its  reliance  on 
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analogy  with  other  known  systems,  we  have  noted  in  the  main  manuscript  that  estimates  of 


diBenedetto  parameters  for  2  are  not  of  the  same  quality  as  those  for  1  and  ESR-255. 

S4.2.3  for  ESR-255  (l,l,l-tris(4-cyanatophenyl)ethane).  Due  to  the  large  available  quantity  of  this 
monomer,  a  different  technique  for  the  estimation  of  diBenedetto  parameters  were  employed.  The 
three  diBenedetto  equation  parameters  Tg0,  Tgoo,  and  X  were  estimated  simultaneously  using  the  Solver 
nonlinear  optimization  routine  (with  default  settings)  found  in  Microsoft®  Excel  2010.  Since  a  detailed 
study  of  the  relative  merits  of  various  nonlinear  techniques  was  judged  to  be  well  beyond  the  required 
scope  of  effort,  and  since  values  for  Tg0  and  X  were  easy  to  check  for  reasonableness,  we  judged  that 
simply  adapting  the  generic  nonlinear  optimization  algorithm  present  in  Solver  was  the  best  choice. 

In  order  to  generate  data  on  partly  cured  systems,  samples  were  heated  at  10  °C  /  min  using  the 
DSC  method  analogous  to  that  described  in  S4.2.1  to  a  chosen  isothermal  temperature  and  cured  for  30 
minutes,  then  cooled  at  10  °C  /  min  to  50  °C.  After  five  minutes  of  isothermal  quiescence  (to  remove 
transients  from  the  DSC  signal),  the  samples  were  subjected  to  modulated  DSC  by  heating  at  5  °C  /  min 
to  350  °C  to  measure  residual  cure  (based  on  the  non-reversing  portion  of  the  signal).  A  modulation 
amplitude  of  1  °C  and  a  modulation  period  of  60  s  were  used,  along  with  a  re-scanned  baseline  for  both 
isothermal  and  modulated  residual  heating.  All  AH  values  were  calibrated  using  the  enthalpy  of  melting 
of  ESR-255  as  an  internal  standard,  as  described  in  Section  S6.  Note  that  because  samples  did  not  cure 
fully  under  these  conditions,  the  conversion  was  determined  based  on  the  measurement  of  the  heat 
evolved  during  the  isothermal  portion  of  the  cure,  using  AH0  =  110  kJ/eq.  (based  on  the  results 
described  in  Section  S3).  In  other  words,  the  formula  used  to  determine  conversion  was 


a  =  AH  iso  /  AH0. 


S4-3 
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As  a  check  on  consistency  of  the  data,  we  did  calculate  the  residual  heat  AHr  based  on  AHr,DSCand  AHxi 
according  to  Eq.  S3-2,  then  checked  that  the  traditional  Eq.  S3-1  gave  the  same  conversion  values  as 
those  obtained  via  Eq.  S4-3. 

The  isothermal  30-minute  temperatures  used  to  determine  values  of  X,  Tg0,  and  Tgo<,  were  150, 
220,  230,  240,  245,  250,  255,  260,  270,  and  290  °C.  For  each  data  point,  the  value  of  Tg  seen  during  the 
residual  cure  scan  (from  the  reversible  portion  of  the  signal)  was  compared  to  the  value  predicted  by  the 
diBenedetto  equation,  and  the  parameters  of  the  diBenedetto  equation  that  minimized  the  sum  of 
squared  differences  between  predicted  and  observed  values  were  determined  using  the  Microsoft® 
Excel  Solver  algorithm.  Note  that  the  lower  temperatures  tended  to  give  low  conversions  (and  thus 
values  of  Tg  not  that  different  from  Tg0,  while  temperatures  of  270  and  290  °C  did  not  produce  a 
measurable  Tg.,  hence  8  data  points  were  used  to  find  the  values  for  the  three  unknown  variables.  Tg0 
was  also  determined  using  a  separate  DSC  scan  using  the  procedure  described  in  S4.2.1  and  compared 
to  the  value  found  by  the  Solver  algorithm. 

S4.3  Results  and  Discussion 

S4.3.1  for  monomer  1  (l,3,5-tris(4-cyanatophenylmethyl)benzene).  The  value  of  Tg0  found  by  DSC  was 

-16  °C.  Table  S4-1  provides  the  observed  actual  and  minimum  possible  observed  Tg  values  for  the  partly 
cured  samples  as  a  function  of  cure  temperature,  based  on  the  reasoning  described  in  S4.2.1,  along  with 
the  resulting  point  estimate  or  range  estimate  for  A,  and  the  observed  value  of  Tg„.  The  estimated  value 
of  V  was  taken  as  the  average  of  these  values,  or  320  ±  4  °C.  The  value  of  X  determined  by  the  non¬ 
linear  least  squares  was  0.300,  in  line  with  the  average  value  of  0.31  ±  0.07  obtained  from  the  three 
point  estimates.  The  linear  regression  using  Eq.  S4-2  is  shown  in  Figure  S4-1.  The  value  of  X  was 
estimated  at  0.29  using  a  weighted  regression  with  weights  proportional  to  the  inverse  of  the 
uncertainty  in  the  Tg  measurement  squared.  An  unweighted  regression  produced  an  estimate  of  0.28  ± 
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0.02  for  X.  Note  that  these  values  are  all  consistent  with  the  constraints  arising  from  the  DSC  scans  in 


which  no  discernible  Tgwas  observed.  Taking  all  the  methods  into  consideration,  a  value  of  0.30  ±  0.05 
for  X  seemed  reasonable.  When  considering  the  uncertainties  in  the  model,  it  should  be  remembered 
that  the  uncertainties  in  many  of  the  input  parameters  (e.g.  DSC  measurements)  are  not  based  on  well- 
studied  statistical  distributions,  and  that  the  errors  in  the  model  parameters  are  correlated  (e.g.  an 
overestimate  of  Tgoo  will  produce  an  underestimate  in  X),  thus  standard  propagation  of  error  calculations 
are  not  likely  to  be  useful.  Instead,  more  direct  methods  can  be  utilized.  For  instance,  the  error  in 
predicting  the  Tgof  partly  cured  systems  may  be  estimated  simply  by  examining  the  residuals  in  the  best 
fit  model,  which  in  this  case  revealed  an  average  deviation  of  just  5  °C. 


Table  S4-1 


Estimated  diBenedetto  equation  parameters  from  DSC  scans  of  monomer  1 


Cure  Temp.  (°C) 

a 

Tg  as-cured  (  C) 

Ts~  (°C) 

X 

230 

0.82 

none  (158-168)* 

320  ±5 

<  0.25  or  >  0.28 

250 

0.92 

246  ±2 

321  ±2 

0.29  ±0.02 

270 

0.96 

273  ±8 

318  ±2 

0.23  ±0.08 

290 

0.98 

304  ±2 

319  ±2 

0.38  ±0.09 

*  the  numbers  indicate  the  temperature  range  over  which  a  Tg  was  clearly  not  observed 
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Figure  S4.1  Regression  plot  using  the  linearized  form  of  the  diBenedetto  equation  (eq.  S4-2) 
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S4.3.2  for  monomer  2  (3,5-Bis[(4-cyanatophenyl)methyl)]phenylcyanate).  The  value  of  Tg0  found  by 
DSC  was  -23  °C,  similar  to  that  found  for  monomer  1.  As  mentioned  in  S4.2.2,  the  values  of  Tgoo  and  X 
were  estimated  simply  by  analogy  with  data  on  comparable  cyanate  esters  as  385  °C  and  0.31, 
respectively.  The  constrains  on  the  values  of  X  found  by  DSC  scans  are  listed  in  Table  S4-2,  and,  in  total, 
they  reveal  only  that  values  of  X  below  about  0.29  are  unlikely.  A  reasonable  upper  bound  for  X  based 
on  comparable  cyanate  esters  would  be  about  0.4.  Hence,  the  uncertainty  in  X  is  asymmetric.  Since 
there  are  no  partly  cured  Tg  values  to  provide  a  direct  estimate  of  the  prediction  error,  we  note  that  if 
the  model  predictions  (listed  in  the  Table)  were  typically  more  than  about  5  °C  too  high,  a  Tg  would  have 
been  visible  where  none  was  observed.  Based  on  the  difference  in  model  predictions  between  X  =  0.3 
and  X=  0.4,  an  underestimate  of  about  20  °C  on  average  was  possible.  The  uncertainty  in  Tg»o  also  results 
in  a  change  of  plus  or  minus  10  °C  on  average.  If,  for  lack  of  a  better  simple  alternative,  the  errors  are 
assumed  to  combine  randomly,  then  the  total  uncertainty  ranges  from  an  overestimate  of  about  12  °C 
to  an  underestimate  of  about  22  °C.  Thus,  the  quality  of  these  estimates  is  lower  than  for  the  other 
monomers,  however  the  values  are  still  quite  reasonable  to  employ  in  studies  of  the  evolution  of  the 
glass  transition  temperature  during  cure  so  long  as  the  errors  and  assumptions  are  noted.  One  such 
assumption  is  that  the  system  cured  fully  when  heated  to  350  °C  in  the  DSC,  which,  given  the  expected 
value  of  Tgo<„  may  not  be  true.  However,  a  major  overestimation  of  Tg  due  to  invalidity  of  this 
assumption  is  unlikely,  as  there  is  a  good  chance  it  would  have  produced  at  least  one  clearly  seen  Tg 
value  where  in  fact  none  was  seen. 


Table  S4-2 


Estimated  diBenedetto  equation  parameters  from  DSC  scans  of  monomer  2 


Cure  Temp.  (°C) 

a 

Tg  as-cured-  (  C) 

T  /°r\** 

1  g  model  V  W 

X 

230 

0.78 

none  (155-185)* 

190 

<0.22  or  >0.29 

250 

0.90 

none  (155-244)* 

279 

<0.08  or  >0.21 

270 

0.95 

none  (155-265)* 

326 

<0.04  or  >0.12 

290 

0.98 

none  (155-255)* 

360 

<0.02  or  >0.04 

*  the  numbers  indicate  the  temperature  range  over  which  a  Tg  was  clearly  not  observed 
**  based  on  TgM  =  385  °C  and  A  =  0.31. 
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S4.3.3  for  ESR-255.  Table  S4-3  shows  the  observed  and  predicted  (based  on  the  best  fit  model)  Tg 
values  for  partly  cured  ESR-255  as  a  function  of  isothermal  cure  temperature.  The  best  fit  model 
parameters  obtained  by  the  Solver  algorithm  were  Tg0  =  -9  °C,  Tg^  =  558  °C,  and  X  =  0.32.  Figure  S4-2 
illustrates  the  agreement  between  the  fitted  model  and  the  data.  Note  that  Tg0  was  separately 
measured  at  8  ±  3  °C,  hence  the  fitted  model  tends  to  underestimate  the  Tg  at  low  conversions. 
Although  the  value  of  X  is  similar  to  what  was  found  for  monomer  1,  the  value  of  Tgo<,  seems  remarkably 
high.  By  heating  a  disc  sample  to  decomposition  at  50  °C/min  in  a  TMA  under  nitrogen,  a  loss  peak  was 
observed  at  419  °C  (see  Figure  S5-7).  Such  an  experiment  should  be  attempted  only  with  caution  and 
with  strong  measures  to  contain  the  gaseous  decomposition  products.  The  results  indicate  that  the  Tgof 
the  fully  cured  network  in  the  absence  of  thermal  degradation,  which  is  what  the  diBenedetto 
parameter  represents,  must  be  at  least  419  °C.  Indeed,  considering  that  the  removal  of  just  one  flexible 
methylene  group  from  monomer  1  increased  the  Tg«  of  the  network  by  50  -  80  °C,  a  further  increase  of 
160  -  190  °C  caused  by  the  removal  of  two  more  such  groups  (which  yields  ESR-255)  is  not 
unreasonable.  The  average  deviation  between  the  predicted  and  observed  glass  transition 
temperatures  using  these  model  parameters  was  10  °C,  hence,  for  Tg  values  below  about  300  °C,  the 
model  provides  a  good  empirical  predictive  capability. 


Table  S4-3 


Estimated  diBenedetto  equation  parameters  from  DSC  scans  of  ESR-255 


Cure  Temp.  (°C) 

a 

Tg  as-cured  (  C) 

Tg  predicted  (  C) 

150 

0 

8  ±  3 

-10 

220 

0.17 

22 

25 

230 

0.19 

24 

30 

240 

0.25 

34 

45 

245 

0.34 

62 

71 

250 

0.50 

122 

129 

255 

0.60 

198 

174 

260 

0.80 

304 

310 
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Figure  S4.2.  Glass  transition  temperature  as  a  function  of  conversion  for  ESR-255 
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S5.  TMA  Data  on  "As  Cured",  "Fully  Cured",  and  "Wet"  Samples 


Figures  S5-1  through  S5-3  show  the  dynamic  TMA  (storage  and  loss  components  of  stiffness 
along  with  their  ratio  (tan  delta),  for  monomer  1,  monomer  2,  and  ESR-255,  respectively.  Note  that  the 
heating  rate  used  in  these  experiments  was  a  very  rapid  50  °C  /  min.,  a  rate  that  is  necessary  to  avoid 
significant  in-situ  cure  of  the  samples  prior  to  reaching  the  glass  transition  temperature.  Above  the  glass 
transition  temperature,  in-situ  cure  is  unavoidable,  and  its  effects  are  evident  in  the  TMA  data  plots. 
For  monomers  1  and  2  the  effects  are  subtle,  such  as  asymmetric  peak  shapes  and  a  less  than  expected 
drop  in  the  storage  component  of  stiffness  on  passing  through  Tg,  indicating  that  in-situ  cure  is  not 
severe.  For  ESR-255,  however,  significant  in-situ  cure  begins  right  at  the  glass  transition  temperature, 
owing  to  the  comparatively  low  conversion.  The  in-situ  cure  causes  the  Tg  to  rise  faster  than  the  sample 
heating  rate,  which  leads  to  a  reversal  of  the  signal  trend,  that  is,  the  storage  actually  increases  rapidly, 
and  the  peak  in  tan  delta  coincides  with,  rather  than  lags,  the  peak  in  the  loss  component  (a  lagging  tan 
delta  indicates  normal  de-vitrification,  a  leading  tan  delta  indicates  vitrification,  and,  typically,  a 
coincident  tan  delta  indicates  de-vitrification  followed  almost  immediately  by  re-vitrification). 
Remarkably,  the  glass  transition  temperature  of  ESR-255  in  the  solid  state  increased  faster  than  50  °C  / 
min.  Through  the  diBenedetto  equation,  one  may  calculate  that  a  conversion  rate  of  5  x  10'4  /  s  is 
sufficient  to  raise  Tg  at  the  required  rate,  a  rate  about  seven  times  higher  than  the  observed  rate  seen  in 
the  kinetic  data  for  ESR-255  at  250  °C  and  a  similar  conversion. 
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Figure  S5-1.  Dynamic  TMA  of  "as  cured"  monomer  1. 
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Figure  S5-2.  Dynamic  TMA  of  "as  cured"  monomer  2. 
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Figure  S5-3.  Dynamic  TMA  of  "as  cured"  ESR-255. 
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Figures  S5-4  through  S5-6  present  the  corresponding  dynamic  TMA  data  for  the  "fully  cured' 


monomer  1,  monomer  2,  and  ESR-255,  respectively.  Note  that  full  cure  is  achieved  by  first  heating  these 
samples  at  a  much  slower  10  °C  /  min  to  350  °C,  in  effect  facilitating  rather  than  preventing,  in-situ  cure. 
Only  in  the  case  of  monomer  1  is  a  Tg  visible,  at  338  °C  (loss  peak).  Figure  S5-7  shows  the  dynamic  TMA 
of  ESR-255  heated  to  decomposition.  When  the  storage  finally  begins  to  drop,  with  a  loss  peak  at  419 
°C,  decomposition  is  already  well  under  way.  (Please  note:  such  an  experiment  should  be  done  only 
with  caution  and  careful  provision  for  the  safe  and  rapid  removal  of  decomposition  products.)  Thus,  it  is 
not  possible  to  determine  the  (theoretical)  softening  point  of  the  network  in  the  absence  of 
degradation;  only  a  lower  limit  for  the  softening  point  was  found.  Although  the  softening  point  of  ESR- 
255  has  been  reported  previously,56  the  extent  of  cure  in  the  samples  was  not  simultaneously  reported. 

Figures  S5-8  through  S5-10  show  the  dynamic  TMA  (first  heating)  of  the  "wet"  monomer  1, 
monomer  2,  and  ESR-255,  respectively.  Note  that  "wet"  refers  to  having  been  immersed  in  85  °C  water 
for  96  hours.  As  is  typical  for  "wet"  samples,  bubbles  formed  soon  after  the  Tg  was  reached,  causing 
significant  disturbances.  Nonetheless,  the  actual  glass  transition  was  visible  just  prior  to  the  onset  of  the 
disturbances. 
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Figure  S5-4.  Dynamic  TMA  of  "fully  cured"  monomer  1. 
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Figure  S5-5.  Dynamic  TMA  of  "fully  cured"  monomer  2. 
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Figure  S5-6.  Dynamic  TMA  of  "fully  cured"  ESR-255. 


Figure  S5-7.  Dynamic  TMA  of  ESR-255  heated  to  decomposition. 
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Figure  S5-8.  Dynamic  TMA  of  "wet"  monomer  1. 
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Figure  S5-9.  Dynamic  TMA  of  "wet"  monomer  2. 
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Figure  S5-10.  Dynamic  TMA  of  "wet"  ESR-255. 
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S6.  Detailed  Information  on  the  Investigation  of  Cure  Kinetics 


56.1  Introduction.  The  analysis  of  cure  kinetics  generally  followed  the  pattern  used  in  our  previous 
work;S7  with  some  notable  extensions  in  order  to  handle  more  complex  data.  For  those  readers  familiar 
with  our  previous  efforts.  Section  S6.2  summarizes  these  efforts  in  terms  of  the  minor  changes  and 
expansions  from  the  previous  work.  Section  S6.3  details  the  treatment  of  the  raw  DSC  data  to  extract 
the  rate  of  conversion  as  a  function  of  time  including  baselines  (S6.3.1),  isothermal  data  (S6.3.2), 
estimation  of  residual  cure  (S6.3.3)  and  combined  error  analysis  (S6.3.4).  Section  S6.4  describes  the 
analytical  methods  employed  to  determine  the  rate  constants,  for  both  individual  curves  (S6.4.2)  and 
sets  of  data  with  temperature-independent  exponents  (S6.4.3).  In  Section  S6.5,  alternative  analyses  are 
utilized  in  order  to  provide  for  an  estimation  of  uncertainty,  since  traditional  propagation  of  error 
techniques  are  too  unwieldy  to  be  practical  for  such  a  complex,  nonlinear  multi-variate  problem.  The 
final  data  sets  are  all  tabulated  in  Section  S6.6.  Note  that  kinetic  parameters  were  calculated  only  for 
samples  that  reached  a  clearly  defined  maximum  in  their  cure  rate  within  the  allotted  time.  In  the  case 
of  ESR-255,  for  which  only  the  highest  temperatures  met  the  above  criterion,  additional  data  sets  with 
longer  cure  times  (4  hours  at  250  °C,  from  the  DSC/IR  studies  reported  in  Section  S3,  and  24  hours  at  210 
°C,  from  the  studies  of  purity  reported  in  Section  S8),  were  used  to  supplement  the  data. 

56.2  Summary  of  Changes  Relative  to  Method  Published  in  Guenthner  et  al.  (Polymer,  2011,  52,  3933- 
3942). 

1.  A  two-level  nested  iteration  is  used  to  find  k/,  k2',  m,  and  n,  relaxing  the  assumption  that  k/  « 
k2'  used  previously.  The  starting  point  for  the  iteration  is  the  same,  namely  the  assumption  that 
m  /  (m  +  n)  =  amax,  but  now  the  quantity  m  /  (m  +  n)  is  adjusted  iteratively  to  match  the  model 
otmax  (which  differs  slightly  from  the  value  of  m  /  (m  +  n)  due  to  k/  being  non-negligible)  to  the 
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experimental  amax.  This  expansion  allows  for  determination  of  model  parameters  when  k/  is 


much  larger  compared  to  k2'. 

2.  The  conversion  ranges  used  in  the  model  are  now  based  on  the  observed  maxima  (and  in  some 
cases  minima)  in  the  conversion  rate  as  a  function  of  conversion  curve,  rather  than  on  fixed 
values  of  conversion.  Note  that  these  new  criteria  are  designed  to  match  the  previous,  fixed 
values  in  cases  where  K/-  «  k2'  and  the  maximum  rate  of  conversion  is  near  a  =  0.5,  as  was  true 
in  all  samples  previously  analyzed,  but  is  not  true  in  all  of  the  cases  reported  herein.  Hence,  the 
new  criteria  allow  for  consistent  analysis  across  a  wider  range  of  parameters,  particularly  it 
allows  for  consistent  analysis  when  k/-  is  non-negligible  and  m  and  n  are  well  above  one. 

3.  The  baselines  for  determining  the  rate  of  conversion  are  now  determined  by  re-scanning  the 
sample  through  the  entire  procedure  once  residual  cure  is  complete.  This  change  was 
implemented  in  order  to  match  methods  previously  reported  by  others,  and  its  validity  was 
confirmed  by  the  analysis  of  baseline  effects  reported  in  Section  S6.3.1.  Note  that  for  non- 
isothermal  segments,  the  re-scanned  data  is  offset  as  described  in  S6.3.1  to  construct  the 
baseline. 

Note  that  the  previously  reported  kinetic  data  can  thus  be  regarded  as  simply  a  special  case  of  the 
current  analysis  technique,  one  in  which  k/  «  k2-,  amax  =  0.5,  and  the  re-scanned  and  internal  baselines 
are  not  significantly  different. 

S6.3  Treatment  of  DSC  Data 

S6.3.1.  Baselines.  As  noted  in  S6.2,  the  DSC  data  were  converted  to  conversion  rates  using  re-scanned 
baselines,  that  is,  the  entire  DSC  method  was  repeated  after  the  completion  of  cure  in  order  to  generate 
the  baseline  values.  Since  this  technique  has  been  shown  to  work  well  by  others/3  we  investigated  the 
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effects  of  using  a  re-scanned  baseline  in  order  to  validate  (and  modify,  if  necessary)  the  technique,  as 


well  as  to  understand  the  measurement  uncertainties  related  to  the  choice  of  baseline. 

We  identified  two  main  concerns  with  the  use  of  re-scanned  baselines.  The  first,  referred  to  as 
baseline  drift,  is  that,  in  a  method  that  takes  a  few  hours  to  run,  a  slow  drift  of  the  baseline  over  time 
may  lead  to  significant  error.  To  investigate  this  effect,  we  examined  the  heat  flow  values  measured 
during  the  initial,  five-minute  isothermal  segment  of  the  DSC  method  at  120  °C,  a  temperature  at  which 
the  sample  was  completely  melted  and  no  cure  took  place.  To  minimize  the  effect  of  transients,  only 
the  last  minute  of  the  segment  was  used.  Table  S6-1  shows  the  differences  in  the  heat  flow  values 
recorded  on  the  initial  scan  compared  to  those  recorded  on  the  re-scan.  A  baseline  drift  would  be 
expected  to  produce  differences  that  varied  in  a  random  way  among  samples.  Another  possible  effect, 
baseline  shift,  however  would  be  expected  to  produce  differences  that  varied  in  systematic  ways.  The 
most  likely  cause  of  what  we  term  "baseline  shift"  is  solidification  of  the  sample.  For  the  lower  cure 
temperatures,  the  sample  does  not  solidify  during  the  isothermal  cure  segment,  but  does  solidify  during 
residual  cure.  For  the  higher  temperatures,  the  sample  solidifies  during  the  isothermal  cure  segment.  If 
a  systematic  baseline  drift  took  place,  it  would  be  present  at  about  the  same  value  for  both  the  120  °C 
comparison  and  the  "end  of  isothermal  cure"  comparison,  for  all  samples.  Flowever,  if  a  shift  due  to 
solidification  were  present,  then  for  samples  cured  at  higher  temperatures,  the  difference  should  be 
larger  in  the  120  °C  comparison  than  in  the  "end  of  isothermal  cure"  comparison,  since  solidification  will 
have  taken  place  prior  to  both  of  the  measurements  used  in  the  "end  of  isothermal  cure"  comparison. 

To  help  the  reader  identify  the  relevant  data,  we  have  marked  samples  that  solidify  between  the 
initial  scan  and  the  re-scan  with  an  "x"  in  Table  S6-1.  When  an  "x"  is  present  in  the  first  but  not  the 
second  column,  then  a  systematic  drift  should  result  in  similar  values  in  both  columns,  whereas  a 
solidification-induced  shift  should  result  in  smaller  values  in  the  second  column.  Table  S6-1  also  includes 
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the  projected  "worst  case"  difference  in  measured  conversions  after  the  isothermal  cure  segment,  as 


well  as  in  the  maximum  conversion  rate,  that  result  from  assuming  that  the  baseline  drift  or  shift,  as 
indicated  by  the  differences  in  the  120  °C  segment,  either  are  or  are  not  valid  corrections  to  the  data. 
Note  that,  during  the  isothermal  segments,  changes  in  sample  heat  capacity  should  not  contribute  to 
differences  in  the  measured  heat  flow. 


Table  S6-1 

Parameters  Related  to  Investigation  of  Re-Scanned  DSC  Baseline  (Isothermal  Segments) 


Sample  ID 

Difference  in 

heat  flow 

(initial 

isothermal 

segment), 

mWa 

Difference  in  heat 
flow  (end  of 
isothermal  cure 
segment),  mWa 

Difference  in 

conversion  after 

isothermal  cure 
(drift/shift  vs.  no 
drift/shift 
assumed)15 

Relative  change  in 
max  conversion 
rate  (drift/shift  vs. 
no  drift/shift 
assumed)3 

Monomer  1-210  °C 

-0.0161  x 

n/ac  x 

-1.8% 

-1.0x10 5  /s  (-3.6%) 

Monomer  1  -  230  °C 

-0.0129  x 

n/ac  x 

-1.1% 

-5.9x10 6  /s  (-0.7%) 

Monomer  1  -  250  °C 

-0.0108  x 

-0.020 

-1.2% 

-6.4x10 6  /s  (-0.3%) 

Monomer  1  -  270  °C 

-0.0124  x 

-0.011 

-1.2% 

-6.6x10 6  /s  (-0.2%) 

Monomer  1  -  290  °C 

-0.0128  x 

-0.007 

-1.2% 

-6.6x10 6  /s  (-0.1%) 

Monomer  2  -  210  °C 

-0.0099  x 

n/ac  x 

-0.7% 

-3.8x10 6  /s  (-5.1%) 

Monomer  2  -  230  °C 

-0.0172  x 

n/ac  x 

-2.4% 

-1.3x10 5  /s  (-2.4%) 

Monomer  2  -  250  °C 

-0.0168  x 

-0.039 

-1.5% 

-8.5x10 6  /s  (-0.7%) 

Monomer  2  -  270  °C 

-0.0309  x 

-0.027 

-2.3% 

-1.3x10 5  /s  (-0.5%) 

Monomer  2  -  290  °C 

-0.0349  x 

-0.024 

-2.4% 

-1.3x10 5  /s  (-0.2%) 

ESR-255  -  230  °C 

-0.0116  x 

n/ac  x 

-0.7% 

-4.1x10 6  /s  (-4.3%) 

ESR-255  -  230  °C 

-0.0033  x 

n/ac  x 

-0.3% 

-1.7x10 6  /s  (-1.0%) 

ESR-255  -  250  °C 

-0.0184  x 

n/ac  x 

-1.1% 

-6.3x10 6  /s  (-1.3%) 

ESR-255  -  270  °C 

-0.0428  x 

-0.040 

-3.3% 

-1.8x10 5  /s  (-1.7%) 

ESR-255  -  290  °C 

-0.0238  x 

-0.086 

-2.1% 

-1.1x10  s  /s  (-0.5%) 

ESR-255  -  250  °C  4h 

-0.0253  x 

-0.007 

-16.7% 

-1.2x10 5  /s  (-2.8%) 

a  value  on  re-scan  minus  value  on  first  scan 

b  (value  assuming  drift/shift  minus  value  assuming  no  drift  /  shift)  /  value  assuming  no  drift/shift 
c  heat  flow  signal  did  not  reach  a  constant  value  during  isothermal  cure  (incomplete  cure) 

"x"  indicates  sample  vitrified  (i.e.  Tg  surpassed  the  isothermal  temperature)  between  measurements 

The  differences  seen  in  Table  S6-1  do  seem  to  show  systematic  effects.  In  particular,  the 
samples  cured  at  the  higher  temperatures  (or  the  longest  times),  which  result  in  the  highest  extent  of 
solidification,  show  the  expected  effects  of  a  baseline  shift  due  to  solidification.  Note  that  the  sample 
cure  rates  at  the  end  of  the  isothermal  segment,  which  are  not  quite  equal  to  zero,  may  also  affect  the 
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results.  The  ESR-255  sample  that  was  cured  isothermally  for  four  hours  as  part  of  the  DSC/IR 


experiments  described  in  Section  S3  also  provided  a  useful  set  of  data.  If  baseline  drift  were  a  major 
issue,  one  would  expect  this  sample  to  show  a  significantly  larger  difference  than  the  others,  whereas  in 
fact  the  differences  are  quite  similar.  Thus,  while  baseline  drift  can  be  ruled  out  as  a  source  of  concern, 
baseline  shift  due  to  solidification  cannot  be  ruled  out. 

Although  it  cannot  be  ruled  out  as  a  source  of  concern,  Table  S6-1  shows  that,  for  30  minute 
isothermal  cure  segments,  the  baseline  shift  issue  is  not  a  major  source  of  error.  The  projected 
difference  in  conversion  results  from  assuming  a  constant  error  in  the  heat  flow  rate,  but,  in  the  case  of 
a  shift  due  to  solidification,  such  a  difference  will  only  affect  the  data  up  to  the  point  of  solidification. 
For  samples  cured  at  the  higher  temperatures,  this  point  is  reached  very  early  in  the  isothermal  cure 
segment,  such  that  the  probable  error  from  neglecting  to  account  for  baseline  shift  is  about  an  order  of 
magnitude  less  than  the  worst-case  estimate  shown  in  Table  S6-1.  For  the  ESR-255  /  four  hour  cure 
sample,  similarly,  the  large  error  in  conversion  will  only  occur  if  the  sample  does  not  solidify  during  the 
entire  four  hour  isothermal  cure.  Since,  in  fact,  the  sample  solidifies  after  only  30-60  minutes,  the  actual 
error  in  conversion  is  more  likely  to  be  just  2-4%.  As  will  be  shown  in  discussion  of  errors  (S6.3.4),  the 
projected  errors  in  Table  S6-1  are  generally  a  small  component  of  the  errors  present  in  the 
measurement.  Thus,  although  baseline  shift  due  to  solidification  cannot  be  ruled  out,  it  can  be  safely 
ignored  without  significant  loss  of  data  quality,  so  long  as  either  1)  isothermal  cure  times  are  limited  to 
30  minutes,  or  2)  solidification  of  the  sample  takes  place  in  about  30  minutes  if  a  longer  isothermal  time 
is  used. 


For  non-isothermal  segments,  not  only  are  baseline  drifts  and  shifts  due  to  solidification 
possible,  but  also  changes  in  sample  mass  or  heat  capacity  can  also  contribute  to  apparent  offsets  in  the 
re-scanned  baseline.  As  with  the  isothermal  data,  the  effects  of  drift  should  be  essentially  independent 
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of  sample  condition,  whereas  baseline  shifts  should  be  more  prominent  in  samples  that  solidify  between 


the  initial  scan  and  the  re-scan.  Though  mass  loss  is  expected  to  be  negligible  in  the  sealed  DSC 
samples,  heat  capacity  changes  might  also  be  expected  to  play  a  role  in  samples  that  exhibit  significant 
residual  cure.  Table  S6-2  lists  the  apparent  offset  in  the  baseline  (determined  by  considering  the 
maximum  amount  the  re-scanned  heat  flow  would  need  to  be  offset  to  avoid  negative  values  of  residual 
cure,  as  illustrated  in  Figure  S6-1)  for  the  different  samples  studied.  Samples  that  solidify  between  the 
initial  scan  and  re-scan  are  indicated  by  an  "x"  in  Table  S6-2  as  in  Table  S6-1.  In  most  cases,  samples 
that  have  solidified  prior  to  the  initial  scan  show  much  smaller  offsets  than  do  those  that  solidify 
between  the  scans. 

The  effect  of  ignoring  the  apparent  offset  in  the  baseline  on  conversion  is  also  shown  in 
parenthesis  in  the  entries  in  Table  S6-2.  Unlike  the  case  for  the  isothermal  segments,  the  apparent 
offsets  do  make  a  significant  difference  in  the  calculation  of  conversions,  potentially  up  to  about  10%. 
Also,  unlike  the  isothermal  cases,  if  a  baseline  shift  does  occur  due  primarily  to  solidification,  for  these 
samples  it  generally  would  have  taken  place  at  least  halfway  through  the  non-isothermal  scan 
(vitrification  is  approximately  synchronous  with  the  peak  of  the  residual  cure  for  samples  not  already 
solidified),  and  thus  most  of  the  potential  error  (which  is  based  on  the  assumption  of  an  incorrect  heat 
flow  throughout  the  entire  segment)  will  likely  be  realized.  Hence,  these  errors  are  too  large  to  ignore, 
and,  therefore,  in  the  treatment  of  the  non-isothermal  data  we  have  assumed  that  a  baseline  shift  does 
occur  and  made  use  of  the  apparent  offset  in  calculating  the  residual  enthalpy  of  cure. 
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Table  S6-2 


Apparent  Offset  in  Re-scanned  Baseline  (Non-isothermal  Segments)  and  Its  Effect  on  Conversion 

Estimates 


Cure  Temperature  (°C) 

Monomer  1 

Monomer  2 

ESR-255 

210 

+0.073  (5.7%)  x 

-0.213  (9.6%)  x 

-0.233  (10.0%)  x 

230 

-0.010  (0.6%)  x 

-0.016  (1.0%)  x 

-0.145  (9.2%)  x 

250 

-0.004  (0.3%) 

+0.049  (2.9%) 

-0.276  (11.5%)  x 

270 

+0.003  (0.2%) 

+0.047  (2.2%) 

-0.090  (4.6%) 

290 

-0.005  (0.3%) 

+0.088  (4.0%) 

-0.038  (2.2%) 

"x"  indicates  the  sample  vitrified  (i.e.  Tg  surpassed  the  cure  temperature)  between  measurements 


Temperature  (°C) 


Figure  S6-1.  Illustration  (using  residual  scan  of  ESR-255  cured  at  250  °C  for  30  minutes)  of  the  apparent 
offset  needed  when  constructing  a  baseline  from  heat  flow  data  generated  during  a  re-scan  of  the 
sample  (non-isothermal  segments  only).  The  offset  is  calculated  by  subtracting  the  re-scanned  data 
from  the  initial  data,  noting  the  minimum  (most  negative  or  least  positive)  value. 


S6.3.2  Issues  with  isothermal  data.  The  primary  issue  affecting  the  isothermal  data  is  the 
determination  of  conversion  rates  at  the  start  of  the  isothermal  segment.  Although  heating  rates 
approaching  100  °C  /  min.  are  attainable  in  the  DSC,  this  means  that  the  "jump"  needed  to  reach  the 
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isothermal  cure  temperature  can  take  up  to  about  two  minutes,  followed  by  a  short  period  of 


equilibration  during  which  transient  signals  will  affect  the  data.  Particularly  during  the  last  30  seconds 
or  so  of  the  "jump"  for  samples  cured  at  high  temperatures,  a  significant  amount  of  cure  does  take  place 
and  cannot  be  ignored  without  significant  loss  of  data  quality.  Fortunately,  the  use  of  the  re-scanned 
baseline  provides  a  convenient  way  to  estimate  the  heat  flow  during  the  initial  "jump"  and  transient 
periods.  However,  because  the  "jump"  involves  very  high  heating  rates,  the  DSC  signal  during  the  jump 
will  be  dominated  by  the  heat  capacity  term,  which  is  proportional  to  the  heating  rate.  Hence,  during 
the  "jump"  and  initial  transient  period,  which  we  defined  as  any  period  in  which  the  heating  rate 
exceeded  0.1  °C  /  min.,  the  baseline  was  constructed  by  finding  the  re-scanned  data  points  with  the 
matching  heating  rate,  rather  than  the  matching  elapsed  time  as  was  done  for  the  remaining  isothermal 
data. 


The  use  of  the  re-scanned  baseline  in  the  manner  just  depicted  can  lead  to  slightly  negative 
apparent  values  for  the  rate  of  cure.  In  order  to  correct  this  problem,  we  identified  the  latest  point  in 
the  data  (always  well  before  the  end  of  the  "jump")  in  which  an  apparent  cure  rate  of  zero  or  less  was 
computed.  We  then  set  the  cure  rate  for  all  times  up  to  this  point  to  zero,  essentially  excluding  it  from 
consideration.  Furthermore,  the  dynamics  of  the  transients  are  likely  not  reproduced  exactly  during  the 
re-scan,  though,  typically  these  transients  produce  sinusoidal  patterns  superimposed  on  the  data,  and 
therefore,  if  averaged  over  time  (as  in  integration  to  determine  conversion)  the  transient  effects  are 
significantly  reduced.  At  lower  cure  rates,  the  transient  effects  are  important  for  a  longer  time  period. 
Hence,  we  believe  the  best  method  for  dealing  with  the  transients  is  to  exclude  from  further  analysis  the 
cure  rate  data  for  the  first  10%  of  conversion  beyond  the  initial  minimum  seen  in  the  rate  (often  at 
zero),  with  the  exception  of  integrating  the  heat  flow  values  in  order  to  estimate  the  conversion  at  the 
start  of  the  useable  rate  data.  Using  a  cut-off  criteria  based  on  conversion  automatically  compensates 
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for  the  aforementioned  effects  of  temperature  on  the  length  of  time  that  transients  are  important,  but 


also  necessitates  the  integration  of  the  heat  flow  during  the  transient  period. 

As  mentioned  in  the  discussion  on  baselines  (S6.3.1),  for  the  isothermal  data,  any  apparent 
baseline  shifts  due  to  solidification  of  the  sample  are  not  taken  into  account.  In  addition  to  subtraction 
of  baseline,  there  are  issues  relating  to  how  to  calibrate  (or  scale)  the  heat  flow  data.  The  weighing  of 
mg  scale  DSC  samples  typically  involves  errors  of  at  least  a  few  per  cent,  and  for  isothermal  cure  studies 
on  cyanate  esters  it  is  important  not  to  use  large  samples  due  to  the  need  to  rapidly  remove  the 
considerable  quantities  of  heat  generated  during  rapid  cure.  Fortunately,  for  the  highly  crystalline 
monomers  under  study,  the  enthalpy  of  melting  can  be  measured  to  within  about  1%,  and,  since  it  was 
measured  many  times  (see  the  tables  in  Section  S7),  it  was  used  as  an  internal  standard  to  correct  for 
small  weighing  errors. 

S6.3.3  Issues  with  non-isothermal  data.  For  the  non-isothermal  data,  the  only  analysis  performed  is 
an  integration  of  the  residual  heat  flow  to  determine  the  residual  enthalpy  of  cure  and,  hence, 
conversion  at  the  end  of  the  isothermal  cure  segment.  For  this  data  collection,  the  initial  transient 
involves  quenching  the  sample  and  holding  at  low  temperatures.  We  did  not  attempt  to  integrate  the 
heat  flows  during  the  quench,  since,  for  samples  cured  at  high  temperature,  significant  cure  tended  to 
cease  long  before  the  end  of  the  isothermal  period,  while,  at  low  cure  temperatures,  significant  cure 
ceased  almost  immediately  once  the  quench  began.  Since,  as  mentioned  in  S6.3.2,  the  main  issue  with 
transients  is  loss  of  accurate  heat  flow  data  rather  than  integrated  errors,  and  since  the  only  operation 
performed  on  non-isothermal  data  was  integration,  we  did  not  investigate  transient  effects  in  detail.  In 
determining  the  offset  needed  to  construct  the  baseline,  there  is,  as  a  consequence  of  the  method  used, 
a  single  point  at  which  the  heat  flow  equals  zero  (see  Figure  S6-1).  We  simply  assumed  that  this  point 
represented  the  start  of  residual  cure  and  did  not  integrate  any  signals  prior  to  this  point  during  the 
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segment.  Since  this  point  typically  was  well  past  the  start  of  the  segment  (and  typically  became  later 


and  later  in  the  segment  with  increasing  cure  temperature),  it  effectively  excludes  portions  of  the  data 
affected  by  transients  from  consideration. 

Note  that,  as  explained  in  Section  S6.3.1,  apparent  baseline  shifts  are  not  ignored  in  the  non- 
isothermal  data,  due  to  the  magnitude  of  the  potential  errors  involved.  Including  the  shifts  does  seem 
appropriate  as,  at  lower  isothermal  cure  temperatures,  it  tends  to  result  in  a  lower  apparent  onset 
temperature  for  residual  cure,  while  at  higher  temperatures,  it  tends  to  result  in  a  higher  apparent  onset 
temperature  for  residual  cure.  In  addition,  as  with  the  isothermal  cure  data,  the  measured  enthalpy  of 
melting  was  used  as  an  internal  standard  to  calibrate  the  heat  flow  values  as  explained  at  the  end  of 
Section  S6.3.2. 

S6.3.4  Error  estimates.  Since,  as  outlined  above,  there  are  significant  potential  issues  with  the  handling 
of  the  raw  DSC  data,  we  believe  that  it  is  important  to  attempt  to  understand  the  potential  errors 
involved.  In  Tables  S6-1  and  S6-2,  we  have  included  "worst-case"  error  estimates  for  the  effects  on 
measured  conversions  at  the  end  of  the  isothermal  cure  step  and  maximum  rates  of  cure  during  the 
isothermal  step.  The  effect  of  our  treatment  of  transients,  though,  is  much  more  difficult  to  assess, 
other  than  to  note  that,  based  on  the  appearance  of  what  seem  to  be  artificial  "bumps"  in  the  otherwise 
smooth  conversion  rate  curves  (see  Figures  S6-10  through  S6-12,  and  note  that  these  data  were 
collected  and  shown  but  not  used  in  subsequent  analysis),  the  effect  of  transients  appears  to  extend 
consistently  through  the  first  4-7%  of  conversion.  The  magnitude  of  these  "bumps"  would  be  consistent 
with  errors  in  conversion  rates  of  around  10%,  and  hence  total  errors  in  integrated  conversion  of  less 
than  1%. 

Although  we  were  able  to  obtain  some  rough  estimates  of  the  effects  of  individual  sources  of 
error,  many  of  these  errors  are  not  likely  to  involve  normal  probability  distributions.  Baseline  shifts,  for 
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instance,  either  do  or  do  not  occur,  and  although  their  estimated  effects  involve  some  random  errors, 


the  actual  error  in  the  baseline  estimate  is  more  likely  to  be  10%  of  the  computed  offset  (for  a  properly 
corrected  shift)  or  90%  of  the  computed  offset  (for  an  improperly  corrected  shift)  than  say,  50%  of  the 
computed  offset,  resulting  in  a  complex  bimodal  probability  distribution.  Therefore,  propagation  of 
error  calculations  cannot  be  performed  in  a  straightforward  way;  a  better  approach  to  estimating  the 
overall  error  is  to  look  at  the  variability  in  composite  measurements. 

Fortunately,  the  total  enthalpy  of  cure  for  each  monomer  should  be  independent  of  cure 
temperature.  Even  for  ESR-255,  which  cannot  be  fully  cured  in  the  DSC  (see  Sections  S2  and  S3),  the 
enthalpy  of  cure  up  to  350  °C  is  still  likely  to  be  independent  of  the  earlier  isothermal  cure  temperature, 
provided  that  the  isothermal  cure  temperatures  used  are  well  below  350  °C.  Even  better,  for  cyanate 
esters,  the  total  enthalpy  of  cure  is  almost  always  in  the  range  of  100  -  110  kJ/eq.  for  uncatalyzed 
samples,  and  is  available  from  non-isothermal  measurements,  thus,  an  assessment  of  accuracy  as  well 
as  precision  is  available  by  comparing  the  total  enthalpy  of  cure  for  these  samples.  Again,  even  for  ESR- 
255,  the  extent  of  cure  for  ramps  to  350  °C  is  known  from  the  IR  data,  and  can  be  combined  with  the 
known  enthalpy  of  cure  of  cyanate  esters  to  provide  an  assessment  of  accuracy.  Lastly,  if  we  wish  to 
determine  whether  or  not  a  particular  factor  contributes  significantly  to  the  total  error,  we  can  in  many 
cases  check  the  effect  on  the  variation  in  the  total  enthalpy  of  cure.  If  a  particular  change  (such  as 
ignoring  a  baseline  shift)  does  not  affect  the  variation  in  the  total  enthalpy  of  cure,  then  it  is  likely  not  a 
significant  source  of  error.  In  addition  to  composite  measurements,  in  some  cases  we  have  data  sets 
that  are,  practically  speaking,  replicates.  For  instance,  the  250  °C  isothermal  cure  data  for  ESR-255  that 
was  collected  as  part  of  the  IR/DSC  determination  of  conversion  (Section  S3)  can  be  compared  to  the 
data  generated  using  the  standard  cure  kinetics  procedure  at  the  same  temperature  to  get  a  sense  of 
the  level  of  reproducibility  of  the  data. 
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Table  S6-3  shows  the  total  enthalpy  of  cure  for  each  of  the  samples  analyzed  for  cure  kinetics, 


plus  the  replicated  ESR-255  data  at  250  °C.  The  first  numerical  column  shows  the  values  for  the  data  set 
as  treated;  subsequent  columns  show  the  effect  of  possible  variations,  namely,  including  a  correction  for 
baseline  shift  in  the  isothermal  data,  not  offsetting  the  re-scanned  baseline  in  the  non-isothermal  data, 
and  not  making  use  of  the  heat  of  melting  as  an  internal  calibration,  respectively.  For  each  monomer 
tested,  the  average  and  standard  deviation  are  shown.  The  average  values  are  also  compared  to  those 
found  from  non-isothermal  DSC  runs  (the  standard  method  for  estimating  enthalpies  of  cure).  For  ESR- 


Table  S6-3 

Enthalpies  of  cure  of  cyanate  ester  monomers  for  various  DSC  data  treatment  scenarios 


Sample  ID 

AHcure  as- 

measured 

(kJ/eq.) 

AHcure  -  if  baseline 
shift  in  isothermal 
segments  is 
included  (kJ/eq.) 

AHcure  -  if  offset  in 
non-isothermal 
segment  baseline 
is  not  included 
(kJ/eq.) 

AHcure-  if  internal 
calibration  is 
ignored  (kJ/eq.) 

Monomer  1-210  °C 

118.9 

116.7 

125.7 

94.0 

Monomer  1  -  230  °C 

123.0 

121.7 

123.7 

154.5 

Monomer  1  -  250  °C 

127.1 

125.6 

126.7 

123.1 

Monomer  1  -  270  °C 

128.7 

127.2 

129.0 

124.0 

Monomer  1  -  290  °C 

127.8 

126.3 

127.4 

122.2 

Average  (Isothermal  runs) 

125.1  ±4.1 

123.5  ±4.3 

126.5  ±2.0 

123.6  ±21.5 

Non-isothermal  DSC 

119.3 

n/a 

106.6 

122.7 

Monomer  2  -  210  °C 

127.1 

126.3 

115.0 

130.5 

Monomer  2  -  230  °C 

123.2 

120.2 

121.9 

117.9 

Monomer  2  -  250  °C 

119.4 

117.6 

122.9 

125.1 

Monomer  2  -  270  °C 

113.1 

110.5 

115.6 

117.3 

Monomer  2  -  290  °C 

128.0 

124.9 

133.1 

129.1 

Average  (Isothermal  runs) 

122.1  ±6.1 

119.9  ±6.3 

121.7  ±7.3 

124.0  ±6.2 

Non-isothermal  DSC 

105.7 

n/a 

145.6 

100.4 

ESR-255  -  230  °C 

102.6 

101.8 

92.3 

101.1) 

ESR-255  -  230  °C 

100.2 

99.9 

91.0 

109.4 

ESR-255  -  250  °C 

99.4 

98.3 

88.0 

91.4 

ESR-255  -  270  °C 

95.2 

92.0 

90.8 

85.2 

ESR-255  -  290  °C 

97.8 

95.7 

95.7 

105.7 

Average  (Isothermal  runs) 

99. 1±  2.8 

97.6  ±3.8 

91.6  ±2.8 

97.9  ±  11.5 

Non-isothermal  DSC 

93.1 

n/a 

80.1 

92.6 

ESR-255 -250  °C  4h  (IR) 

95.5 

79.5 

84.4 

85.0 

*  from  the  data  set  used  to  determine  diBenedetto  equation  parameters 
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255,  the  enthalpy  shown  is  for  cure  to  350  °C  rather  than  complete  cure  (the  DSC/IR  technique  in 


Section  S3  indicates  that  the  total  enthalpy  of  cure  should  be  about  14  kJ/eq.  higher  than  the  listed 
enthalpies). 

The  key  result  indicated  by  the  data  in  Table  S6-3  is  that,  of  the  four  methods,  the  as-treated 
data  provided  the  least  variation  between  isothermal  and  non-isothermal  methods,  the  least  variation 
among  the  different  monomers,  which  should  have  similar  enthalpies  of  cure,  and  either  the  lowest  or 
second  to  lowest  variation  within  each  set  of  isothermal  kinetics  data  for  a  given  monomer.  Of  the  19 
samples  investigated,  the  only  sample  for  which  an  alternate  method  appears  to  give  a  more  reliable 
estimate  of  the  enthalpy  of  cure  is  monomer  1  cured  isothermally  at  210  °C.  This  particular  sample  was 
notable  (see  Table  S6-2)  for  having  the  largest  positive  offset  in  the  non-isothermal  baseline  of  any  of 
the  samples  listed  in  the  table,  an  offset  that  is  an  order  of  magnitude  larger  than  any  of  the  others  for 
monomer  1.  This  sample  features  the  shallowest  quench  depth  of  any  sample  (87  °C).  Moreover,  in  this 
sample,  the  conditions  at  the  end  of  isothermal  cure  (conversion  near  50%  and  hence  auto-catalysis 
having  maximum  influence)  are  the  most  difficult  for  readily  quenching  the  reaction.  It  thus  appears 
that  in  this  particular  case,  complete  quenching  was  not  achieved,  meaning  that  the  minimum  observed 
DSC  signal  with  respect  to  the  re-scanned  value  did  not  reflect  true  baseline  (no  residual  cure) 
conditions,  since  cure  was  likely  not  completely  stopped  at  the  start  of  the  residual  heating  scan.  As  a 
result,  the  computed  offset  for  the  baseline  was  artificially  positive,  leading  to  a  residual  enthalpy  of 
cure  that  was  too  low  and  an  error  of  approximately  5%  in  the  total  enthalpy  of  cure.  Since  this 
particular  shortcoming  in  the  as-used  baseline  construction  method  affected  only  1  of  19  samples,  was 
associated  with  clear  indicators  of  anomalous  conditions,  and  resulted  in  an  error  that  was  similar  in 
magnitude  to  the  apparent  uncertainty  of  the  measurement  in  general,  we  believe  it  is  more 
appropriate  to  retain  an  identical  data  treatment  method  for  all  samples  rather  than  to  employ  complex 
criteria  in  order  to  eliminate  this  occasional  source  of  error. 
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From  Table  S6-3  it  can  also  be  seen  that  the  isothermal  methods  (which  include  a  non- 


isothermal  segment  with  the  same  temperature  ramp  rate  and  range  as  the  purely  non-isothermal 
methods)  yield  slightly  larger  values  of  the  enthalpy  of  cure  using  the  as-treated  data.  This  result  makes 
sense  given  that  the  isothermal  method  always  involves  additional  opportunities  for  the  sample  to  cure 
compared  to  the  purely  non-isothermal  data,  although  the  same  result  could  also  be  due  to  systematic 
errors  such  as  an  overestimation  of  heat  flow  when  transient  signals  affect  the  isothermal  data.  In 
contrast,  any  data  sets  that  show  a  lower  enthalpy  of  cure  for  the  isothermal  method  most  likely  involve 
errors  in  data  treatment.  Such  data  sets  are,  however,  only  generated  when  the  offset  in  the  non- 
isothermal  segment  baseline  is  ignored.  Further,  for  the  purely  non-isothermal  methods,  ignoring  the 
offset  in  the  baseline  results  in  large  differences  in  cure  enthalpies  for  the  three  monomers  studied 
(even  accounting  for  the  incomplete  cure  of  ESR-255),  a  situation  which  is  highly  unlikely  in  reality. 
Ignoring  the  internal  calibration  of  heat  flows  by  measurement  of  the  enthalpy  of  melting  at  times 
results  in  large  variations  within  monomers,  but,  on  average,  yields  results  very  similar  to  the  other 
methods.  This  feature  is  a  good  indicator  that  the  errors  introduced  in  weighing  the  samples  (that  are 
corrected  via  the  internal  calibration)  are  indeed  unbiased,  and  hence,  in  non-crystalline  monomers  that 
do  not  afford  such  internal  calibration,  replication  of  samples  can  be  used  effectively  to  reduce  this 
source  of  error,  which  appears  to  average  around  10%. 

In  summary  the  total  enthalpy  of  cure  data  indicate  a  precision  of  around  3-5%  in 
measurements  obtained  directly  from  DSC  data  without  the  use  of  kinetic  models.  Comparisons  among 
different  monomers  also  indicate  a  relative  error  of  around  5%  in  the  measured  data,  while  comparisons 
of  the  measured  enthalpies  of  cure  to  the  expected  value  of  110  kJ/eq.  indicate  that,  in  terms  of 
absolute  accuracy  (which  is  difficult  to  gauge  with  newly  synthesized  materials),  the  errors  may  be  as 
high  as  10%  on  average.  These  results  illustrate,  as  previously  pointed  out  by  FlamertonS8  and  others. 
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the  importance  of  carefully  considering  issues  relating  to  baseline  construction  in  the  treatment  of  DSC 


chemical  kinetics  data. 

S6.3.5  Step-by-step  implementation  of  DSC  cure  kinetics  method.  Figure  S6-2  illustrates  the  steps  in 
the  DSC  program,  which  are  as  follows:  A)  ramp  at  5  °C/min  to  120  °C  in  order  to  melt  the  monomer  B) 
isothermal  for  5  minutes  (check  for  baseline  drifts  and  shifts)  C)  jump  to  isothermal  cure  temperature 
(approx.  100  °C  /  min.)  D)  isothermal  cure  for  30  minutes  E)  quench  to  120  °C  (note  that  in  reality  the 
quench,  at  approx.  100  °C  /  min  only  reaches  minimum  temperatures  near  135  °C  before  the  program 
set-point  again  induces  heating)  F)  heat  at  10  °C  /  min.  to  350  °C  (residual  cure).  Following  this  program 
the  DSC  is  cooled  (step  G)  to  120  °C  and  steps  B  through  F  (labeled  B'  through  F')  are  repeated  in  order 
to  generate  the  re-scanned  baseline. 
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Figure  S6-2.  Illustration  of  DSC  isothermal  kinetics  program  (see  text  at  start  of  S6.3.5  for  detailed 
description  of  steps)  showing  temperature  (solid  line,  left  axis)  and  measured  heat  flow  (dashed  line, 
right  axis)  for  the  ESR-255  sample  with  an  isothermal  cure  temperature  of  250  °C. 
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To  prepare  the  data  for  kinetics  analysis,  perform  the  following  steps: 


1.  Subtract  baseline.  Subtract  heat  flow  values  found  in  steps  C,  D,  and  F  from  their  matching  re¬ 
scanned  baseline  values  (C',  D',  and  F',  respectively).  For  step  C,  match  values  on  the  basis  of  heating 
rate,  for  step  D,  match  values  on  the  basis  of  elapsed  time,  and  for  step  F,  match  values  on  the  basis  of 
temperature.  Note  that  the  boundary  between  sections  C  and  D  (as  well  as  C'  and  D')  is  defined  as  the 
point  in  time  when  the  heating  rate  falls  below  0.1  °C  /  min.  Herein  the  resultant  sets  of  heat  flow 
difference  data  are  denoted  AC,  AD,  and  AF. 

2.  Offset  non-isothermal  data  to  correct  for  baseline  shift.  Find  the  minimum  (most  negative)  value  in 
AF  (referred  to  as  AFmin).  Subtract  AFmiri  from  all  values  in  AF.  The  resultant  set  of  values  is  referred  to 
as  AF0ffSet.  Note  that  the  minimum  value  of  AF0ffset  should  be  zero. 

3.  Discard  unreliable  data.  Find  the  latest  point  in  time  for  which  an  associated  value  in  AC  <  0.  Set  all 
values  in  AC  associated  with  a  prior  time  to  zero.  Set  any  values  remaining  in  AC  or  in  AD  that  are  <  0  to 
zero.  Find  the  earliest  point  in  time  for  which  an  associated  value  in  AFoffset  equals  zero.  Set  all  values  in 
AFoffset  associated  with  a  prior  time  to  zero. 

4.  Calibrate  heat  flow  data.  Measure  the  enthalpy  of  melting  of  the  sample,  and  multiply  all  the  newly 
adjusted  values  in  AC,  AD,  and  AFoffset  by  the  ratio  of  the  known  enthalpy  of  melting  for  the  material  to 
the  measured  enthalpy  of  melting.  The  resulting  heat  flows  (as  functions  of  time,  and  normalized  by 
weight)  are  referred  to  as  qc(t),  qD(t),  and  qFoffset(t). 

5.  Integrate  heat  flows.  Integrate  the  heat  flows  qc(t)  and  qD(t),  as  a  function  of  time  (define  the 
reference  time  t'  =  t  -t0,  where  t0  is  the  time  associated  with  the  first  non-zero  value  of  qc)  to  generate 
the  isothermal  enthalpy  of  cure  Hiso(t').  Integrate  the  heat  flow  qF0ffset(t)  over  its  entire  range  to 
compute  the  residual  enthalpy  of  cure  Hres.  For  fully  curable  samples,  define  the  total  enthalpy  of  cure 
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(H0)  as  the  maximum  value  of  Hiso(t')  +  Hres,  or,  for  samples  such  as  ESR-255  that  cannot  be  fully  cured 


without  thermal  degradation,  use  a  previously  determined  value  for  H0  (for  ESR-255,  110  kJ/eq.,  that  is, 
866  J/g). 

6.  Generate  conversion  and  conversion  rate  data  sets.  At  intervals  of  every  0.25  minutes  (or  0.125 
minutes  for  samples  that  reach  50%  conversion  in  under  3  minutes)  starting  at  t'  =  0,  average  the  values 
of  qc(t')  and/or  qD(t'),  divided  by  H0,  over  a  window  of  0.1  minute  centered  on  t'  to  generate  conversion 
rate  (da/dt)  data.  At  the  same  intervals,  divide  (without  averaging  over  a  window),  the  last  available 
(prior  to  t')  value  of  Hiso(t')  by  H0  to  obtain  conversion  (a)  data.  The  result  should  be  a  set  of  values  of 
da/dt  and  a  at  regular  intervals  oft',  referred  to  herein  as  "isothermal  kinetics  data" 

S6.4  Analysis  of  Isothermal  Kinetics  Data  using  Kamal  Model 

S6.4.1  Introduction.  In  previous  workS7  we  have  analyzed  cyanate  ester  isothermal  kinetics  data  using 
the  Kamal  model,S9  namely 

^  =  kx' (1  -  a)n  +  k2'am(  1  -  a)n  (S6-1) 

in  combination  with  a  modified  version  of  Kenny's  graphical  method/10  and  found  an  excellent  and 
reliable  fit  to  the  data.  The  previously  developed  method,  however,  was  limited  to  data  in  which  the 
catalytic  rate  constant  k/  was  much  smaller  than  the  auto-catalytic  rate  constant  k2'.  We  then  found 
parameters  for  a  version  of  the  Kamal  model  in  which  the  exponents  are  assumed  to  remain  constant 
(at  the  average  values  of  m  and  n  obtained  over  all  temperatures  investigated),  namely 

^  =  /cx  (1  —  a)n  +  k2am{l  —  a)n  (S6-2) 

(note  the  differences  in  notation)  and  utilized  it  to  compute  the  activation  energy  of  the  auto-catalytic 
component  for  the  monomers  under  study.  Note  that  because  of  correlations  among  model  parameters 
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and  actual  reaction  rates  (i.e.  a  different  value  of  m  or  n  will  give  a  different  reaction  rate  for  the  same 


value  of  k2'),  meaningful  values  of  activation  energy  require  exponents  that  are  constant  with  respect  to 
temperature. 

In  what  follows,  we  expand  the  previously  developed  method  to  include  the  more  general  case 
in  which  the  auto-catalytic  parameter  is  merely  significant,  that  is,  ki'  <  k2'  or  ki'  ~  k2'.  The  goal  of  the 
analysis,  however,  remains  the  determination  of  significant  comparative  parameters  in  order  to  better 
understand  structure-property  relationships.  Because  catalysis  of  cyanate  esters  depends  on  impurities 
rather  than  on  network  structure,  and  may  involve  many  different  reaction  mechanisms,  the  catalytic 
rate  constant  conveys  information  of  limited  usefulness  for  structure  property  relationships.  However, 
because  cyanate  esters  (when  cured  under  appropriate  conditions)  form  networks  with  a  single,  well- 
defined  network  structure  based  on  cyanurate  linkages,  and  because  it  is  the  catalytic  effect  of  these 
cyanurate  groups  that  is  responsible  for  auto-catalysis,  the  auto-catalytic  rate  constant  provides  a  useful 
means  of  comparing  the  effect  of  monomer  structure  on  cure  chemistry  involving  a  mechanism  that  is 
common  to  all  monomers.  Our  approach  is  thus  designed  to  emphasize  the  extraction  of  meaningful 
data  relating  to  rates  of  auto-catalysis. 

Section  S6.4.2  covers  the  analysis  of  individual  isothermal  kinetics  data  sets  to  find  the  best  fit 
parameters  of  Eq.  S6-1,  while  in  Section  S6.4.3,  the  determination  of  the  rate  constants  in  the  Kamal 
model  with  temperature-independent  exponents  (Eq.  S6-2)  is  described. 

S6.4.2  Kamal  model  parameters  for  individual  isothermal  kinetics  data  curves.  In  the  Kamal  model, 
there  are  four  parameters  to  be  determined.  In  Kenny's  graphical  method,  two  of  these  parameters  are 
determined  via  linear  regression,  while  a  third  (the  catalytic  rate  constant)  is  determined  by  matching  a 
key  feature  of  the  data  (the  initial  reaction  rate)  to  the  model  prediction  (which  depends  only  on  a  single 
parameter),  with  the  fourth  parameter  determined  by  visual  judgment  of  which  value  yields  the 
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straightest  line  for  linear  regression.  Since  Kenny's  graphical  method  was  developed  before  the  advent 


of  powerful  desktop  computers,  the  visual  method,  despite  its  subjectivity,  was  needed  to  enable  the 
analysis  to  be  completed  in  a  reasonable  period  of  time  in  practical  settings.  In  the  previously 
developed  approach,  we  replaced  the  subjective  visual  analysis  with  a  single-parameter  iteration  that  is 
now  easily  carried  out  in  a  few  seconds  on  a  desktop  computer.  The  iteration  relied  on  matching 
another  key  feature  of  the  data,  the  conversion  at  which  the  reaction  rate  is  maximum  (amax),  which 
when  k/  «  k2',  is  used  as  the  basis  of  a  simple  constraint  among  two  model  parameters  that  replaces 
the  "straightest  line"  constraint. 

In  expanding  the  model  to  more  general  situations,  there  are  two  key  complications.  The  first  is 
that  the  parameter  amax  is  no  longer  given  by  a  simple  analytical  formula,  however,  a  more  complex 
procedure  can  still  be  used  in  which  the  model  amax  is  matched  to  the  experimental  value  via  a  second 
level  of  iteration  (the  resulting  procedure  thus  involves  a  two-level,  nested  iteration  using  two  data 
features  to  determine  two  model  parameters,  with  two  others  determined  by  linear  regression).  The 
second  is  that  the  shape  of  the  reaction  rate  curve  as  a  function  of  conversion  may  be  more  complex, 
starting  at  a  non-zero  value  at  zero  conversion,  and  then  either  rising  through  a  single  maximum,  or,  if  m 
~  n  »  1,  falling  to  an  initial  local  minimum  then  rising  through  a  maximum  as  conversion  increases.  In 
the  latter  case,  the  range  of  conversions  in  which  auto-catalysis  dominates  is  between  the  early  local 
minimum  (herein  referred  to  as  amin)  and  amax.  Our  method  is  therefore  designed  to  focus  on  modeling 
reaction  rates  for  conversions  between  amin  and  amax,  (with  amin  =  0  in  many  cases).  The  flow  chart  in 
Figure  S6-3  summarizes  the  method. 

The  outer  iteration  loop  is  primarily  designed  to  match  the  model  value  of  amax  to  the 
experimental  value.  To  accomplish  this  task,  an  intermediate  parameter  amax.k2  is  defined  as  equal  to 
the  ratio  m  /  (m+n).  This  parameter  represents  the  conversion  at  which  the  reaction  rate  due  to  auto- 
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Figure  S6-3.  Flow  chart  illustrating  the  procedure  for  fitting  individual  isothermal  kinetics  data  to  Kamal 
model  parameters  using  a  modified  version  of  Kenny's  graphical  method. 


catalysis  is  maximum.  When  ki'  =  0,  amax_k2  and  amax  coincide  exactly.  The  outer  loop  iteration  works  on 
the  principle  that  if  k/  is  not  too  large,  amax.k2  and  amax  will  be  fairly  close  to  one  another,  and 
approximately  related  by  amax.k2  =  amax  +  Ckl,  where  the  offset  Ckl  is  likely  most  strongly  a  function  of 
kZ/k/.  To  initialize  the  outer  loop,  the  initial  guess  Cki  =  0  (that  is,  amax.k2  =  amax)  is  used.  The  inner  loop 
then  operates,  as  described  later)  by  choosing  a  value  of  m,  using  the  guess  for  amax_k2  to  compute  n, 
then  computing  k/,  k2',  and  the  updated  guess  for  m,  until  the  inner  loop  converges.  After  convergence, 
the  reaction  rate  is  calculated  for  all  values  of  a  with  a  resolution  of  0.001  based  on  the  converged  inner 
loop  values  of  k/,  k2',  m,  and  n.  From  the  set  of  reaction  rate  values,  the  value  of  a  associated  with  the 
largest  computed  reaction  rate,  that  is,  amax.modei,  is  selected.  The  model  and  observed  values  of  amaxare 
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then  compared,  and  a  new  estimate  for  Cki  is  computed  according  to  the  simple  formula,  Cki  =  Cki  +  amax 


-  otmax -model-  The  updated  value  of  Ckl  is  then  used  to  generate  a  new  value  of  amax_k2  and  repeat  the  inner 
loop  based  on  the  updated  constraint  between  m  and  n.  The  process  continues  until  the  absolute  value 
of  Cki  is  less  than  0.0005,  a  criterion  that  was  always  met  within  four  iterations  of  the  outer  loop.  Note 
that,  in  order  to  improve  the  estimation  of  amax,  the  experimental  value  is  refined  by  fitting  a  cubic 
polynomial  to  the  five  data  points  spanning  the  experimentally  observed  maximum  at  the  mid-point  and 
using  the  maximum  of  the  interpolating  polynomial. 

Figure  S6-4  illustrates  the  process  using  the  data  for  ESR-255  cured  at  270  °C.  Figure  S6-4  shows 
the  polynomial  interpolation,  which  is  used  to  refine  the  estimate  of  amaxfrorn  0.549  (with  resolution  of 
about  0.016)  to  0.546  (with  a  resolution  of  about  0.001).  For  the  initial  guess,  amax_k2  =  0.546,  which, 
after  operation  of  the  inner  loop,  generates  the  model  curve,  shown  compared  to  experimental  data  in 
Figure  S6-5.  After  the  first  operation  of  the  inner  loop,  amax_m0dei  =  0.523,  hence  Ckl  is  then  updated  from 
zero  (the  initial  guess)  to  0.023.  This  generates  amax_k2  =  0.569  for  the  next  iteration,  which,  after 
operation  of  the  inner  loop  with  the  updated  constraint  between  m  and  n,  generates  slightly  different 
values  for  k/,  k2',  m,  and  n,  and  amax_modei=  0.543.  Cki  is  then  incremented  from  0.023  to  0.026,  and  the 
process  repeated.  The  outer  loop  converges  at  Ckl  =  0.027,  or  amax_k2  =  0.573.  The  set  of  k/,  k2',  m,  and 
n  for  which  m  /  (m+n)  =  0.573  are  then  used  as  the  "best  fit"  model  parameters. 
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Figure  S6-4.  Isothermal  kinetics  data  near  the  maximum  reaction  rate  for  ESR-255  cured  at  270  °C, 
illustrating  the  use  of  polynomial  interpretation  to  improve  the  experimental  estimate  of  amax. 


o  Exp 
- Model 


Figure  S6-5.  Model  and  experimental  isothermal  kinetics  data  near  the  maximum  reaction  rate  for  ESR- 
255  cured  at  270  °C,  after  the  first  iteration  of  the  outer  loop,  amax  =  0.546,  amax.k2  =  0.546,  and  the 
resulting  amax_modei  =  0.523.  The  updated  estimate  for  Ckl  (initially  zero)  is  thus  0.023,  and  the  updated 
estimate  for  amax.k2  =  0.569.  After  further  iterations,  the  loop  converges  at  amax.k2  =  0.573. 
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For  the  inner  loop,  the  main  consequence  of  having  k/  values  comparable  to  k2'  values  is  the 


effect  on  curve  shape.  As  mentioned  previously,  the  key  conversion  range  for  fitting  data  lies  between 
the  parameters  amin  and  amax,  with  amin  =  0  in  some  cases,  rather  than  always  starting  at  zero 
conversion.  To  reflect  this,  the  procedure  developed  previously  for  finding  k/,  k2',  m,  and  n  is  updated 
so  that,  rather  than  performing  regressions  over  fixed  conversions,  the  conversion  ranges  are  always 
computed  relative  to  amin  and  amax.  The  determination  of  k/  also  becomes  more  complex,  and  changes 
with  each  iteration  of  the  loop  when  amin  does  not  equal  zero.  In  practice,  however,  the  inner  loop  is 
easily  automated. 

To  start  the  inner  loop,  a  value  of  the  parameter  m  is  chosen  as  a  starting  guess  (in  these  cases, 
m  =  1  was  used).  Based  on  the  current  value  of  amax_k2  from  the  outer  loop,  the  relation  amai(.k2  =  m  /  (m 
+  n)  is  used  to  compute  the  value  of  n.  To  compute  k/,  a  linear  regression  of  all  data  points  within  and 
adjacent  to  the  interval  amin  +  0.1  to  amin  +  0.15  is  performed.  The  predicted  value  of  the  regression  at 
amin,  that  is,  (da/dt)min-modei,  is  then  used  to  compute  k/  using  the  assumption  that  the  effect  of  auto¬ 
catalysis  at  the  low  conversions  involved  is  small  and  positive,  and  approximately  large  enough  to  just 
compensate  for  the  slight  underestimate  of  the  experimental,  (da/dt)min  due  to  performing  the 
regression  on  data  with  a  slight  upward  curvature.  In  practice,  this  assumption  works  very  well.  Hence, 
the  auto-catalytic  part  of  the  rate  formula  is  assumed  negligible  and  the  formula  is  thus  k/  =  (da/dt)miri  / 
(1  -  amin)n.  Since  this  formula  does  not  involve  k2',  it  can  be  updated  along  with  m  and  n  at  the  start  of 
each  pass  through  the  inner  loop.  Note  that  when  amin  =  0,  k/  becomes  simply  the  intercept  of  a  linear 
regression  of  the  data  from  a  =  0.10  to  a=  0.15,  thereby  guaranteeing  that  data  from  a  <  0.1,  which  is 
subject  to  the  effect  of  transients,  is  never  used  in  the  computation  of  k/.  Figure  S6-6  illustrates  how  k/ 
is  computed  when  amin  *  0,  showing  how  the  selected  method  more  accurately  estimates  the  scale  of 
the  auto-catalytic  effect  compared  to  alternative  approaches. 
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Figure  S6-6.  Comparison  of  proposed  methods  for  determination  of  k/  based  on  the  resultant  model 
curves  (solid  lines)  compared  to  experiment  (symbols)  for  ESR-255  cured  at  270  °C.  Method  1,  which  is 
described  in  detail  in  the  text,  provides  the  best  agreement  between  the  model  and  experiment  with 
regard  to  the  auto-catalytic  reaction  rate,  as  quantified  by  the  difference  between  maximum  (near  50% 
conversion)  and  minimum  (18.5%  conversion)  reaction  rates.  Since  the  maximum  rate  is  the  same  for  all 
models,  the  value  of  the  reaction  rate  (da/dt)  at  the  minimum  represents  the  discriminating 
characteristic.  The  difference  between  model  and  experiment  was  5%  for  method  1,  24%  for  method  2, 
and  11%  for  method  3. 

Once  m,  n,  and  k/  are  computed,  the  second  part  of  the  inner  loop  computes  k2'  and  updates 
the  estimate  of  m.  To  do  this,  a  linear  regression  based  on  the  following  rearrangement  of  Eq.  S6-1  is 
used: 


da 


In  ,  —  fei'  =  In  k2'  +  rn  In  a 

(1  -a)n  1  z 


(S6-3) 
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A  linear  regression  of  the  left-hand  side  of  Eq.  (S6-3),  which  contains  only  parameters  available  from  the 


earlier  part  of  the  inner  loop,  as  a  function  of  In  a  yields  an  estimate  of  k2'  from  the  intercept  and  a 
second  estimate  of  m  (herein  referred  to  as  mreg).  The  regression  is  carried  out  over  the  conversion 
range  from  amin  +  0.1  to  amax.  To  complete  the  loop,  the  objective  function  mreg  -  m  is  computed,  and 
the  Newton-Raphson  algorithm  is  utilized  (as  part  of  the  Microsoft®  Excel  Goal  Seek  function)  to 
generate  an  updated  estimate  for  m.  The  updated  estimate  for  m  is  then  fed  back  into  the  inner  loop, 
and  the  procedure  repeated  until  the  objective  function  converges  to  within  0.0001.  Once  converged, 
the  parameters  k/,  k2',  m,  and  n  are  passed  back  to  the  outer  loop  to  compute  amax_modei  and  continue  on 
with  the  outer  loop  as  described  earlier.  Figure  S6-7  shows  an  example  of  the  regression  used  in  the 
inner  loop.  The  resultant  model  parameters  for  each  of  the  individual  isothermal  kinetics  data  sets 
collected  as  part  of  this  work  are  tabulated  in  Section  S6.6. 
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Figure  S6-7.  Linear  regression  used  to  determine  In  k2'  and  mreg,  for  ESR-255  cured  at  270  °C.  The  filled 
circles  represent  the  subset  of  experimental  data  points  used  to  compute  the  regression  parameters. 
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S6.4.3  Kamal  model  with  temperature-independent  exponents.  In  order  to  compute 


activation  energies,  it  is  important  to  have  a  rate  constant  that  allows  for  quantitative  comparisons  of 
the  actual  reaction  rate.  Although  such  a  point  might  seem  trivial,  in  fact,  some  care  must  be  exercised 
to  make  sure  that  comparability  is  maintained.  In  analyzing  individual  isothermal  kinetics  data,  for 
instance,  the  value  of  the  constant  k2'  that  best  describes  the  data  (found  through  the  procedure 
described  in  Section  S6.4.2)  changes  based  on  the  value  of  the  reaction  exponents  m  and  n  that  are 
simultaneously  determined  to  best  describe  the  data.  As  a  result,  in  two  data  sets  where  m  and  n 
differ,  the  ratio  of  the  rate  constants  does  not  provide  a  direct  measure  of  the  ratio  of  reaction  rates.  To 
avoid  this  issue  when  calculating  activation  energies  from  a  traditional  Arrhenius  plot,  it  is  necessary  to 
use  temperature-independent  rate  constants,  which  are  herein  denoted  as  m  and  n.  The  values  of 
these  exponents  were  determined  by  averaging  the  values  of  m  and  n  found  in  all  available  sets  of 
isothermal  kinetics  data  for  a  given  compound  within  the  temperature  range  of  210  -  290  °C.  With 
these  two  parameters  fixed,  the  individual  isothermal  kinetics  data  sets  were  then  re-analyzed  to 
determine  a  new  set  of  rate  constants  (herein  denoted  ki  and  k2)  that  were  suitable  for  computation  of 
activation  energies. 

Since  only  two  parameters  needed  to  be  determined,  a  single  linear  regression  based  on  the 
following  form  of  Eq.  (S6-2)  was  utilized: 

da 

T^-K  =  k1  +  k2a™  (S6-4) 

A  plot  of  the  left-hand  side  as  a  function  of  am,  using  the  data  over  a  conversion  range  from  amin  +  0.1 
to  amax  yields  a  slope  of  k2  and  an  intercept  of  ki.  Figure  S6-8  illustrates  this  regression.  In  Figure  S6-9, 
the  resultant  model  isothermal  kinetics  curve  is  compared  to  the  data  as  well  as  the  individually  fitted 
curve  found  by  the  method  in  Section  S6.4.2,  for  the  ESR-255  sample  cured  at  270  °C.  As  expected. 
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having  only  two  adjustable  parameters  reduces  the  quality  of  the  fit  slightly,  however,  the  key  features 


of  the  data,  including,  most  importantly,  the  scale  of  the  auto-catalytic  effect,  are  still  modeled  well. 
The  resultant  data  sets  for  all  of  the  isothermal  kinetics  data  collected  as  part  of  this  work  are  tabulated 
in  Section  S6.6,  along  with  the  resultant  activation  energies  (the  plots  are  shown  in  Figure  6  of  the  main 
manuscript). 


Figure  S6-8.  Regression  used  to  calculate  ka  and  k2  from  experimental  data  on  ESR-255  cured  at  270  °C. 
Note  that  only  the  filled  symbols  are  used  as  experimental  input  data  points  for  the  regression. 


61 


0.0012 


0.001 

_  0.0008 

T— 

0.0006 
"3 

73  0.0004 
0.0002 
0 

0 


o  Experimental 

- Fixed  Exponents 

- Individual  Fit 


- 1 - 1 - 1 - 1 - 1 

0.2  0.4  0.6  0.8  1 

a 


Figure  S6-9.  Comparison  of  experimental  and  model  kinetic  data  (both  individual  curve  fit  and  fit  with 
temperature-independent  exponents)  for  ESR-255  cured  at  270  °C. 


S6.5  Alternative  Analyses  for  Error  Estimates.  Although  the  regressions  used  enable  estimation  of  the 
errors  in  individual  parameters  for  a  given  regression,  these  errors  (which  we  have  published  for  some 
data  sets  previously),511  do  not  reflect  the  propagation  of  error  through  the  complex  process  of 
generating  the  nonlinear  model.  As  with  the  treatment  of  DSC  data,  a  more  informative  approach  to 
estimating  the  errors  is  to  examine  the  variability  of  results  that,  if  error  were  not  present,  would  be 
expected  to  show  no  variation.  For  the  model  exponents  m  and  n,  the  variability  in  the  individual  data 
sets,  which  was  10-20%,  represents  a  reasonable  estimation  of  error.  For  the  rate  constants,  the 
standard  error  in  the  Arrhenius  plot  provides  one  means  of  estimation,  and  appears  to  be  around  5%  as 
shown  in  Figure  6  of  the  main  manuscript. 

An  additional  technique  that  we  have  found  useful  in  the  pastS7  has  been  to  overlay  the 
individual  rate  curves  with  a  vertical  scaling  factor.  This  scaling  factor  represents  the  relative  rate  of 
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reaction  as  a  function  of  temperature  without  the  need  to  invoke  any  particular  model.  An  Arrhenius 


plot  of  the  scale  factor  then  provides  both  an  alternative  estimate  of  the  activation  energy,  which  allows 
for  an  error  estimate  of  this  parameter,  which  the  standard  error  of  the  plot  provides  another  indication 
of  the  relative  error  in  the  rate  constants.  Figures  S6-10  through  S6-12  show  the  overlaid  data  for  the 
three  monomers  studied.  It  is  clear  that  at  both  low  conversions  (at  low  temperatures)  as  well  as  high 
conversions  beyond  amax,  the  curves  do  not  overlay  well.  Within  the  analysis  region,  the  overlay  is  best 
for  ESR-255,  followed  by  monomer  1,  then  monomer  2.  The  Arrhenius  plots  of  the  scale  factors  are 
shown  in  Figure  S6-13.  The  activation  energies  from  the  scale  factors,  like  the  rate  constants,  exhibit  a 
standard  error  of  around  5%.  These  activation  energies,  however,  differ  from  those  found  by  modeling 
by  only  1%  for  ESR-255,  6%  for  monomer  1,  and  11%  for  monomer  2.  Thus  the  level  of  disagreement 
appears  to  correlate  with  the  lack  of  superposition  in  the  relative  rate  curves.  Note  that  these  error 
estimates  include  both  the  error  in  extracting  the  conversion  rates  from  the  raw  DSC  data  (3-5%)  as  well 
as  additional  errors  introduced  by  the  assumption  of  constant  activation  energy.  The  similarity  of  error 
magnitudes  suggests  that  the  latter  assumption  introduces  a  relatively  small  error  into  the  process, 
while  the  differences  in  reported  activation  energy  using  the  two  techniques  suggest  that  the  model- 
based  approach  also  introduces  a  comparatively  small  source  of  error. 
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Figure  S6.10.  Relative  reaction  rate  overlay  for  isothermal  cure  of  monomer  1. 


Figure  S6.ll.  Relative  reaction  rate  overlay  for  isothermal  cure  of  monomer  2. 
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Figure  S6.12.  Relative  reaction  rate  overlay  for  isothermal  cure  of  ESR-255. 


1  /  T  (1  /  K) 


Figure  S6-13.  Arrhenius  plot  based  on  maximum  observed  conversion  rates  (no  model)  in  isothermal 
kinetics  data. 
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S6.6  Tables  of  Kinetic  Model  Parameters 


Table  S6-4 

Individual  isothermal  kinetic  model  parameters  for  monomer  1 


Cure  Temperature 

ki'  (1/s) 

k2'  (1/s) 

m 

n 

210  °C  * 

0.000202 

0.00528 

1.588 

3.535 

230  °C 

0.000472 

0.00271 

1.023 

1.251 

250  °C 

0.000889 

0.00888 

1.162 

1.485 

270  °C 

0.00170 

0.0179 

1.191 

1.625 

290  °C 

0.00224 

0.0291 

1.031 

1.356 

*not  used  in  computation  of  average  exponents  due  to  the  large  discrepancy  with  other  data  sets. 


Table  S6-5 

Individual  isothermal  kinetic  model  parameters  for  monomer  2 


Cure  Temperature 

ki'  (1/s) 

k2'  (1/s) 

m 

n 

210  °C  * 

n/a 

n/a 

n/a 

n/a 

230  °C 

0.000438 

0.00228 

1.326 

1.249 

250  °C 

0.000794 

0.00522 

1.433 

1.206 

270  °C 

0.00125 

0.00836 

1.314 

0.940 

290  °C 

0.00189 

0.0221 

1.234 

0.978 

*  sample  did  not  reach  maximum  cure  rate  during  test,  therefore,  no  parameters  were  calculated 


Table  S6-6 


Individual  isothermal  kinetic  model  parameters  for  ESR-255 


Cure  Temperature 

ki'  (1/s) 

k2'  (1/s) 

m 

n 

210  °C  * 

0.0000496 

0.00124 

2.165 

2.198 

230  °C** 

n/a 

n/a 

n/a 

n/a 

250  °C  *** 

0.000227 

0.0216 

3.333 

2.566 

270  °C 

0.000729 

0.0437 

3.197 

2.382 

290  °C 

0.00115 

0.0390 

2.551 

1.773 

*  sample  did  not  reach 

maximum  cure  rate  during  30  minute  test, 

therefore  no  parameters  were 

calculated  from  the  30  minute  data  set. 

but  data  from  a  24  hour  isothermal  test  was  available  and  used. 

though  not  for  calculation  of  average  exponents  due  to  the  large  discrepancy  with  other  data  sets 

**  sample  did  not  reach  maximum  cure  rate  during  30  minute  test,  therefore  no  parameters  were 

calculated 

***  sample  did  not  reach  maximum  cure  rate  during  30  minute  test,  therefore  no  parameters  were 
calculated  from  the  30  minute  data  set,  but  data  from  a  4  hour  isothermal  test  was  available  and  used 
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Table  S6-7 


Isothermal  kinetic  model  (temperature-independent  exponent)  parameters  for  monomer  1 


Cure  Temperature 

kx (1/s) 

k2  (1/s) 

in 

n 

210  °C 

0.000194 

0.00110 

1.102 

1.429 

230  °C 

0.000494 

0.00319 

1.102 

1.429 

250  °C 

0.000834 

0.00823 

1.102 

1.429 

270  °C 

0.00165 

0.0147 

1.102 

1.429 

290  °C 

0.00256 

0.0315 

1.102 

1.429 

Table  S6-8 

Isothermal  kinetic  model  (temperature-independent  exponent)  parameters  for  monomer  2 


Cure  Temperature 

ki (1/s) 

Mi/s) 

in 

a 

210  °C  * 

n/a 

n/a 

n/a 

n/a 

230  °C 

0.000441 

0.00203 

1.327 

1.093 

250  °C 

0.000739 

0.00458 

1.327 

1.093 

270  °C 

0.00115 

0.0100 

1.327 

1.093 

290  °C 

0.00194 

0.0260 

1.327 

1.093 

*  sample  did  not  reach  maximum  cure  rate  within  the  30  minute  test  time,  therefore  no  parameters 
were  calculated 


Table  S6-9 

Isothermal  kinetic  model  (temperature-independent  exponent)  parameters  for  ESR-255 


Cure  Temperature 

ki  (1/s) 

Mi/s) 

m 

a 

210  °C  * 

0.0000753 

0.00241 

3.027 

2.240 

230  °C** 

n/a 

n/a 

n/a 

n/a 

250  °C  *** 

0.000189 

0.0139 

3.027 

2.240 

270  °C 

0.000668 

0.0351 

3.027 

2.240 

290  °C 

0.00133 

0.0796 

3.027 

2.240 

*  sample  did  not  reach  maximum  cure  rate  during  30  minute  test,  therefore  no  parameters  were 
calculated  from  the  30  minute  data  set,  but  data  from  a  24  hour  isothermal  test  was  available  and  used 

**  sample  did  not  reach  maximum  cure  rate  during  30  minute  test,  therefore  no  parameters  were 
calculated 

***  sample  did  not  reach  maximum  cure  rate  during  30  minute  test,  therefore  no  parameters  were 
calculated  from  the  30  minute  data  set,  but  data  from  a  4  hour  isothermal  test  was  available  and  used 
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S7.  Van't  Hoff  Determinations  of  Purity  by  DSC 


Tables  S7-1  through  S7-3  contain  the  melting  endotherm  characteristics  for  all  DSC  runs  of  monomer  1, 
monomer  2,  and  ESR-255,  respectively.  These  runs  include  the  isothermal  kinetics  investigations  for  all 
monomers  as  well  as  the  runs  used  to  determine  the  diBenedetto  parameters  for  ESR-255.  The 
tabulated  data  include  the  melting  point,  enthalpy  of  melting,  and  van't  Hoff  purity  (as  determined  by 
the  TA  Instruments  Universal  Analysis  software  package).  The  averages  are  also  shown.  Note  that,  as 
mentioned  in  Section  S6,  the  enthalpy  of  melting  was  used  as  an  internal  standard  to  calibrate  the  heat 
flow  data,  because,  as  the  tables  show,  small  weighing  errors  resulted  in  variations  in  the  heat  flow 
values.  The  tables  show  clearly  that  ESR-255  exhibited  the  highest  van't  Hoff  purity,  followed  by 
monomer  2,  and,  at  a  slightly  lower  value,  monomer  1.  The  van't  Hoff  method  is  based  on  melting  point 
depression,  and  therefore  may  not  quantitatively  represent  the  actual  impurity  levels  present. 
Nonetheless,  measurements  of  both  as-received  (99.4%  purity)  commercial  BADCy  as  well  as  BADCy 
with  2%  added  impurity  show  van't  Hoff  purities  near  99.5%  and  98%,  respectively. 


Table  S7-1 

Melting  endotherm  characteristics  of  monomer  1  samples 


Sample  Type 

Melting  Point  (°C) 

AHm  (J/g) 

Van't  Hoff  purity 

Kinetics  -  210  °C  iso 

100.6 

57.3 

96.17% 

Kinetics  -  230  °C  iso 

100.6 

91.0 

93.25% 

Kinetics  -  250  °C  iso 

100.6 

70.2 

95.28% 

Kinetics  -  270  °C  iso 

100.6 

69.8 

95.13% 

Kinetics  -  290  °C  iso 

100.4 

69.3 

95.60% 

Standard  non-isothermal 

101.1 

75.2 

94.90% 

Average* 

100.7  ±0.1 

72.1  ±4.5 

95.06  ±  0.40% 

*Shown  with  the  uncertainty  in  the  estimate  of  the  mean,  not  the  standard  deviation 
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Table  S7-2 


Melting  endotherm  characteristics  of  monomer  2  samples 


Sample  Type 

Melting  Point  (°C) 

AHm  (J/g) 

Van't  Hoff  purity 

Kinetics  -  210  °C  iso 

122.4 

129.7 

96.08% 

Kinetics  -  230  °C  iso 

122.0 

121.0 

96.17% 

Kinetics  -  250  °C  iso 

121.6 

132.5 

96.50% 

Kinetics  -  270  °C  iso 

122.2 

124.2 

96.61% 

Kinetics  -  290  °C  iso 

122.0 

131.1 

94.85% 

Standard  non-isothermal 

123.3 

120.0 

97.38% 

Average* 

122.310.2 

126.412.2 

96.27  1  0.34% 

*Shown  with  the  uncertainty  in  the  estimate  of  the  mean,  not  the  standard  deviation 


Table  S7-3 

Melting  endotherm  characteristics  of  ESR-255 


Sample  Type 

Melting  Point  (°C) 

AHm  (J/g) 

Van't  Hoff  purity 

DiBenedetto  fit  - 150  °C  iso 

116.0 

74.0 

99.51% 

DiBenedetto  fit  -  220  °C  iso 

115.2 

69.9 

99.44% 

DiBenedetto  fit  -  230  °C  iso 

115.3 

71.6 

99.64% 

DiBenedetto  fit  -  240  °C  iso 

115.9 

66.9 

99.29% 

DiBenedetto  fit  -  245  °C  iso 

115.2 

74.3 

99.41% 

DiBenedetto  fit  -  250  °C  iso 

115.3 

61.7 

99.14% 

DiBenedetto  fit  -  255  °C  iso 

115.5 

73.1 

99.34% 

DiBenedetto  fit  -  260  °C  iso 

115.5 

71.4 

99.17% 

DiBenedetto  fit  -  270  °C  iso 

115.5 

70.2 

99.14% 

DiBenedetto  fit  -  290  °C  iso 

115.2 

69.3 

99.41% 

Kinetics  -  210  °C  iso 

115.0 

73.5 

99.61% 

Kinetics  -  230  °C  iso 

114.9 

80.1 

99.65% 

Kinetics  -  250  °C  iso 

115.2 

67.9 

99.54% 

Kinetics  -  270  °C  iso 

115.1 

65.9 

99.46% 

Kinetics  -  290  °C  iso 

114.9 

79.1 

99.44% 

Kinetics  -  210  °C  iso  24  hrs 

115.6 

71.7 

99.38% 

Standard  non-isothermal 

115.5 

71.5 

99.36% 

Monomer  Tg  check 

115.2 

74.9 

99.38% 

Monomer  /  full  cure  Tg  check 

114.9 

76.2 

99.41% 

Average* 

115.310.1 

71.71 1.0 

99.41 1  0.03% 

*Shown  with  the  uncertainty  in  the  estimate  of  the  mean,  not  the  standard  deviation 
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S8.  Data  on  Monomer  1  Purified  at  Small  Scale  by  Flash  Chromatography 


58.1  Introduction.  Since,  for  cyanate  esters,  the  cure  kinetics  depend  strongly  on  the  concentration  of 
catalytic  impurities,  and  since  the  cure  kinetics  may  affect  the  ultimate  conversion  (and  hence  the  glass 
transition  temperature)  achieved  during  cure,  it  is  reasonable  to  examine  what  role,  if  any,  the  known 
differences  in  impurity  levels  play  in  determining  the  cure  kinetics  and  properties  of  the  cyanate  esters 
under  study.  Given  that  the  impurity  profiles  of  cyanate  esters  synthesized  via  unique  multi-step  routes 
likely  differ  in  many  different  ways,  a  comprehensive  investigation  was  judged  beyond  the  scope  of  this 
work,  and  a  simpler  approach  was  sought.  For  the  simplified  approach,  we  took  a  small  quantity  of 
monomer  1,  which  had  the  lowest  van't  Hoff  purity,  and  passed  it  through  a  flash  chromatography  unit 
using  a  mixture  of  ethyl  acetate  and  hexane.  The  main  fraction  collected  exhibited  a  van't  Hoff  purity  of 
99.6%,  which  is  almost  identical  to  that  of  ESR-255,  the  monomer  with  the  highest  measured  van't  Hoff 
purity.  We  could  thus  compare  monomer  1  and  ESR-255  having  equivalent  van't  Hoff  purity  values,  as 
well  as  check  the  effect  of  purification  on  the  properties  of  monomer  1. 

58.2  Results.  Table  S8-1  compares  the  melting  endotherm  and  isothermal  cure  kinetics  parameters  at 
210  °C  of  monomer  1  at  regular  and  high  purity,  along  with  ESR-255,  with  the  latter  two  samples  cured 
using  an  isothermal  period  of  24  hours  and  analyzed  according  to  the  methods  presented  in  Section  S6. 
A  24  hour  period,  rather  than  30  minutes,  was  chosen  in  order  to  allow  for  analysis  over  the  widest 
possible  range  of  conversion,  as  well  as  to  provide  insight  into  the  physical  properties  of  the  monomers, 
which  were  measured  on  samples  cured  at  210  °C  for  24  hours.  The  total  heat  flow  data  from  the 
residual  heating  scans,  which  were  completed  using  modulated  DSC,  are  shown  in  Figure  S8-1.  In  this 
case,  however,  no  features  that  could  clearly  be  assigned  to  a  glass  transition  temperature  were  present 
in  the  reversible  heating  curves,  and  the  irreversible  heating  curves  were  no  different  in  appearance 
than  the  total  heating  curves  shown  in  Figure  S8-1.  In  these  situations,  Kasehagen  et  al.,su  for  example. 
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have  pointed  out  that  the  onset  of  significant  residual  cure  is  a  good  indicator  of  the  "as  cured"  glass 


transition  temperature,  since  the  rate  of  reaction  is  expected  to  jump  significantly  once  de-vitrification 
takes  place. 


Table  S8-1 

DSC  Parameters  of  Monomer  1,  Highly  Purified  Monomer  1,  and  ESR-255 


Parameter 

Monomer  1  (pre-column) 

Monomer  1  (post-column) 

ESR-255  (no  column) 

Tm  (°C) 

100.7* 

102.8 

115.3* 

AHm  (J/g) 

72.1* 

71.2 

71.7* 

van't  Hoff  purity 

95.06%* 

99.61% 

99.41%* 

ki'  (1/s) 

0.000202 

0.000207 

0.0000496 

ki'  (1/s) 

0.00528 

0.000500 

0.00124 

m 

1.588 

0.713 

2.165 

n 

3.535 

0.841 

2.198 

ki (1/s) 

0.000194 

0.000230 

0.0000753 

k2  (1/s) 

0.00110 

0.00110 

0.00241 

m 

1.102 

1.102 

3.027 

n 

1.429 

1.429 

2.240 

*average  of  all  DSC  runs  (see  Tables  S7-1  and  S7-3  for  error  estimates) 


Figure  S8-1.  Modulated  DSC  of  highly  purified  1  and  ESR-255  after  24  hours  at  210  °C. 
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Interestingly,  Table  S8-1  shows  that  the  regular  and  high  purity  versions  of  monomer  1  differ 


only  modestly  in  cure  kinetics.  Note  that  the  rate  constants  are  only  comparable  when  the  exponents 
are  the  same,  thus  the  catalytic  rate  constant  is  only  15%  higher  after  purification  while  the  auto- 
catalytic  rate  constant  is  unchanged.  The  very  large  difference  in  the  individual  fit  exponents,  though, 
does  indicate  that  some  aspects  of  the  cure  reaction  are  affected.  On  the  other  hand,  the  "as  cured" 
glass  transition  temperature  of  the  highly  purified  1  is  little  different  from  the  regular  purity  version 
(both  are  275-280  °C).  The  onset  temperature  of  residual  cure  is  more  difficult  to  discern  for  ESR-255, 
however,  in  this  case  the  temperature  of  245  °C  indicated  by  TMA  data  (see  Section  S5)  seems 
reasonable.  Had  the  impurities  that  were  removed  been  highly  catalytic,  a  significant  drop  in  the 
catalytic  rate  constant  ka  after  further  purification  would  have  been  expected.  In  a  multi-step  synthesis, 
there  are  many  possible  impurities  that  may  have  little  catalytic  activity.  Conversely,  it  has  been  shown 
that  some  impurities  at  the  level  of  tens  of  ppm  can  have  significant  effects  on  cyanate  ester  cure,  and 
these  may  not  have  been  removed  by  flash  chromatography  of  the  sample.  Hence  the  van't  Hoff  purity 
is  not  always  an  indication  of  the  degree  of  catalysis  that  can  be  expected. 

Whatever  the  differences  in  cure  kinetics  that  resulted  from  further  purification  of  1,  the  glass 
transition  temperatures  of  samples  cured  for  24  hours  at  210  °C  were  not  significantly  affected.  Figure 
S8-2  provides  insight  into  why  this  is  the  case.  The  figure  shows  the  estimated  glass  transition 
temperature  (using  the  diBenedetto  equation)  as  a  function  of  cure  time  obtained  from  the  DSC  scans  of 
the  highly  purified  1  and  ESR-255.  The  vast  majority  of  the  cure  for  even  the  highly  pure  1  takes  place  in 
the  first  hour  or  so,  and  the  minor  differences  between  the  regular  and  high  purity  versions  of  1  do  not 
affect  this  time  significantly  (the  times  to  40%  cure  for  instance,  which  can  be  compared  using  the 
standard  30-minute  isothermal  data  for  1  at  regular  purity,  show  a  difference  of  just  2  minutes,  or  about 
8%).  Both  the  regular  and  highly  purified  versions  of  1  thus  vitrify  at  about  the  same  time,  and  have 
similar  amounts  of  time  available  after  vitrification  for  the  Tg  to  slowly  increase. 
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In  contrast,  the  development  of  a  high  glass  transition  temperature  is  significantly  slower  in  ESR- 


255  according  to  Figure  S8-2.  In  the  main  manuscript,  it  was  noted  that  "as  cured"  (in  an  oven,  not  the 
DSC)  samples  of  ESR-255  had  a  lower  glass  transition  temperature  than  1  when  cured  for  24  hours  at 
210  °C,  and  the  estimated  glass  transition  temperature  data  appears  to  support  the  conclusion  that  the 
highly  purified  sample  also  exhibits  an  "as  cured"  glass  transition  temperature  that  is  as  high  or  higher 
compared  to  ESR-255,  in  accordance  with  the  modulated  DSC  data  seen  in  Figure  S8-1.  It  should  be 
noted  that  the  data  series  in  Figure  S8-2  do  not  cover  the  entire  range  of  cure  times  because,  as  cure 
proceeds,  the  reaction  rate  becomes  so  slow  that  it  is  indistinguishable  from  fluctuations  in  the  DSC 
baseline.  Note  also  that  the  data  in  Figure  S8-2  does  not  depend  on  any  particular  kinetic  model,  simply 
the  estimate  of  conversion  from  the  DSC  data  and  the  diBenedetto  parameters  determined  as  described 
in  Section  S4. 


Figure  S8-2.  Estimated  glass  transition  temperature  as  a  function  of  cure  time  at  210  °C  for  highly 
purified  monomer  1  compared  to  ESR-255. 
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